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INTRODUCTIONARY  ARTICLE 


Dear  Colleagues, 

We  would  like  to  introduce  you  to  the  journal  «Materials  Physics  and  Mechanics)),  which  is  launched  this  year  by 
Institute  of  Problems  of  Mechanical  Engineering  (Russian  Academy  of  Sciences,  St.Petersburg,  Russia)  and  Ad¬ 
vanced  Study  Center  Co.Ltd.  (St.Petersburg,  Russia). 

If  any  one  technology  can  be  said  to  shape  the  face  of  today  society,  it  is  materials  technology.  New  materials, 
materials  processes  and  applications  are  of  major  importance  in  boosting  many  high  technological  sectors.  As  dis¬ 
cussions  and  intererst  in  materials  science  and  technology  grow,  the  journal  “Materials  Physics  and  Mechanics” 
will  became  the  place  to  disseminate  knowledge  of  the  physical  and  mechanical  phenomena  and  processes  in 
materials. 

The  journal  specializes  in  the  rapid  publication  of  international  conference  proceedings  in  every  area  of  materials 
physics  and  mechanics.  All  papers  will  be  subject  to  rigourous  peer  review. 

The  editors  of  the  journal  “Materials  Physics  and  Mechanics”  look  forward  to  receiving  your  comments  and 
contributions  to  the  journal. 

Have  a  useful  enjoyable  read! 

Yours  sincerely, 


Vladimir  Bulatov 
Principal  Editor 


and 


Ilya  Ovid’ko 
Principal  Editor 
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PREFACE 


This  volume  is  a  collection  of  selected  original  papers  presented  at  the  International  Workshop  on  Applied  As¬ 
pects  of  Interfaces  Science  (St.Petersburg,  Russia,  September  22-24,  1999).  The  workshop  successfully  provided  a 
critical,  up-to-date  tutorial  review  and  discussion  on  applied  science  of  nanostructured  materials  and  polycrystal¬ 
line  high-Tc  superconductors.  Focuses  were  placed  on  synthesis,  proccessing,  characterization,  properties,  and 
applications  of  nanostructured  materials  and  polycrystalline  superconductors.  A  total  of  33  scientists  from  8 
countries  (Belgium,  France,  Germany,  Russia,  Singapore,  South  Africa,  UK,  USA)  participated  in  the  workshop, 
making  it  a  really  international  event. 

About  equal  number  of  papers  in  this  volume  cover  theoretical  and  experimental  research  in  the  area  of  interface 
science.  All  papers  have  been  peer-reviewed  prior  to  publication. 

Support  of  the  workshop  from  the  Office  of  the  US  Naval  Research,  International  Field  Office,  Europe;  the 
European  Office  of  the  US  Aerospace  Research  and  Development,  UK;  and  the  European  Research  Office  of  the 
US  Army,  UK  is  greatly  appreciated.  Administrative  support  of  the  Institute  of  Problems  of  Mechanical  Engineering 
(Russian  Academy  of  Sciences,  St.Petersburg,  Russia)  is  also  acknowledged. 


Ilya  Ovid’ko 

Chairman,  Organizing  Committee 
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AC  MAGNETIC  PROPERTIES  OF  COMPACTED  FeCo 

NANOCOMPOSITES 


Anit  K.  Giri,  Krishna  M.  Chowdary  and  Sara  A.  Majetich 

Department  of  Physics.  Carnegie  Mellon  University 
Pittsburgh,  PA  15213-3890.  USA 

Received:  September  23,  1999 

Abstract.  Here  we  report  the  AC  magnetic  properties  of  soft  magnetic  nanocomposites  made  from  compaction 
ofFe  Co,,,  nanoparticles.  Following  a  discussion  of  previous  work  on  soft  magnetic  nanocomposites,  the 
sample  preparation  and  experimental  characterization  by  AC  permeametry  are  described.  The  permeability 
is  constant  and  equal  to  the  DC  value  for  low  frequencies,  but  drops  off  sharply  above  a  characteristic 
frequency  f  A  model  is  developed  to  explore  the  relation  between  f  and  material  parameters  including  the 
effective  anisotropy,  exchange  coupled  volume,  temperature,  and  saturation  magnetization. 


1.  INTRODUCTION 

Many  recent  developments  in  soft  magnets,  which  are 
used  in  power  applications,  have  been  based  on 
amorphous  and  nanocrystalline  materials  [1-5],  In 
amorphous  solids  the  magnetocrystalline  anisotropy  K 
is  by  definition  equal  to  zero,  eliminating  the  main 
material-dependent  contribution  to  the  coercivity.  The 
effective  anisotropy  AT  can  also  be  low  in 
nanocrystalline  materials  if  conditions  specified  by  the 
Random  Anisotropy  Model  are  met  [1-3,  6].  The  grain 
size  must  be  smaller  than  the  magnetic  exchange  length 
L  .,  and  the  grains  must  be  exchange  coupled  with  their 
easy  axes  for  magnetization  randomly  oriented.  The 
combination  of  random  orientation  and  averaging  over 
multiple  grains  makes  the  preference  for  magnetiza¬ 
tion  in  a  particular  direction,  and  therefore  Kc/p  very 
small.  Theoretically,  low  K  materials  become  even 
softer  if  the  grain  size  D  is  small  enough.  The  DC 
magnetic  properties  of  nanocomposites  have  been  mod¬ 
eled  in  terms  of  the  Random  Anisotropy  Model  [1-3]. 
The  magnetization  reversal  field  or  coercivity  Hc  is 
predicted  to  be  proportional  to  D 6,  while  the  perme¬ 
ability  | x=B/H  is  predicted  to  be  proportional  to  D  f‘. 
Here  we  extend  this  approach  to  understand  the  AC 
magnetic  properties  of  nanocomposites. 

Previous  work  on  the  AC  magnetic  properties  of 
bulk  magnetic  materials  has  shown  that  the  perme¬ 
ability  can  be  split  into  in-phase  (p’)  and  out-of-phase 


(p”)  components,  just  like  the  permittivity  of  dielec¬ 
tric  materials: 

(1) 

These  components  obey  a  Kramers-Kronig 
relationship  [7],  so  that  by  knowing  the  value  of  one 
for  all  frequencies,  the  other  is  uniquely  determined. 
The  real  part  p’  is  typically  flat  up  to  a  cutoff  frequency, 
and  then  drops  off,  while  the  imaginary  part  p”  shows 
a  peak  at  the  same  frequency.  In  the  best  high  frequency 
ferrite  materials,  the  cutoff  frequency  is  in  the  GHz 
range,  and  is  associated  with  the  precession  frequency 
of  spins  about  the  applied  magnetic  field.  While  high 
frequency  magnetic  materials  are  highly  insulating  to 
reduce  eddy  current  losses,  their  properties  have  not 
been  analyzed  in  terms  of  their  microstructure. 

A  benchmark  for  soft  magnetic  materials  is  the  to¬ 
tal  power  loss,  P  ,  as  the  material  is  magnetized  and 
demagnetized.  It  depends  on  the  hysteresis  loss,  losses 
due  to  the  eddy  currents  created  by  the  changing  mag¬ 
netic  field,  and  anomalous  losses  associated  with  do¬ 
main  wall  motion  [8],  These  contributions  have  differ¬ 
ent  frequency  dependencies: 

P  =  P  +P„  +  P  = 

tot  hvs  eddy  current  oil  out 

(2) 

W.  f  +  Af/p+P  (/)• 

hys  J  Jr  a  noin  v  J  ' 

P  varies  linearly  with  the  frequency / and  the  area 
within  the  hysteresis  curve  Whys,  P ciUy  cumM  is  quadratic 


Corresponding  author:  Sara  A.  Majetich,  e-mail:  sm70+@andrew.cmu.edu 
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in  frequency  and  inversely  proportional  to  the  resistivity' 
p.  A  is  a  constant.  The  anomalous  loss  P  has  no 
fixed  power  law  dependence  on  the  frequency. 

The  AC  magnetic  properties  of  ferrofluids.  which 
are  small  (<  1 0  nm)  magnetic  particles  dispersed  in  a 
liquid,  show  distinctly  different  behavior  from  that  of 
bulk  materials.  There  arc  two  characteristic  relaxation 
times,  one  due  to  Brownian  relaxation  caused  by 
rotation  of  the  particles  in  the  liquid,  and  the  other  due 
to  Neel  relaxation  caused  by  coherent  rotation  of  the 
atomic  spins  within  the  particle  [9,  10],  Both  processes 
occur  at  frequencies  well  below  the  Larmor  frequency 

U 

Nanocomposites  have  potential  for  two  niche  ap¬ 
plications  of  soft  magnetic  materials.  In  low  frequency 
applications,  amorphous  materials  should  continue  to 
dominate,  since  they  have  the  smallest  hysteretic  losses. 
Without  significant  eddy  currents,  only  the  small  per¬ 
centage  of  atoms  required  for  glass  forming  do  not 
contribute  to  the  magnetization.  However,  at  higher 
temperatures  amorphous  materials  crystallize.  In  the 
current  nanocrystalline  materials  [1-5],  grain  growth 
is  limited  by  precipitation  of  additional  phases,  and  a 
larger  fraction  of  nonmagnetic  atoms  is  needed.  At 
higher  frequencies,  nanocomposites  have  potential 
advantages  because  they  can  have  higher  resistivity  to 
reduce  eddy  current  losses.  While  ferrites  are  currently 
used  for  high  frequency  applications  because  of  their 
high  resistivities,  they  have  much  lower  magnetiza¬ 
tions  than  iron-based  alloys. 

Iron-cobalt  alloys  are  of  interest  for  high  tempera¬ 
ture  soft  magnets  because  they  have  a  high  Curie  tem¬ 
perature,  high  magnetization,  and  low  coercivity.  The 
value  of  the  Curie  temperature  depends  on  the  relative 
Co  abundance.  Iron-cobalt  alloys  have  the  largest  mag¬ 
netic  moment  per  atom,  and  therefore  large  saturation 
magnetizations  [11].  They  have  a  low  coercivity  due  to 
the  low  magnetocrystalline  anisotropy  K  of  the  body- 
centered  cubic  structure.  There  is  a  zero  crossing  of 
the  anisotropy  near  the  equiatomic  composition  [12]. 

We  investigate  nanolaminates  formed  by 
compaction  of  FeCo  nanoparticles  with  a  very  thin 
coating  of  a  protective  carbon  or  oxide.  By  choosing  a 
particle  or  grain  diameter  well  below  the  maximum 
monodomain  size,  the  coercivity  and  hysteretic  power 
losses  can  be  minimized  [1—3].  The  coatings  act  as  a 
barrier  to  eddy  currents,  reducing  the  eddy  current 
power  losses  at  high  frequency.  The  goals  are  to 
minimize  the  portion  of  nonmagnetic  atoms  while 
retaining  a  stable  nanocrystalllne  structure  at  high  tem¬ 
perature,  and  to  determine  the  degree  of  coating  which 
retains  significant  exchange  coupling  for  minimum  K 
while  creating  significant  banders  to  eddy  currents. 


2.  EXPERIMENTAL 

Particle  Synthesis.  Here  the  alloy  nanoparticles  were 
synthesized  by  the  polyol  method  [  1 3],  a  chemical  route 
which  leads  to  highly  monodisperse  metal  particles. 
Iron  chloride  tetrahydrate  and  sodium  hydroxide  were 
dissolved  in  ethylene  glycol  and  heated  to  1 10  °C  while 
stirring.  A  second  solution  of  cobalt  hydroxide  in  eth¬ 
ylene  glycol  at  the  same  temperature  was  added,  and 
the  mixture  was  heated  to  195  °C,  where  a  precipitate 
fonned.  The  total  metal  ion  concentration  was  0.2  M, 
and  the  hydroxide  concentration  was  2.0  M.  Water  and 
other  reaction  products  were  distilled  off,  and  the  sol¬ 
vent  was  refluxed  for  one  hour  at  the  maximum  tem¬ 
perature.  The  precipitate  was  thoroughly  washed  with 
methanol,  and  dried. 

While  the  particle  size  can  be  varied  by  changing 
the  metal  salt  concentration  and  reaction  time,  the  av¬ 
erage  grain  diameter  was  roughly  constant  and  approxi¬ 
mately  20  nm.  The  composition  was  determined  from 
electron  energy  loss  spectroscopy  [14],  The  average 
grain  size  was  obtained  from  Scherrer  analysis  of  the 
x-ray  diffraction  peaks,  and  the  particle  size  was  found 
from  transmission  electron  microscopy.  The  exchange 
length  in  magnetically  soft  materials  like  FeCo  is  large, 
and  here  is  comparable  to  tbe  grain  size.  Adjacent  grains 
can  then  be  exchange  coupled  to  each  other,  leading  to 
a  further  reduction  in  the  anisotropy  and  coercivity,  as 
explained  by  the  random  anisotropy  model  [1-3,  6], 
Compaction.  A  nanocomposite  has  nanoparticle-sized 
grains,  which  may  be  dispersed  in  a  matrix  that  medi¬ 
ates  their  coupling.  The  coupling  strength  depends  on 
the  magnetic  behavior  of  the  matrix  and  on  the  degree 
of  coherence  at  the  interface  between  the  nanoparticles 
and  the  matrix.  Exchange  interactions  may  couple 
nanoparticles  together.  We  worked  with  Dr.  S. 
Sudarshan  and  Mr.  Sang  Yoo  of  Materials  Modifica¬ 
tion,  Inc.,  on  plasma  pressure  compaction,  which  has 
the  advantage  of  being  able  to  reduce  or  remove  the 
thin  coating  surrounding  FeCo  nanoparticles  under 
ambient  conditions.  Powders  were  compacted  by  a 
combination  of  pressure  and  an  electric  arc  in  a  two 
stage  process.  First  a  voltage  pulse  established  a 
continuous  current  path  across  the  sample  in  a  graph¬ 
ite  mold  under  a  pressure  of  1 0  MPa.  This  generated  a 
plasma  in  the  voids,  and  partly  removed  surface 
impurities  such  as  oxides.  In  the  second  stage  a 
continuous  current  path  led  to  resistive  sintering.  The 
pressure  was  increased  to  70  MPa  to  deform  and 
compact  the  particles.  The  entire  process  required  only 
a  few  minutes,  thereby  limiting  grain  growth.  The 
resistive  heating  concentrated  in  the  “necks”  forming 
as  the  particles  sinter,  making  them  deform  more 
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Fig.  1.  a.)  Comparison  of  Fe^Co^,  nanoparticles  made  by  the  polyol  method  (open  circles),  and  a  compact  made  from  them 
(filled  circles);  b.)  An  expanded  region  shows  the  coercivity  of  the  pellet  is  on  the  order  of  3  Oe,  much  less  than  that  of  the 
precursor  particles,  which  is  approximately  130  Oc. 


readily.  Compacts  with  over  90%  theoretical  density 
were  made  from  the  nanoparticle  powders. 

Magnetic  Properties.  There  was  a  substantial  amount 
of  magnetic  characterization,  including  SQUID  magne- 
tometry  (with  a  Quantum  Design  MPMS)  to  determine 
the  coercivity,  and  AC  permeametry  (with  a  Walker  Sci¬ 
entific  AC  Hysteresisgraph)  on  toroidal  compacts  to 
measure  the  frequency  dependence  of  the  permeability 
m  and  the  power  loss.  For  the  AC  measurements,  the 
toroid  outer  diameter  was  2.54  cm,  and  it  had  a  square 
cross  section  with  dimensions  3.0x4. 1  mm.  There  were 
154  primary  turns  and  11  turns  in  the  pickup  coil.  The 
amplitude  of  the  AC  magnetic  field  applied  was  125  Oe, 
and  the  frequency  ranged  from  100  to  100,000  Hz.  The 
temperature  of  the  toroid  was  controlled  by  placing  it  in 
thermal  contact  with  an  ice  bath  (0°C),  or  with  it  inserted 
in  a  furnace  (100  °C  -  500  °C).  At  elevated  tempera¬ 
tures  the  coils  were  sheathed  with  fiberglass  sleeving 
(Omega  Engineering,  Inc.). 

3.  EXPERIMENTAL  RESULTS 

Fig.  1  shows  a  comparison  of  the  DC  hysteresis  loops 
for  the  powder  and  compacted  pellet  of  a  typical  sample. 
Compaction  significantly  reduces  the  coercivity,  indi¬ 
cating  that  the  exchange  coupling  is  increased,  as  pre¬ 
dicted  by  the  random  anisotropy  model. 

The  AC  hysteresis  loops  were  measured  for  toroids 
cut  from  the  same  compressed  pellets,  for  a  range  of 
frequencies.  With  AC  excitation,  the  pickup  loop 
measures  the  magnetic  induction  B~  4kM  \  H.  Fig.  2 
illustrates  the  characteristic  behavior  as  a  function  of 
frequency.  At  low  frequencies  (Fig,  2a)  the  applied  field 
H  was  sufficient  to  saturate  the  sample,  reaching  the 
single-valued  region  of  B{H).  As  the  frequency  was 
raised,  the  width  of  the  hysteresis  loop  and  therefore 
the  coercivity  increased,  but  the  maximum  induction 
remained  approximately  the  same,  up  to  a  characteristic 


frequency  (Fig.  2b).  Approaching  this  frequency,  the 
sample  was  no  longer  saturated  at  the  maximum  ap¬ 
plied  field,  and  the  maximum  magnetic  induction 
dropped  (Fig.  2c). 

The  permeability  as  a  function  of  the  maximum 
applied  field  is  illustrated  in  Fig.  3,  for  several 
representative  frequencies.  Defining  the  permeability 
as  the  maximum  value  of  B/H,  it  is  expected  to  vary 
depending  on  the  magnitude  of  H  .  If  H  is  small 
and  the  sample  is  not  saturated,  B  increases  approxi¬ 
mately  linearly  with  H,  and  p  is  large.  If  Hmax\ s  large 
and  the  sample  is  never  saturated,  B  increases  more 
slowly  with  increasing  //and  p  is  reduced.  In  order  to 
understand  the  fundamental  physics  of  the  AC  excita¬ 
tion  of  the  nanocomposites,  a  value  of  Hmax  =  125  Oe 
was  selected  for  further  measurements  so  that  the 
samples  were  saturated  at  low  frequency. 

Several  important  quantities  for  soft  magnetic  ma¬ 
terials  were  measured  as  a  function  of  frequency.  The 
real  (in-phase)  and  imaginary  (out-of-phase)  parts  of 
the  permeability  \.v-B/H  were  also  measured  as  a  func¬ 
tion  of  frequency  (Fig.  4).  In  the  low  frequency  range, 
the  phase  lag  between  B  and  H  is  small,  and  p”  is  neg¬ 
ligible.  p’  is  roughly  constant  up  to  the  characteristic 
frequency,  where  it  drops  off  sharply,  p”  peaks  as  p’ 
drops  off,  and  the  critical  frequency  fc  is  determined 
from  the  position  of  the  peak  in  p”.  The  characteristic 
frequency  of  the  permeability  varied  with  temperature 
(Fig.  5).  While  the  low  frequency  magnitude  of  m 
dropped  off  slightly  with  increasing  temperature,  due 
to  a  reduction  in  the  magnetization,  fc  rose  substan¬ 
tially,  from  1400  Hz  at  0  °C  to  6000  Hz  at  500  °C.  As 
noted  in  the  hysteresis  loops,  the  coercivity  rose  with 
increasing  frequency  (Fig.  6).  It  peaked  at  a  frequency 
approximately  2-3  times  f  ,  and  decreased  somewhat 
at  high  frequencies.  The  overall  power  loss  per  cycle  is 
shown  in  Fig.  7.  At  low  frequencies  it  rises  linearly,  as 
expected  from  Eq.  2  if  hysteretic  losses  are  dominant. 
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At  frequencies  above  f  a  weaker  dependence  is  ob¬ 
served,  closer  to  f05.  As  the  temperature  is  increased 
the  power  loss  per  cycle  decreases  somewhat,  due 
mainly  to  the  reduction  in  the  coercivity. 

4.  MODEL  OF  AC  MAGNETIC 
PROPERTIES 

Our  model  of  the  AC  properties  of  nanocomposites  is  based 
on  a  monodomain  magnetic  particle  in  an  oscillating 
magnetic  field.  For  a  uniaxial  particle  with  the  applied 
field  parallel  to  the  easy  axis,  the  energy  E  is  given  by 

E  =  AT  sin2  9-  M //cos (9,  (3) 

where  K  is  the  magnetocrystalline  anisotropy,  V  is  the 
particle  volume,  0  is  the  angle  between  the  magneti¬ 
zation  direction  and  the  applied  field  H,  and  M  is  the 
saturation  magnetization.  Stoner-Wohlfarth  theory  [  1 5- 
17]  has  shown  that  energy  minima  occur  at  0  =  0°  or 
180°,  and  a  maximum  exists  where  cos0  =-HM/2K. 


In  order  to  reverse  the  magnetization  direction,  the 
energy  barrier  must  be  overcome.  If  this  does  not  oc¬ 
cur  within  the  measurement  time,  then  hysteresis  is 
observed.  The  rate  of  going  over  the  barrier  is  given 
by 


r 


=  T,  exP  - 


(4) 


where  r~'  =f  is  the  Larmor  frequency,  A E  is  the  mag¬ 
nitude  of  the  energy  barrier,  k  is  the  Boltzmann 
constant,  and  T  is  the  temperature. 

If  the  magnetic  field  oscillates  in  time, 

H(t)  =  Hn  cos(2 nft\  (5) 

then  the  positions  of  the  energy  minima  and  maxima 
will  also  change  with  time,  as  shown  in  Fig.  8.  Early 
in  the  cycle  state  1  is  the  global  minimum  and  state  2 
is  a  local  minimum,  but  later  the  roles  are  reversed. 
The  height  of  the  energy  barrier  which  must  be  over- 


Fig.  3.  Overall  permeability  for  Fc]()Co(J||,  as  a  function  of 
H  ,  at  0  °C,  for  500,  900,  and  1400  Hz.  The  peak 
permeability  at  22  Oc  reflects  where  the  slope  of  the  B(H) 
starts  to  decrease. 


Fig.  4.  The  real  (p’)  and  imaginary  (|i”)  parts  of  the 
permeability  for  Fe|()Co,J(),  at  0  °C,  as  a  function  of  frequency. 
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/(Hz)  /(Hz) 

Fig.  5.  Overall  permeability  as  a  function  of  frequency  for  Fig.  6.  Coercivity  as  a  function  of  frequency  for  Fe|0Co90,  at 
Fe10Co90,  at  0  °C  and  500  °C.  0  °C  and  500  ”C. 


come  for  the  particle  to  switch  from  the  local  to  the 
global  minimum  also  oscillates  in  time  and  is  greatest 
at  zero  applied  field. 

Given  an  ensemble  of  particles  and  monitoring  the 
population  N  in  state  1  and  N2  in  state  2  as  a  function 
of  time,  the  magnetization  can  be  determined: 

M(t)  =  =  M,(2Af(0-l)-  (6) 

Here  we  assume  that  population  is  conserved,  and 
that  at  all  times 

Nt  +  N\  =  1.  (7) 

The  population  in  state  1  as  a  function  of  time  de¬ 
pends  on  the  rate  ,  at  which  particles  get  over  the 
barrier  from  state  2,  and  on  the  rate  r'  that  they  leave 
state  1  in  the  forward  direction  over  the  barrier: 


dNt  _  dN  2 
dt  dt 


the  differential  equation  becomes 


-  =  (1- At,) exp 


-KV  HM , 

-  1  + - 1 

kT  2  K 


/(Hz) 


Fig.  7.  Power  Loss  per  unit  volume  as  a  function  of  frequency 
for  FeinCoQ0,  at  0  °C. 


(11) 


In  an  oscillating  magnetic  field  7/(7),  Eq.  1 1  must  in 
general  be  solved  numerically. 

This  problem  has  many  similarities  to  that  of  an 
electric  dipole  in  an  AC  electric  field,  which  was 
originally  solved  by  Peter  Debye  [18,19],  with  two  main 
differences.  In  the  electric  dipole  model,  there  was  no 
equivalent  to  magnetocrystalline  anisotropy,  so  the  so¬ 
lution  corresponds  to  the  superparamagnetic  limit. 
Also,  in  the  electric  dipole  model,  the  factor  equiva¬ 
lent  to  x  =  M  VH/kT  is  much  less  than  one,  and  when 
ex  is  approximated  by  1  +x,  the  equation  corresponding 
to  Eq.  1 1  can  be  solved  analytically. 
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Fig.  8.  Energy  levels  as  a  function  of  time,  for  different  parts  of  the  cycle  of  H(t ).  Here  £(  is  the  energy  of  a  particle  in  state  1, 
E2  is  the  energy  of  a  particle  in  state  2,  and  A E  is  the  magnitude  of  the  barrier  between  the  local  minimum  and  the  maximum 
energy  states. 


Equation  1 1  can  be  solved  analytically  in  the  limit 
of  large  K,  where  A Ejwd  ~  &Ekkwd-  With  the  substitu¬ 
tion 


KV 


esis  curve  during  each  half-cycle.  In  applications  of 
soft  magnetic  materials,  this  is  undesirable  since  the 
permeability  p  =B/H  is  reduced.  Therefore  the  model 
can  be  used  to  understand  the  AC  magnetic  behavior 
in  situations  of  greatest  technological  importance. 


Eq.  1 1  then  simplifies  to 

dN.  i  , 

-^-*r„V(l-2fV,).  (13) 

dt 

The  solution  to  Eq.  13  is  given  by 

W,(0  =  ^[l  +  exp(-2r>)].  (14) 

There  will  be  a  characteristic  or  resonant  frequency  /  , 
where 


-i  -y  f  -i 

T  a  1  ~ 

/;  =  ^ 


n 


r_^ 
V  * 


exp. 


-KV 

kT 


(15) 


where  the  population  oscillates  in  resonance  with  the 
driving  field  H(t).  Above  this  frequency,  the  popula¬ 
tion  shifts  cannot  keep  up  with  the  rapidly  changing 
driving  field. 

The  model  breaks  down  for  cases  where  HMJ2K  > 
1,  since  this  requires  that  |cos0  =-HMJ2K\  be  greater 
than  one.  To  investigate  the  small  K  limit,  which  is  of 
critical  importance  for  soft  magnetic  materials,  we 
reduced  the  amplitude  of  the  magnetic  field  H{)  where 
necessary.  In  experiments  it  is  possible  to  have  large 
amplitudes  of  the  AC  magnetic  field  when  studying 
materials  with  very  low  anisotropy.  Here  the  sample  is 
being  driven  past  the  high  slope  portion  of  the  hyster¬ 


5.  MODELING  RESULTS  AND 
COMPARISON  WITH  EXPERIMENT 

The  simulation  results  show  the  same  characteristics 
as  the  experimental  data,  and  provide  insight  about 
the  mechanism  responsible  for  the  frequency  depen¬ 
dent  behavior. 

Fig.  9  illustrates  the  connection  between  the  time- 
dependence  of  M  and  the  B(H)  curves  for  driving  fields 
of  different  frequencies  H  H cos  (2n  ft) .  Ifthe  frequency 
is  low  enough  (Fig.  9a)  the  sample  has  time  to  equili¬ 
brate  with  the  slowly  varying  field  and 
superparamagnetism  is  observed.  There  is  a  negligible 
phase  lag  0  between  M(l)  and  H(t).  Because  of  this  M(t) 
reaches  its  maximum  value  well  below  the  maximum 
amplitude  of  the  driving  field  Hf).  At  somewhat  higher 
frequencies  M(t)  starts  to  lag  behind  H(t)  significantly, 
leading  to  hysteresis  (Fig.  9b).  However  M(t)  still 
reaches  its  maximum  value  before  H(t)  starts  to  de¬ 
crease.  At  the  critical  frequency  f  this  is  no  longer 
true  (Fig.  9c),  and  the  shape  of  M(t)  and  the  hysteresis 
loop  change  noticeably.  B(H)  no  longer  contains  single 
valued  regions,  and  the  maximum  value  of  B  starts  to 
decrease  relative  to  its  value  at  lower  frequencies.  At 
the  highest  frequencies  M(t)  and  H(t)  are  almost 
completely  out  of  phase  and  the  maximum  values  of 
M(t)  and  B(H)  are  very  small  (Fig.  9d).  The  model  simu¬ 
lates  characteristics  of  the  experimental  results  shown 
in  Fig.  2,  and  shows  that  the  shape  and  magnitude  of 
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Fig.  9.  M( t)  and  FI(t ),  and  the  corresponding  B(H )  curves  for  a  sphere  of  exchange  coupled  volume  V  and  effective  anisotropy 
K,  a.)  for  a  superparamagnetic  particle,  b.)  for  a  monodomain  but  not  superparamagnetic  particle  at  low  frequency  /,  c.)  at  med. 
/,  and  d.)  at  high  /.  Note  that  in  a.)  and  b.)  B(H)  appears  to  saturate,  but  this  is  because  4v.  M  is  so  much  larger  than  H.  These 
simulations  used  the  following  parameters:  Ms  =  1000  cmu/cc,  H0  —  125  Oe,  K  =  1.95'105  erg/cc,  T=273  K,  and  V  = 
1.47T0'18  cc.  The  critical  frequency  /(  was  1400  Hz. 

B(H)  depend  on  whether  the  magnetization  of  the  ex-  in  the  high  frequency  limit.  This  is  in  contrast  to  the 

change  coupled  volume  within  the  sample  can  equili-  experimental  results  shown  in  Fig.  5,  where  p  decays 

brate  with  the  changing  external  magnetic  field.  more  slowly.  Such  a  difference  would  arise  if  the  sample 

Figure  10  shows  the  corresponding  trend  in  the  contained  a  distribution  of  characteristic  frequencies 

permeability  p  as  a  function  of  frequency.  At  low  frc-  fc,  perhaps  related  to  a  distribution  of  grain  sizes, 
quencies  it  is  constant  and  equal  to  the  DC  value,  but  Fig.  1 1  reveals  the  behavior  of  the  simulated 

as  soon  as  M(t)  is  unable  to  reach  its  maximum  value,  coercivity  as  a  function  of  frequency.  At  low  ffequen- 

B  and  therefore  m  begin  to  drop,  approaching  zero  cies  it  follows  Sharrock’s  Law  [20]: 
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Freguency  (Hz) 


Fig.  10.  The  permeability  (.1  —B  „  /H  j  as  a  function  of 
frequency,  for  the  same  parameters  as  in  Fig.  9. 


Fig.  11.  The  coercivity  H  as  a  function  of  frequency,  for  the 
same  parameters  as  in  Fig.  9. 
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frequency  /)..  The  dependence  of f  on  the  anisotropy, 
volume,  and  temperature  is  summarized  in  Fig.  13. 
From  the  simulated  data,  an  empirical  functional  form 
for  fJK,  V,T)  was  obtained: 


f,  =  10"  T~a'  exp 


-(5.88-1  O'5  K I  erg) 


KV 


-  (4.33  •  1  0“  cc  /  erg)K  +  (4.1  M  02"  K  /  cc )  (1 7) 


+  18.98 


rising  with  the  log  of  the  frequency.  The  peak  in  H 
occurs  between  two  and  three  times  f .  and  then  H 
begins  to  decrease.  This  occurs  when  the  magnetiza¬ 
tion  due  to  the  driving  field  is  low  enough,  and  the 
value  of  H(l)  needed  to  make  M(t)  equal  to  zero  is 
reduced.  The  experimental  coercivity  shown  in  Fig.  6 
also  rises  and  peaks  at  roughly  the  same  frequency,  but 
the  rise  is  not  proportional  to  In  f.  While  there  may  be 
a  region  where  //.  obeys  this  functional  dependence, 
there  appears  to  be  a  low  frequency  contribution  to  the 
coercivity  from  a  factor  not  included  in  the  model.  This 
could  arise  from  particles  not  exchange  coupled  to  their 
neighbors,  possibly  from  a  thicker  coating. 

Fig.  12  shows  the  hysteretic  power  loss  P  found 
from  the  are  of  the  hysteresis  loop  B(H),  as  a  function 
of  frequency.  This  rises  linearly  at  low  frequencies,  as 
expected  from  Eq.  2,  but  then  it  levels  off  at  f  At  very 
high  frequency  it  should  decrease  when  both  Mmax  and 
//.  are  decreasing. 

In  comparison,  Fig.  7  shows  the  total  experimental 
power  loss.  At  frequencies  below/).,  hysteretic  power 
losses  dominate.  This  is  expected  for  frequencies  be¬ 
low  10,000  Hz. 

The  qualitative  behavior  illustrated  in  Fig.  9-12  is 
similar  for  a  wide  range  of  input  parameters;  only  the 
magnitudes  of  p,  //.  and  Ph,s  vary,  along  with  the  critical 


Note  that  this  formula  is  considerably  more  complex 
than  Eq.  15,  which  was  valid  only  in  the  limit  of  large 
K.  For  soft  magnetic  materials  K  is  small,  and 
numerical  simulations  are  the  only  was  to  predict  the 
cutoff  frequency. 

For  the  experimental  grain  size,  r  =10  nm,  and  an 
effective  anisotropy  equal  to  1.5-105  erg/cc  would  be 
needed  to  fit  the  observed  room  temperature  value  of 
f  .  The  is  a  much  larger  value  of  the  anisotropy  esti¬ 
mated  from  the  experimental  coercivity  and  Sharrock’s 
Law.  Presumably  the  physically  relevant  volume  V  is 
not  that  of  a  single  grain,  but  of  a  small  number  of 
coupled  grains. 


Freguency  (Hz) 

Fig.  12.  The  hysteretic  power  loss  as  a  function  of 
frequency,  for  the  same  parameters  as  in  Fig.  9. 


Cutoff  Freguency  (Hz)  Cutoff  ^eguency  (Hz)  ^  Cutoff  Freguency  (Hz) 
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Fig.  13  a.)  fr  as  a  function  of  the  anisotropy  for  spheres  of 
different  radii.  Here  V=(47t/3 fr3-,  b.)  /  as  a  function  of  the 
volume;  c.)  /  as  a  function  of  the  inverse  temperature. 


6.  CONCLUSIONS 

The  AC  magnetic  properties  of  soft  magnetic 
nanocomposites  made  from  compaction  of  Fe|f|Co9(1 
nanoparticles  can  be  understood  in  tenns  of  an  effective 
anisotropy  and  a  related  exchange  coupled  volume.  The 
permeability  is  constant  and  equal  to  the  DC  value  for 
low  frequencies,  but  drops  off  above  a  characteristic 
frequency  /).  For  bulk  materials  such  as  ferrites  the 
cutoff  frequency  can  be  extremely  high,  and  the  be¬ 
havior  is  understood  in  terms  of  the  Landau-Lifshitz- 
Gilbert  equation  [21-23],  However  for  compacted 
nanocomposites  the  cutoff  frequency  is  much  lower, 


and  depends  on  the  sample  microstructure.  The  value 
of  this  frequency  can  be  increased  by  raising  the  tem¬ 
perature,  or  by  reducing  the  anisotropy  and  grain  size. 
Experiments  are  in  progress  to  prepare  samples  to  fur¬ 
ther  test  this  model. 
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Abstract.  Nanostructured  Ni  films  were  deposited  on  Cu  by  reducing  nickel  acetate  tetrahydrate  in  refluxing 
ethylene  glycol  at  about  1 94  °C.  The  effects  of  deposition  time  ( t )  and  Pt  as  nucleating  aid  were  investigated 
in  this  polyol  process.  The  solution  chemistry  varied  with  t.  At  longer  t,  corrosive  attack  of  the  solvent  on 
deposited  films  led  to  dissolution  and  re-deposition  of  Ni.  Saturation  magnetization  increased  with  t  and 
decreased  with  increasing  Pt  precursors.  The  deposited  Ni  films  were  magnetically  saturated.  They  showed 
in-plane  magnetization  anisotropy  and  perpendicular  cocrcivity  anisotropy.  The  use  of  Pt  caused  rapid 
powder  precipitation  but  did  not  enhance  film  deposition.  The  Pt-catalyzed  Ni  films  showed  a  (200)  texture 
and  a  higher  squareness  ratio. 


1.  INTRODUCTION 

Nanostructured  films  have  attracted  much  attention 
because  of  their  unique  properties  that  arise  from  ef¬ 
fects  of  size  reduction  and  large  amount  of  interfaces 
[1  ].  These  films  can  be  prepared  by  many  methods,  for 
example,  sputtering,  laser  ablation,  electrodeposition 
and  electroless  deposition.  For  uniform  metallization 
of  non-planar  and/or  hidden  surfaces,  solution 
chemistry  is  better  suited  than  the  physical  vapor 
methods,  since  it  is  not  limited  by  the  line-of-sight  depo¬ 
sition. 

An  electroless  polyol  method,  which  has  been 
mainly  used  to  synthesize  micron,  submicron  and 
nanostructured  powders  [1-8],  has  been  recently  ap¬ 
plied  for  deposition  of  nanostructured  metal  films,  for 
example,  Cu  film  on  AIN  substrate  [9]  and  NixCo|00x 
film  on  Cu  substrate  [101.  The  Ni  Co...  films  exhib- 
ited  good  magnetic  properties  with  perpendicular 
coercivity  anisotropy  [10].  The  structure  of  these  films 
have  been  recently  investigated  using  anomalous  x- 
ray  scattering  to  determine  if  a  solid  solution  or  phase- 
separated  composite  was  formed  and  the  element-se¬ 
lected  long  range  order  [11].  To  date,  the  film  deposi¬ 
tion  chemistry  and  mechanisms  of  this  polyol  process 
are  not  yet  well  understood.  Unlike  physical  vapor  depo¬ 
sition  where  the  flux  of  atoms  can  be  maintained 
constant,  the  solution  chemistry  (thus  deposition 


chemistry)  may  vary  as  a  function  of  deposition 
(reaction)  time. 

In  this  paper,  we  report  our  investigation  of  the 
magnetic  properties  of  polyol-deposited  Ni  films  and 
their  dependence  on  the  deposition  time  (/)•  The  ef¬ 
fects  of  increasing  t  on  crystallite  size,  film  thickness, 
microstructures  and  magnetic  properties  were  studied. 
The  influence  of  Pt  as  nucleating  aid,  which  has  been 
used  in  polyol  synthesis  of  powders  [12],  on  Ni  film 
deposition  was  also  studied. 

2.  EXPERIMENTAL 

Polycrystalline,  coarse-grained  copper  with  (200)  tex¬ 
ture  was  used  as  substrates  without  any  pre-deposition 
surface  treatment  or  polishing.  0.1  M  nickel  acetate 
tetrahydrate  was  suspended  in  200  ml  ethylene  glycol. 
Hexachloroplatinic  acid  in  small  quantity  was  added 
in  reactions  where  the  effects  of  nucleating  aid  were 
investigated.  The  ratio  of  concentration  of  Pt  to  Ni 
([Pt]/[Ni])  was  estimated  from  the  starting  precursors. 
The  substrate  was  vertically  suspended  in  the  mixture, 
which  was  brought  to  reflux  at  about  194  °C.  Deposi¬ 
tion  time  (/  =  0.5,  1,  2  and  3  h)  was  taken  on  the  base 
of  the  time  required  to  reach  the  the  refluxing  tem¬ 
perature.  After  each  deposition,  the  coated  sample  was 
immediately  removed  from  the  mixture,  rinsed  with 
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Fig.  1.  Dependence  of  Ni  film  thickness  on  deposition  time. 

distilled  water  and  acetone  or  2-propanol,  and  subse¬ 
quently  air-dried  for  characterization. 

The  structure  of  the  films  was  studied  using  x-ray 
diffraction  (XRD)  in  the  9-29  mode.  The  diffraction 
peaks  were  analyzed  using  a  profile  fit  routine  and  the 
FWHM  of  the  specular  peaks  were  used  for  estimation 
of  average  crystallite  size  and  strain  of  these  films.  It 
was  assumed  that  the  total  line  broadening  J  was 
the  linear  sum  of  that  due  to  size  and  strain.  A  plot  of 
(7ijM  i/cos0)  v.v  sin  0  yielded  the  strain  and  the  size  from 
the  slope  and  intercept,  respectively.  The  microstruc¬ 
ture  was  investigated  using  scanning  electron  micros¬ 
copy  (SEM)  and  transmission  electron  microscopy 
(TEM).  Cross-sectioned  samples  were  checked  with 
SEM  for  film  thickness  determination.  The  hysteresis 
loops  were  measured  using  vibrating  sample  magne- 
tometry  (VSM)  in  the  parallel  (//)  and  perpendicular 
(1)  directions  to  the  film  plane  using  a  maximum  mag¬ 
netic  field  of  ±  50kOe. 

3.  RESULTS  AND  DISCUSSIONS 

The  dependence  of  film  thickness  on  1  and  the  [Pt]/ 
[Ni]  is  shown  in  Fig.  1 .  For  /  =  3  h,  the  decrease  in  film 
thickness  may  be  related  to  the  corrosive  attack  of  eth¬ 
ylene  glycol  on  deposited  Ni  films.  It  has  been  reported 
for  other  systems  involving  chemical  oxidation  of  eth¬ 
ylene  glycol  at  elevated  temperatures,  oxalic  acid  was 
identified  as  a  product.  Both  it  and  various  organic 
acids  that  were  formed  as  intermediates  are  corrosive 
to  metals  [13],  Nickel  and  Pt  are  miscible  below  200 
°C.  When  Pt  was  used  as  nucleating  aid,  the  reaction 
completed  at  about  0.5  h  for  [Pt]/[Ni]  of  1.5- 10'2.  When 
[Pt]/[Ni]  <  2- 1  O'5,  the  reaction  took  about  1  h  to  finish. 
The  use  of  Pt  precursors  for  aiding  nucleation  of  Ni 
nanocrystals  has  been  reported  [12],  The  similar  results 
were  reported  for  Pd(Il)  and  Cu(ll)  acetates  in  2- 
ethoxyethanol  at  135  °C,  in  which  initial  Pd  reduction 
assisted  the  reduction  of  copper  [14],  In  this  work,  Pt 


Fig.  2.  Dependence  of  crystallite  size  and  strain  on  deposition 
time. 


precursor  only  promoted  rapid  precipitation  of  Ni  pow¬ 
der  but  it  did  not  enhance  film  deposition,  since  thin¬ 
ner  films  were  obtained  compared  to  those  deposited 
without  Pt  (Fig.  1).  The  undesirable  competition  be¬ 
tween  powder  precipitation  and  film  formation  can  be 
minimized  by  application  of  electric  field  during  depo¬ 
sition  and  lower  deposition  temperature  [15], 

The  XRD  results  showed  that  Ni  was  polycrystal¬ 
line  fee.  The  average  crystallite  size  and  strain  of  Ni 
films  are  shown  in  Fig.  2.  The  crystallite  size  reached 
a  maximum  of  about  138  nm  at  /  =  2  h,  but  it  signifi¬ 
cantly  decreased  to  40  nm  at  /  =  3h.  The  decrease  in 
crystallite  size  and  film  thickness  at  longer  t  was  the 
result  of  corrosive  attack  of  the  solvent  on  the  depos¬ 
ited  films,  leading  to  dissolution  and  re-deposition  of 
finer  particles.  The  micro-strain  of  these  films  changed 
from  tensile  to  compressive  with  increasing  /.  Fig.  3 
shows  the  SEM  micrographs  of  Ni  films  deposited  at  t 
=  2  and  3  h,  respectively.  It  can  be  seen  that  sintering 
of  films  occurred  (thus  largest  ciy  stall  ile  size  as  shown 
in  Fig.  2)  at  t  =  2  h,  whereas  finer  particles  were  ob¬ 
served  in  film  deposited  at  t  =  3  h.  Note  that  the  par¬ 
ticles  seen  in  SEM  consisted  of  nanostructured 
crystallites  (measured  by  x-ray  line  broadening).  Fig. 
4  shows  a  TEM  micrograph  of  Ni  crystallites  on  Cu 
substrate  deposited  at  /  =  2h. 

The  XRD  patterns  of  Ni  films  prepared  using  Pt 
catalysts  only  showed  only  Ni,  without  any  evidence  of 
Pt  as  a  separate  phase.  The  average  crystallite  size  (es¬ 
timated  from  (111)  line  broadening)  increased  with 
decreasing  [Pt]/[Ni]  precursor  molar  ratio  from  33  to 
about  70  nm.  The  similar  results  of  dependence  of  Ni 
particle  size  on  Pt  nucleating  aid  was  reported  in  the 
polyol  synthesis  of  Ni  powders  [12], 

The  Ni  films  were  magnetically  saturated.  Both  Ni 
and  Pt-catalyzed  Ni  films  exhibited  in-plane  magneti¬ 
zation  anisotropy.  Fig.  5  shows  the  A/  and  MiL  for  Ni 
films  deposited  at  0.5  and  3  h.  Fig.  6  shows  the  depen- 
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Fig.  3.  SEM  micrographs  of  Ni  films  deposited  at  :  a)  2  h. 
b)  3  h. 


100  nm 

Fig.  4.  TEM  bright  field  micrograph  of  Ni  film  on  Cu 
substrate  deposited  at  2  h.  Dark  particles  arc  Ni. 

dence  of  M  .  FI  and  H  .  of  Ni  films  on  t.  The  M  of 

s .  •  C  t’l  s 

bulk  fee  Ni  is  about  M  is  484  611111/0111’.  The  values  of 
M  increased  with  t  and  reached  405  emu/cm’  for  the 
Ni  film  deposited  3  h.  This  may  be  due  to  the  increased 
film  density  at  longer  t.  For  films  consisting  of 
randomly  oriented  crystallites,  the  squareness  ratio 
(ratio  of  remanent  magnetization.  A/,  to  saturation 
magnetization,  M)  is  typically  about  0.5.  This  ratio 
was  between  0.4  and  0.5  for  Ni  films,  except  it  was 
higher  for  that  deposited  at  1  h  (about  0.6).  Fig.  7  shows 
M  and  the  squareness  ratio  of  Ni  films  deposited 
with  Pt  catalysts.  M  .  increased  whereas  squareness  ratio 


Fig.  5.  Hysteresis  loops  of  Ni  films  deposited  at  0.5  h  and  3 
h,  showing  magnetic  saturation  and  in-plane  anisotropy. 


Fig.  6.  Dependence  of  saturation  magnetization  and  coercivity 
of  Ni  films  on  deposition  time. 


decreased  with  decreasing  [Pt]/[Ni].  The  deposition  time 
forNi  films  using  [Pt]/[Ni]  of2T0'5  and  2-10'6  was  about 
1  h,  and  was  0.5  h  for  [Pt]/[Ni]  of  1.5-1  O'2.  Pt  had  a 
predominant  effect  on  Ms//and  a  lower  [Pt]/[Ni]  ratio 
led  to  a  higher  Ms  as  expected.  Lower  Af  may  be  related 
to  alloying  of  Pt  with  Ni.  Examination  of  XRD  data  of 
these  films  showed  a  (200)  texture.  This  is  consistent 
with  the  higher  squareness  ratio  (between  0.56  to  0.65) 
of  these  films.  Fig.  8  shows  that  the  films  deposited 
with  Pt  catalysts  also  had  higher  perpendicular 
coercivity.  The  coercivity  anisotropy  may  be  related  to 
the  microstructures  of  deposited  film  [15],  The  work  is 
currently  underway  to  investigate  the  structure  and 
microstructure  of  these  films  by  transmission  electron 
microscopy. 

4.  SUMMARY 

Nanostructured  Ni  films  were  deposited  by  the  polyol 
method  on  Cu  substrates  with  or  without  Pt  as  nucleat¬ 
ing  aid.  The  solution  chemistry  varied  with  t.  At  longer 
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Fig.  7.  Dependence  of  saturation  magnetization  and  square 
ratio  of  Ni  films  on  [Pt]/[Ni]  precursor  molar  ratio. 


I,  corrosive  attack  of  the  solvent  on  deposited  films  led 
to  dissolution  and  re-deposition  of  Ni.  Saturation  mag¬ 
netization  increased  with  /  and  decreased  with  increas¬ 
ing  Pt  precursor  concentration.  The  magnetically  satu¬ 
rated  films  showed  in-plane  magnetization  anisotropy 
and  perpendicular  coercivity  anisotropy.  The  use  of  Pt 
favored  rapid  powder  precipitation  but  did  not  enhance 
film  deposition.  The  Pt-catalyzed  Ni  films  showed  a 
(200)  texture  and  a  higher  squareness  ratio.  Further 
work  is  needed  to  understand  the  solution  chemistry 
and  film  deposition  chemistry  in  the  polyol  process. 
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Abstract.  Under  certain  growth  conditions  InAs/Si  heteroepitaxial  growth  proceeds  via  Stranski-Krastanow 
or  Volmer- Weber  growth  modes  depending  on  the  growth  parameters.  The  critical  thickness  at  which  three 
dimensional  InAs  islands  start  to  appear  at  the  Si(100)  surface  is  within  the  range  of  0.7-4. 0  monolayers 
(substrate  temperature  range  is  350  °C  -  430  °C).  Their  size  depends  critically  on  the  growth  conditions 
and  is  between  5  nm  and  80  nm  (uncapped  islands).  Critical  lateral  size  of  the  coherent  (Si  capped)  dislocation- 
free  island  is  equal  to  2-5  nm  depending  on  the  island  height.  Islands  having  larger  size  are  dislocated. 

Optical  properties  of  InAs  nanoscale  islands  capped  with  Si  reveal  a  luminescence  band  in  the  1.3  pm 
region. 


1.  INTRODUCTION 

Self-organisation  effects  on  semiconductor  surfaces 
which  occure  in  the  course  of  molecular  beam  epitaxy 
(MBE)have  attracted  much  attention  over  the  last  ten 
years.  A  particular  interest  was  shown  in  the  quantum 
dots  (QDs)  in  III/V  -  HI/V(e.g.,  InAs  -  GaAs)  [1,2], 
II/VI  -  II/VI  (e.g,  CdSe  -  ZnSe)  [3]  or  IV  -  IV  (Si  - 
Ge)  [4]  materials  systems.  This  system  stimulates  much 
investigations  because  silicon  remains  a  key  material 
in  the  modem  microelectronics  industry.  For  optoelec¬ 
tronic  applications  silicon  is  not  well  suitable  because 
of  its  indirect  band  gap  nature.  Attempts  to  improve 
the  situation  using  SiGe  -  Si  QDs  did  not  lead  to  a 
significant  progress  as  these  nanostructures  provide 
indirect  band  alignment  both  in  k  -  and  in  real  space. 
In  order  to  increase  luminescence  efficiency  of  silicon- 
based  structures  we  propose  to  insert  direct  bandgap 
InAs  QDs  in  a  Si  matrix  using  MBE  growth  [5].  We 
emphasize  that  this  approach  is  different  from  the 
growth  utilizing  InAs  QDs  on  thick  GaAs  buffers  [6], 
Small  coherent  InAs  islands  additionally  enable  a  pos¬ 
sibility  to  overcome  the  problems  of  strong  lattice 
mismatch  and  formation  of  antiphase  domains  inherent 
for  the  growth  of  thick  III  -  V  epilayers  (e.g.  GaAs)  on  Si. 

In  this  paper  we  report  on  the  MBE  growth  and 
properties  of  the  InAs  nanoscale  islands  formed  on  sili¬ 


con.  We  have  found  that  under  certain  growth 
conditions  InAs/Si  heteroepitaxial  growth  proceeds  via 
Stranski-Krastanow  or  Volmer- Weber  growth  modes 
depending  on  the  growth  parameters  leading  to  the 
formation  of  nanoscale  islands.  The  lateral  size  depends 
critically  on  the  growth  conditions  and  is  in  the  range 
of  2-80  nm.  Optical  properties  of  InAs  QDs  capped 
with  Si  reveal  a  luminescence  band  in  the  1.3  pm 
region.  A  pronounced  excitation  density  dependence 
of  the  photoluminescence  (PL)  peak  position  and  a 
decay  time  of  approximitely  400  ns  were  observed. 

2.  EXPERIMENTAL  DETAILS 

The  growth  experiments  were  carried  out  using  EP1203 
(Russia)  or  Riber  Supra  (France)  MBE  machines  on 
exactly  oriented  Si(100)  substrates.  The  Si(100)  sur¬ 
face  was  prepared  in  a  way  similar  to  that  described  in 
[7].  Thermal  desorption  of  silicon  native  oxide  was 
performed  in  a  growth  chamber  at  the  substrate  tem¬ 
perature  of  830  °C  -  870  °C  during  1 5  min.  Afterwords, 
well-resolved  (2x1)  or  mixed  (1x2)  and  (2x1)  surface 
reconstructions  typical  for  cleaved  Si(l  00)  surface  were 
observed.  Then  the  substrate  temperature  was  smoothly 
reduced  to  the  desired  value  and  the  InAs  deposition 
was  initiated  at  a  conventional  MBE  mode.  The  InAs 
deposition  rate  was  typically  0.1  monolayer  (ML)  per 
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Fig.  1.  Critical  thickness  dependence  on  substrate  temperature 
and  Ja.4/Ji„  fluxes  ratio  during  InAs/Si  growth. 

second.  When  the  deposition  of  the  desired  average 
thickness  of  InAs  had  been  completed,  the  sample  was 
immediately  quenched  to  the  room  temperature  and 
removed  from  the  growth  chamber  in  order  to  perform 
scanning  tunnelling  microscopy  (STM)  and  scanning 
electron  microscopy  (SEM)  studies.  SEM 
measurements  were  performed  using  CamScan  setup. 
Pieces  for  STM  studies  were  covered  with  silicon 
vacuum  oil  immediately  after  exposure  to  the  atmo¬ 
sphere,  a  procedure  used  before  in  our  studies  of  InAs 
QDs  on  GaAs  surfaces  [8],  A  calibration  of  the  growth 
rate,  III  -  V  flux  ratio,  and  monitoring  of  the  surface 
morphology  during  growth  was  performed  using 
reflection  high  energy  electron  diffraction  (RHEED) 
system  composed  of  a  high  sensitivity  video  camera,  a 
video  tape  recorder  and  a  computer,  all  interconnected 
via  specially-designed  interface  [9].  Transmission  elec¬ 
tron  microscopy  (TEM)  images  were  obtained  on  the 
samples  when  InAs  islands  were  covered  with  30-50 
nm  of  silicon  using  JEOL  4000EX  microscope.  The 
same  structures  were  used  for  luminescence  experi¬ 
ments.  A  514.5  nm  line  of  the  Ar+  laser  was  used  as  an 
excitation  source  and  a  Ge  photodiode  as  a  detector. 
The  excitation  density  was  5-100  W/cm2.  The  struc¬ 
tures  for  TEM  and  PL  studies  were  grown  in  a  follow¬ 
ing  manner.  Just  after  InAs  deposition,  a  10-20  nm  Si 
cap  layer  was  grown  at  the  same  substrate  temperature 
T  as  for  InAs  deposition  followed  by  a  1 0  min  annealing 
procedure  at  650  °C  —  700  °C.  Then  20-40  nm  of  Si 
was  grown  at  the  same  temperature  with  optional 


consequent  10  min  annealing  at  700  °C  -  800  °C  in 
order  to  smooth  the  resulting  surface. 

3.  RESULTS  AND  DISCUSSION 

In  Fig.  1  we  summarise  RHEED  data  of  the  dependen¬ 
cies  of  the  substrate  temperatures  and  As4/In  fluxes  ratio 
in  the  ranges  350  °C  -450  °C  and  2-12,  respectively, 
on  the  critical  thickness  dcrji  (mean  thickness  at  which 
3D  nanoscale  islands  start  to  appear  at  the  surface). 
The  lowest  limit  of  Ts  was  chosen  for  two  reasons: 
first,  determination  of  7)  at  very  low  temperatures  is 
difficult  and  second,  because  of  poor  quality  of  Si  cap 
layer  grown  at  too  low  temperatures.  Poor  quality  of 
the  Si  cap  layer  may  hinder  the  device  applications. 
On  the  other  hand,  we  have  found  that  at  temperatures 
higher  then  450  °C  no  nanoislanding  is  observed  on 
the  surface  independently  on  the  V/III  flux  ratio,  total 
arsenic  overpressure,  etc.  limiting  the  range  for  the 
highest  T  as  indicated  in  Fig.l.  (At  higher  Ts  only 
mesoscopic  dislocated  clusters  with  lateral  size  of~400 
nm  appeared  on  the  Si(100)  surface  [10]).  Following 
the  data  presented  in  Fig.l  one  can  conclude  that  de¬ 
pending  on  the  growth  conditions,  2D-3D  transition 
can  be  tuned  in  the  thickness  range  of  0.7-4. 0  ML.  On 
the  other  hand,  we  have  observed  two  different  growth 
modes  for  the  InAs/Si  heteroepitaxial  growth.  First, 
Vo lmer- Weber  growth  mechanism  could  be  realised 
when  InAs  nanoscale  islands  are  formed  on  the  bare 
Si  substrate  (i.e.  for  dir<  1  ML).  Typical  RHEED  pat¬ 
tern  for  this  case  is  presented  in  Fig.  2, a.  At  higher  Ts 
(400  °C  -  450  °C)  we  observed  Stranski-Krastanow 
growth  mechanism  when  dmi  exceeded  1.0  ML  and 
the  formation  of  InAs  nanoscale  islands  occurs  on  the 
top  of  the  wetting  layer.  RHEED  pattern  for  this  par¬ 
ticular  case  is  presented  in  Fig.  2,b.  These  results  are 
in  agreement  with  the  2D  -  3D  transition  intervals  ob¬ 
served  by  other  groups  [11,  12]. 

Typical  surface  morphology  extracted  from  the  STM 
images  for  these  cases  revealed  high  density  arrays  of 
nanoscale  islands  (with  surface  density  of  the  order  of 
((l-r6)x  10"  cm'2)  and  the  lateral  sizes  of  10-20  nm 
[13,  14],  These  islands  are  basically  pyramidal  in  shape 
with  rectangular  or  triangular  base. 

We  have  found  that  the  growth  conditions  (the 
amount  of  InAs  deposited,  fluxes  ratio  influenced 
greately  the  arrangement,  sizes  and  surface  density  of 
InAs/Si  islands  [15].  In  order  to  clarify  how  growth 
parameters  influence  the  surface  morphology  we  grew 
three  samples,  labelled  1—3.  Their  MBE  growth 
conditions  and  geometrical  characteristics  measured 
with  SEM  are  presented  in  the  Table  1 . 

The  main  trends  for  the  uncapped  InAs/Si  island 
formation  are  the  following.  At  higher  As  pressure  and 
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Fig.  2.  Typical  RHEED  patterns  in  [Oil]  direction  taken 
after  the  formation  of  InAs  3D  islands  in  Si(100)  surface  via 
(a)  Volmcr-Weber  and  (b)  Stranski-Krastanow  growth 
modes. 

the  same  InAs  mean  thickness  the  island  size  drops  as 
compared  to  the  case  of  the  decreased  As/In  fluxes  ratio 
(samples  1  and  3).  We  believe  that  this  is  due  to  the 
suppressing  of  the  surface  migration  length.  The  den¬ 
sity  of  the  islands  for  the  sample  3  is  four  times  higher 


than  that  of  the  sample  1.  The  islands  size  distribution 
becames  narrower,  too.  With  increasing  the  average 
InAs  thickness  (sample  2)  the  islands  start  to  form  large 
conglomerates.  These  clusters  are  elongated  and  ex¬ 
hibit  crystallographic  shape  with  the  base  oriented 
along  [01 1]  and  [OlT]  directions.  The  similar  situation 
was  observed  in  [11]  where  such  clusters  were  formed 
during  the  InAs  growth  on  Si(IOO)  passivated  (hydro¬ 
gen-terminated)  surfaces  with  approximately  same  size 
and  surface  density.  We  note  that  for  the  sample  2  the 
critical  thickness  is  larger  by  the  value  of  1 .8  ML,  i.e. 
2.5  times  higher.  In  InAs/GaAs  heteroepitaxial  system 
such  large  sizes  should  lead  to  the  formation  of  the 
dislocations.  In  order  to  check  the  crystalline  quality 
we  studied  the  same  uncapped  samples  with  cross — 
sectional  TEM. 

Most  of  the  islands  exhibited  considerably  good 
crystalline  quality  although  misfit  dislocations  appeared 
at  the  InAs/Si  interface.  In  Fig.  3  typical  cross-section 
TEM  image  for  single  InAs  nanoisland  is  presented. 
Using  high  resolution  TEM  data  taken  from  the  Si 
capped  samples  we  have  found  that  the  critical  lateral 
size  of  the  coherent  dislocation — free  islands  is  equal 
to  2-5  nm  depending  on  the  island  height.  Islands  of 


Fig.  3.  Typical  high  resolution  cross-section  TEM  image  of 
the  single  InAs  nanoscalc  island  at  Si(100)  surface. 


Tabic.  MBE  growth  conditions  and  geometrical  characteristics  of  the  InAs /Si(100)  quantum  dot  samples 


Sample 

No. 

InAs  thickness. 

ML 

As/In  fluxes 
ratio 

Island  lateral 
size ,  nm 

Surface  density  of 
islands,  10'°  cm2 

I 

1.2 

3 

15-50 

4 

2 

2.5 

3 

20-80 

1.6 

3 

1.2 

10 

3-15 

16 

18 
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Fig.  4.  PL  data  for  the  InAs  nanostructures  embedded  in  Si 
via  Volmcr- Weber  growth  modes. 


larger  size  are  defected  due  to  misfit  dislocations  or 
twins.  We  also  suggested  that  most  of  the  InAs 
reevaporated  during  the  deposition  (or  annealing)  of 
the  cap  layer.  Fig.  3.  shows  typical  high  resolution 
cross-section  TEM  image  of  the  single  InAs  nanoscale 
island  at  Si(100)  surface. 

The  samples  with  such  InAs  quantum  dots  capped 
with  30-50  nm  Si  layer  show  a  luminescence  band  in 
the  1.3  pm  region  for  temperatures  up  to  -170K.  The 
PL  exhibited  a  pronounced  blue  shift  with  increasing 
the  excitation  density  and  decays  at  the  time  constant 
of -400  ns  [16].  Fig.4  compares  low  temperature  PL 
spectra  of  the  4.8  ML  with  the  1.6  ML  InAs  samples 
deposited  in  different  growth  modes  (Stranski- 
Krastanow  for  4.8  ML  and  Volmer- Weber  for  1 .6  ML'. 
These  spectra  reveal  TO-phonon  assisted  exciton 
emission  of  the  Si  substrate  and  barrier  at  ~1.11  eV. 
The  broad  emission  band  in  the  0.9  to  1 .05  eV  region 
was  observed  for  both  InAs  deposition  thickness  sug¬ 
gesting  its  origin  to  be  associated  with  the  InAs 
nanoobjects  embedded  in  Si  matrix.  This  broad  line 
was  not  observed  for  the  samples  with  InAs  insertions 
smaller  than  the  critical  thickness  or  on  a  Si  substrate. 
However,  the  size  of  these  insertions  evaluated  using 
TEM  measurements  does  not  exceed  5  nm  [16]  which 


*  InAs  critical  thickness  for  both  samples  was  exceeded 
by  the  value  of  0.8  ML. 
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rix:  (a)  InAs  grown  via  Stranski-Krastanow,  (b)  InAs  grown 

is  definitely  not  enough  for  charge  localisation  and 
quantum  confinement. 

Our  very  recent  detailed  TEM  investigation  of 
these  samples  shows  that  besides  a  random  distribu¬ 
tion  of  InAs  molecules,  two  kinds  of  ordering  were 
detected.  Such  regions  (<()  <  6  nm)  are  characterised  by 
larger  InAs  concentration.  The  first  kind  shows  an 
ordering  of  InAs  molecules  in  (101),  (101)  planes  in¬ 
clined  to  the  [001]  growth  direction,  while  in  the  sec¬ 
ond  one  the  ordered  (110),  (llO)  planes  are  parallel  to 
the  [001]  growth  direction.  The  formation  of  such  a 
structure  results  in  energy  lowering  due  to  the  decrease 
of  the  number  of  mixed  Si-As  and  Si-In  bonds.  We 
believe  that  the  appearance  of  the  broad  PL  peak  is  the 
result  of  solid  solution  formation  and/or  ordering,  be¬ 
cause  the  incorporation  of  InAs  molecules  into  Si  ma¬ 
trix  changes  the  band  gap  of  material  and  possibly  shifts 
the  relative  positions  of  conduction  and  valence  bands. 
Detailed  results  on  TEM  study  will  be  presented  in  a 
separate  paper  [17]. 

4.  CONCLUSIONS 

We  have  demonstrated  that  InAs/S  i(  100) 
heteroepitaxial  growth  might  proceed  via  Stranski  - 
Krastanow  or  Volmer-Weber  growth  modes  depend¬ 
ing  on  the  MBE  growth  parameters.  Critical  thickness 
at  which  three  dimensional  InAs  islands  start  to  ap- 
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pear  at  the  Si(lOO)  surface  varies  within  the  range  of 
0. 7-4.0  monolayers  (substrate  temperature  is  350  °C  - 
430  °C,  As/In  fluxes  ratio  is  2-12).  Their  size  depends 
critically  on  the  growth  conditions  and  is  between  5 
nm  and  80  nm  (uncapped  islands).  Critical  lateral  size 
for  Si  capped  InAs  coherent  dislocation-free  islands  is 
equal  to  2-5  nm  depending  on  the  island  height.  Is¬ 
lands  of  larger  size  are  dislocated.  Optical  properties 
of  the  InAs  nanoscale  islands  capped  with  Si  reveal  a 
luminescence  band  in  the  1.3  pm  region  independently 
on  the  growth  mode  used  for  InAs  deposition.  The  ap¬ 
pearance  of  this  peak  can  be  caused  by  the  solid  solu¬ 
tion  formation  and/or  ordering  due  to  the  incorpora¬ 
tion  of  InAs  molecules  into  Si  matrix  during  post 
growth  annealing  which  changes  the  band  gap  of  ma¬ 
terial  and  possibly  shifts  the  relative  positions  of 
conduction  and  valence  bands.  InAs/Si  heteroepitaxial 
system  seems  to  be  promising  for  Si-based  optoelec¬ 
tronic  applications,  in  particular  for  short  distance  fi¬ 
ber  optics  and  integration  of  microelectronics  and  op¬ 
tical  parts  at  the  same  silicon  wafer. 

ACKNOWLEDGEMENTS 

The  authors  thanks  to  N.P.  Korneeva,  V.N.  Demidov, 
S.A.  Masalov  and  V.M.  Busov  for  their  kind  assistance 
during  experiments.  This  work  was  supported  by 
INTAS  Grant  No.  96-0242,  Russian  Foundation  for 
Basic  Research  Grant  No.  98-02- 1 83 1 7  and  Grant  No. 
99-02-16799,  Russian  Scientific  Programmes  “Phys¬ 
ics  of  Solid  State  nanostructures”  Grant  No.  98-2029 
and  “Methods  and  Devices  of  Micro-  and 
Nanoelectronics”  Grant  No  02.04.5. 1 .40. E. 46. 

REFERENCES 

[1]  L.  Goldstein,  F.  Glas,  J.Y.  Marzin,  M.N.  Charasse, 
and  G.Le  Roux  //  App/.Phys. Lett.  47  (1985)  1099. 

[2]  D.  Leonard,  M.  Krishnamurthy,  C.M.  Reaves, 

S.P.  Denbaars  and  P.M.  Petroff  //  Appl.Phys.Lett. 
63  (1993)  3203. 

[3]  M.  Strassburg,  V.  Kutzer,  U.W.  Pohl,  A. 
Floffmann,  I.  Broser,  N.N.  Ledentsov,  D. 

Bimberg,  A.  Rosenauer,  U.  Fischer,  D.  Gerthsen, 
I.L.  Krestnikov,  M.V.  Maximov,  P.S.  Kop’ev  and 
Zh.I.  Alfero v  II  Appl.  Phys.  Lett.  72  (1998)  942. 

[4]  Feng  Liu  and  M.G.Lagally  //  Surf.  Sci.  386 
(1997)  169. 

[5]  N.N.  Ledentsov  in:  Proceedings  of  the  23rd 
International  Conference  on  the  Physics  of 
Semiconductors,  Berlin,  Germany,  July  21-26, 


1996,  edited  by  M.  Scheffler  and  R.  Zimmermann 
(World  Scientific,  Singapoure,  1996)  1  p.  19. 

[6]  J.M.  Gerard,  O.  Cabrol  and  B.  Sermage  // 
Appl.Phys.Lett.  68  (1996)  3123. 

[7]  A.  lshisaka  and  Y.  Shiraki  //  J.  Electrochem.  Soc. 
133  (1986)  666. 

[8]  G.E.  Cirlin,  G.M.  Guryanov,  A.O.  Golubok, 

S.Ya.  Tipissev,  N.N.  Ledentsov,  P.S.  Kop’ev,  M. 
Grundmann  and  D.  Bimberg  II  Appl.  Phys.  Lett. 
67  (1995)  97. 

[9]  G.M.  Gur’yanov,  V.N.  Demidov,  N.P.  Korneeva, 
V.N.  Petrov,  Yu.B.  Samsonenko  and  G.E.  Tsyrlin 
//  Tech.  Phys.  42  (1997)  956, 

[10]  G.E.  Cirlin,  V.N.  Petrov,  V.G.  Dubrovskii,  S.A. 
Masalov,  A.O.  Golubok,  N.I.  Komyak,  N.N. 
Ledentsov,  Zh.I.  Alferov  and  D.  Bimberg  // 

Tech.  Phys.  Lett.  24  (1998)  290. 

[1 1]  T.  Mano,  H.  Fujioka,  K.  Ono,  Y.  Watanabe  and 
M.  Oshima  It  Appl.  Surf.  Sci.  130 — 132  (1998) 
760. 

[12]  P.C.  Sharma,  K.W.  Alt,  D.Y.  Yeh  and  K.L. 

Wang  //  Appl.  Phys.  Lett.  75  (1999)  1273. 

[13]  G.E.  Cirlin,  V.G.  Dubrovskii,  V.N.  Petrov,  N.K. 
Polyakov,  N.P.  Korneeva,  V.N.  Demidov,  A.O. 
Golubok,  S.A.  Masalov,  D.V.  Kurochkin,  O.M. 
Gorbenko,  N.I.  Komyak,  V.M.  Ustinov,  A. Yu. 
Egorov,  A.R.  Kovsh,  M.V.  Maximov,  A.F. 
Tsatsul’nikov,  B.V.  Volovik,  A.E.  Zhukov,  P.S. 
Kop’ev,  Zh.I.  Alferov,  N.N.  Ledentsov,  M. 
Grundmann  and  D.  Bimberg  //  Semicond.  Sci. 
Technol.  13  (1998)  1262. 

[14]  G.E.  Cirlin,  V.N.  Petrov,  V.G.  Dubrovskii,  Yu.B. 
Samsonenko,  N.K.  Polyakov,  A.O.  Golubok, 

S.A.  Masalov,  N.I.  Komyak,  V.M.  Ustinov, 

A. Yu.  Egorov,  A.R.  Kovsh,  M.V.  Maximov,  A.F. 
Tsatsul’nikov,  B.V.  Volovik,  A.E.  Zhukov,  P.S. 
Kop’ev,  N.N.  Ledentsov,  Zh.I.  Alferov  and  D. 
Bimberg  //  Semiconductors  33  (1999)  972. 

[15]  G.E.  Cirlin,  N.K.  Polyakov,  V.N.  Petrov,  V.A. 
Egorov,  Yu.B.  Samsonenko,  B.V.  Denisov,  V.M. 
Busov,  B.V.  Volovik,  V.M.  Ustinov,  Zh.I. 
Alferov,  N.N.  Ledentsov,  D.  Bimberg,  N.D. 
Zakharov  and  P.  Werner  //  Czech.  J.  Phys.  49 
(1999)  1547. 

[16]  R.  Heitz,  N.N.  Ledensov,  D.  Bimberg,  A.Yu. 
Egorov,  M.V.  Maximov,  V.M.  Ustinov,  A.E. 
Zhukov,  Zh.I.  Alferov,  G.E.  Cirlin,  I.P. 
Soshnikov,  N.D.  Zakharov,  P.  Werner  and  U. 
Gosele  //  Appl.  Phys.  Lett.  74  (1999)  1701. 

[17]  N.D.  Zakharov  et.al.,  in:  MRS’99  Fall  meeting 
proceedings,  in  press. 


Mater.Phys.Mech.l  (2000)  20-26 


FABRICATION  OF  2-D  AND  3-D  PHOTONIC  BAND-GAP 
CRYSTALS  IN  THE  GHz  AND  THz  REGIONS 


G.  Kiriakidis  and  N.  Katsarakis 

Institute  of  Electronic  Structure  and  Laser  (IESL).  Foundation  for  Research  and 
Tcchnology-I  Iellas  (FO.R.T.H.).  P.O.  Box  1527.  Heraklion  711  10,  Greece 


Received:  October  1 999 

Abstract.  I\vo-  and  three-dimensional  dielectric  and  metallic  structures  exhibiting  photonic  band-gaps  in 
a  broad  frequency  range  were  fabricated  by  deep  X-ray  lithography  and  laser  micro-machining.  These 
techniques  seem  promising  for  mass  production  of  photonic  structures  with  full  band-gaps  in  a  spectrum 
ranging  from  the  millimeter  wave  up  to  the  far  infrared  regime.  Deep  X-ray  lithography  was  applied  to 
produce  periodic  3-D  photonic  band-gap  structures  based  on  the  “three-cylinder"  model  with  mid-gap 
frequencies  up  to  2.4  THz.  Metallic  3-D  structures  with  incorporated  point  and  linear  defects  arc  currently 
under  development.  Layered  metallic  and  mctallo-dielectric  structures  exhibiting  a  cutoff  frequency  in  the 
microwave  regime  were  fabricated  by  laser  precision  machining.  The  observed  cutoff  frequency  can  be 
easily  tuned  by  varying  the  interlayer  distance  or  the  filling  fraction  of  the  metal.  Combinations  of  layers 
with  different  metal  filling  fractions  create  defect  modes  with  relatively  sharp  peaks,  which  arc  also  tunable. 
The  mctallo-dielectric  structures  arc  significantly  smaller  than  the  simple  metallic  ones.  The  experimental 
measurements  seem  to  be  in  good  agreement  with  theoretical  calculations. 


1.  INTRODUCTION 

Over  the  last  decade,  a  new  brand  of  artificial  struc¬ 
tures  has  raised  the  excitement  of  scientists.  In  anal¬ 
ogy  to  electrons  in  a  crystal,  the  propagation  of  elec¬ 
tromagnetic  waves  travelling  in  periodic  dielectric 
structures  can  lead  to  the  formation  of  stop  bands  [1,2]. 
When  these  stop  bands  are  wide  enough  to  overlap  for 
both  polarization  states  along  all  crystal  directions,  the 
material  displays  a  complete  Photonic  Band-Gap 
(PBG),  that  is,  a  frequency  range  where  the  propaga¬ 
tion  of  electromagnetic  waves  is  forbidden,  irrespective 
of  the  propagation  direction.  These  effects  lead  to  some 
novel  applications  in  quantum  optics  in  the  frequency 
range  from  microwave  to  the  optical  regime  such  as 
low-cost  efficient  microwave  antennas,  low-loss  cavities 
and  resonators  and  zero-threshold  lasers  [3-7].  Photo¬ 
nic  crystals  represent  also  a  novel  and  promising 
improvement  on  the  operating  efficiency  of  either  Laser 
Diodes  or  Light  Emitting  Diodes  (LEDs).  They  have 
already  been  recognized  for  their  ability  to  control  the 
flow  of  light  and  their  capacity  to  modulate,  suppress, 
concentrate  and  direct  it,  in  an  enhanced  light-matter 
interaction  manner. 

The  potential  of  these  artificial  structures  has  been 
exploited  by  theoretical  and  experimental  studies 
mainly  on  2-D  and  3-D  photonic  crystals  operating  from 


the  microwave  to  the  optical  regime  utilizing  a  number 
of  diverse  fabrication  techniques  and  materials.  It  is 
the  purpose  of  this  work  to  present  an  inside  of  two 
selected  fabrication  techniques  for  2-D  and  3-D 
crystals,  namely  LIGA  [8]  and  laser  machining  [9], 
and  explore  their  potentials,  limitations  and  applica¬ 
tions  utilizing  both  dielectric  and  metallic  materials. 

2.  2-D  PHOTONIC  CRYSTALS 

The  critical  factors  to  control  the  properties  of  2-D 
photonic  crystals  are  the  refractive  index  contrast,  the 
arrangement  of  the  lattice  elements  and  the  fraction  of 
high  and  low  index  materials  in  the  lattice.  In  prin¬ 
ciple,  the  higher  the  refractive  index  contrast  the  stron¬ 
ger  the  expected  effects  and  the  larger  the  band-gap 
achievable.  Fabrication  of2-D  structures  has  flourished 
due  to  the  fact  that  they  are  relatively  easy  to  make 
utilizing  tools  and  techniques  borrowed  from  the  sili¬ 
con  microelectronics  industry.  The  most  significant 
progress  so  far  has  been  reported  with  planar  struc¬ 
tures  patterned  by  standard  photolithography  or 
electron-beam  lithography,  generating  patterns  in  the 
visible  and  near-IR  ranges  [10-12],  More  recently, 
multiple  exposure  holography,  a  technique  well  estab¬ 
lished  for  some  time,  has  been  applied  successfully  to 
fabricate  advanced  structures  [13],  The  importance  of 
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producing  low  dimensionality  structures  relates  to  the 
ease  of  manipulating  for  the  photonic  lattice  particu¬ 
larly  in  as  far  as  interfacing  with  optical  elements  such 
as  fibres,  light  sources,  detectors  and  waveguides  are 
concerned.  A  typical  2-D  photonic  lattice  consists  of 
holes  etched  into  a  semiconductor  substrate.  One  type 
of  lattice  is  triangular,  self-similar  upon  60°  rotations 
with  two  main  directions  of  symmetry  in  the  “honey¬ 
comb”  configuration  [14].  This,  along  with  the  “graph¬ 
ite”  type  have  demonstrated  the  widest  band-gaps  when 
regions  of  semiconductors  are  surrounded  by  large 
regions  of  air,  i.e.  when  the  area  fill-factor  is  low  [15]. 
Under  the  above  configurations  the  criterion  for  maxi¬ 
mum  effect  is  achieved  when  the  optical  path  length  in 
both  materials  is  roughly  equal.  However,  an  additional 
important  factor  that  determines  the  properties  of  2-D 
PBGs  is  that  of  light  confinement  in  the  third  dimen¬ 
sion.  Most  successful  structures  presented  to  date  uti¬ 
lize  some  forms  of  waveguide  confinement.  In  such 
configurations  there  are  a  number  of  obvious  advan¬ 
tages  related  with  the  application  of  epitaxial  layer 
growth  techniques  and  compatibility  with  established 
planar  optoelectronic  elements  [16].  Advanced 
configuration  designs  have  recently  demonstrated  PBG 
structures  with  high  transmission,  reflection  and  very 
sharp  filtering  characteristics,  promising  interesting 

future  optoelectronic  applications  [17,18]. 

3.  3-D  PHOTONIC  CRYSTALS 

The  first  experimental  evidence  of  a  3-D  photonic 
crystal  was  presented  by  E.  Yabionovitch  and 
co-workers  [4,5]  following  the  theoretical  work  of  K. 
M.  Ho  and  co-workers  [6]  in  1990.  This  was  an  array 
of  holes  drilled  into  high  refractive  index  material 
(since  called  the  “Yablonovite”)  showing  a  stop  band 
for  the  transmission  of  microwave  radiation  that  ex¬ 
tended  from  10  to  13  GHz  irrespective  of  the  direction 
of  propagation.  Furthermore,  it  was  demonstrated  that 
by  including  a  defect  either  by  adding  or  removing 
material,  localized  states  in  the  band-gap  were  created, 
in  direct  analogy  to  “doping”  effects  for  semiconduc¬ 
tors.  For  the  microwave  regime  the  immediate  appli¬ 
cation  for  such  an  effect  could  be  related  to  couple  a 
PBG  crystal  with  an  antenna  and  allowing  more  flex¬ 
ibility  and  enhancement  of  the  radiation  process  in 
point-to-point  communications.  The  above 
proof-of-concept  demonstration  by  E.  Yabionovitch  is 
known  today  as  the  “three-cylinder”  model  and  pre¬ 
sents  a  structure  with  the  symmetry  of  the  diamond.  A 
corresponding  diamond  structure  was  also  proposed 
by  C.  T.  Chan  and  co-workers  [19].  Later,  the  same 
group  from  Iowa  State  University  designed  a  novel 
three-dimensional  “layer-by-layer”  structure  demon¬ 


strating  again  a  full  3-D  PBG  over  a  wide  range  of 
structural  parameters.  This  structure  consisted  of  layers 
of  one-dimensional  dielectric  rods  of  rectangular,  cir¬ 
cular  or  elliptical  cross  section.  The  band-gap  is  opti¬ 
mized  by  controlling  the  dielectric  contrast  of  the  ma¬ 
terials  been  used,  the  filling  factor  of  the  dielectric  and, 
of  course,  the  lattice  constant  of  the  PBG  structure. 
Although  it  was  fundamentally  important  to  demon¬ 
strate  the  existence  of  a  full  3-D  PBG  in  the  micro- 
wave  regime,  the  drive  was  soon  focused  towards  the 
fabrication  of  truly  “photonic”  3-D  structures  that  op¬ 
erate  in  the  near  IR  and  eventually  the  visible  regime 
of  the  electromagnetic  spectrum.  The  obvious  approach 
was  to  scale  down  fabrication  techniques  successfully 
employed  in  the  microwave  regime.  Indeed,  a  signifi¬ 
cant  improvement  was  first  demonstrated  by  G.  Feiertag 
and  co-workers  [8],  showing  the  potential  of  combining 
deep  X-ray  lithography  (LIGA  technique)  with  the 
“three-cylinder”  model  to  produce  3-D  structures  with 
lattice  constants  of  227  and  114  pm,  corresponding  to 
mid-gap  frequencies  of  0.75  and  1 .5  THz  respectively. 
Detailed  account  of  the  LIGA  processing  results  will 
be  presented  in  the  following  section.  Since  then,  the 
“layer-by-layer”  model  was  matured  further  to  a  point 
of  reaching  into  the  far-IR  producing  structures  dis¬ 
playing  a  stop-band  between  1.35  and  1.95  pm  [20]. 
The  use  of  microfabrication  techniques  to  produce  struc¬ 
tures  in  the  1.5  pm  window  has  failed  to  show  yet  the 
expected  sharp  spectral  response.  However,  the  tech¬ 
nique  has  been  very  useful  in  the  fabrication  of  Metal¬ 
lic  Photonic  Band-Gap  (MPBG)  structures  and  is  pre¬ 
sented  in  details  in  the  subsequent  chapters.  Finally, 
novel  concepts  that  are  based  upon  the  self-organization 
of  structures  show  the  promise  of  the  approach  rather 
than  providing  a  fully  convincing  demonstration. 

4.  The  LIGA  technique  for  PBG  structures 
in  the  far  infrared  frequency  range 

LIGA  is  a  well  established  technique  for  the  produc¬ 
tion  of  high  aspect  ratio  microstructures  [21,22],  It 
combines  deep  X-ray  lithography  (DXRL)  with  syn¬ 
chrotron  radiation,  electroforming  and  ceramic  moul¬ 
ding.  Deep  X-ray  lithography  with  synchrotron 
radiation  is  used  to  transfer  an  absorber  pattern  from 
an  X-ray  mask  into  a  thick  resist  material  deposited  on 
a  substrate.  The  mask  consists  of  an  X-ray  transparent 
membrane  (e.g.  Be,  Si,  Ti)  with  absorber  structures 
made  of  gold.  By  shadow  printing  the  mask  pattern  is 
transferred  into  the  resist  using  X-rays.  Synchrotron 
radiation  is  applied  because  of  its  high  intensity  and 
small  divergence.  In  this  way,  extremely  precise  mi¬ 
crostructures  with  high  aspect  ratios  can  be  realized. 
Due  to  the  extremely  high  depth  of  focus  the  lateral 
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Fig.  1.  3-D  photonic  band-gap  structure  fabricated  by 


LIGA  according  to  the  “three-cylinder”  model. 


deviations  of  the  resist  walls  are  as  small  as  0.05  pm 
per  100  pm  thickness.  PMMA  is  usually  applied  as 
resist  material.  It  has  very  good  properties  to  achieve 
vertical  walls  with  a  roughness  of  less  than  50  nm. 
Through  exposure,  the  resist  undergoes  a  decrease  in 
molecular  weight  and  becomes  soluble  in  an  organic 
developer,  if  the  deposited  dose  is  above  2  kJ/cm3.  Af¬ 
ter  development  of  the  exposed  resist,  the  gaps  of  the 
resist  relief  are  filled  with  metal  via  electroforming. 
The  extremely  accurate  metal  form  can  be  used,  e.g., 
for  plastic  moulding,  offering  inexpensive  microstruc¬ 
ture  replication.  A  melted  thermoplastic  or  a  casting 
resin  is  injected  through  holes  in  a  conducting  gate 
plate,  yielding  a  corresponding  plastic  structure  which 
is  fixed  on  by  the  remaining  resin  in  the  holes  of  the 
gate  plate  after  cooling  and  separating  the  mold  from 
the  mold  insert.  The  plastic  structure  obtained  in  the 
moulding  process  can  be  used  as  well  as  a  “lost  form” 
for  the  fabrication  of  ceramic  microstructures.  After 
filling  the  plastic  form  with  a  ceramic  paste,  a  dense 
ceramic  structure  is  formed  by  foliowing  conventional 
drying  and  firing  processing.  LIGA  is  a  generic  tech¬ 
nology  and  therefore  has  a  wide  range  of  applications 
in  fields  such  as  sensors  and  actuators,  optical 
engineering  (integrated  optics,  nonlinear  optics,  fibre 
optics),  electrical  connectors,  instrumentation  for  mini¬ 
mal-invasive  surgery,  biomedical  and  fluid-dynamical 
engineering  and  micro-filtration  systems. 

The  “three-cylinder”  mode!  was  selected  for  the  fab¬ 
rication  of  3-D  PBGs  using  LIGA.  It  is  easily  realized 
by  three  tilted  irradiations  of  a  PMMA  slab  using  an 
X-ray  mask  with  a  triangular  array  of  holes.  Mask  and 
resist  arc  tilted  35°  with  respect  to  the  synchrotron 
radiation  beam.  Between  the  irradiations,  the  tilted 


arrangement  of  mask  and  resist  must  be  rotated  each 
time  by  120°. 

A  process  that  comprises  two  main  steps  has  been 
chosen  for  the  fabrication  of  three-dimensional  PBG 
structures  in  the  far-IR  [8].  Deep  X-ray  lithography  is 
used  to  produce  a  mold  which  is  then  filled  with  a  ce¬ 
ramic  or  pre-ceramic  material  (Fig.  1).  It  was  first  tried 
to  fill  the  holes  with  a  ceramic  paste  but  the  resist  mold 
could  not  be  filled  completely.  Therefore,  it  was  de¬ 
cided  to  use  a  pre-ceramic  polymer  (polyvinylsilazane) 
that  can  be  transformed  into  a  silicon-carbonitride  ce¬ 
ramic  by  a  subsequent  pyrolysis  process. 

The  transmission  characteristics  of  a  structure  with 
a  lattice  constant  of  85  pm  and  rod  diameter  of  22  pm 
have  been  measured.  The  filling  ratio  of  the  ceramic 
rods  in  this  structure  was  35%  and  the  total  thickness 
of  the  structure  450  pm.  There  is  a  constant  drop  of 
the  transmission  which  is  attributed  to  the  absorption 
of  the  ceramic.  The  experimental  results  are  in  agree¬ 
ment  with  theoretical  calculations  (Fig.  2).  There  is 
also  a  well-defined  drop  of  the  transmission  at  around 
80  cm'1  (2.4  THz)  which  is  related  with  the  first 
band-gap  created  due  to  the  periodicity  of  the  struc¬ 
ture.  The  gap  over  the  mid-gap  ratio  for  this  stop  band 
is  Ad  o)  =0.35.  Transmission  measurements  for  a  simi- 
lar  photonic  structure  with  rods  diameter  of  31  pm, 
which  corresponds  to  a  filling  ratio  of  57%,  show  a 
gap  at  around  70  cm'1  (2.1  THz)  with  Ad w=0. 32. 

Recent  efforts  are  targeting  the  fabrication  of  me¬ 
tallic  2-D  and  3-D  structures  with  different  kind  of 
defects  in  the  far-IR  range  of  the  electromagnetic  spec¬ 
trum.  New  photosensitive  materials,  like  SU-8,  are 
applied  in  order  to  achieve  even  higher  aspect  ratios 
for  both  2-D  and  3-D  structures. 
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Fig.  2.  Measured  (solid  line)  and  calculated  results  for  a  photonic 
structure  with  85  }.lm  lattice  constant  and  22  |im  rods  diameter. 
In  the  calculations,  a  dielectric  constant  of  the  rods  with  real 
part  equal  to  3  and  imaginary  part  0.1  and  0.2  was  used  (doted 
and  dashed  lines  respectively). 


Fig.  3.  A  three-dimensional  image  of  a  layer-by-layer  metallic 
photonic  band-gap  structure. 


5.  LASER-MACHINED  METALLIC 
AND  METALLO-DIELECTRIC  PBG 
STRUCTURES 

Most  of  the  earlier  research  work  was  concentrated  on 
the  development  of  PBG  crystals  built  from 
frequency-independent  dielectrics.  At  lower  microwave 
and  millimeter  wave  frequencies,  however,  metals  act 
like  nearly  perfect  reflectors,  no  absorption  problems 
occur,  and  there  are  certain  advantages  of  introducing 
metals  to  photonic  crystals.  These  include  reduced  size 
and  weight,  easier  fabrication  and  lower  costs. 

Metallic  photonic  band-gap  crystals  (MPBG)  op¬ 
erating  in  the  microwave  frequency  range  were  fabri¬ 
cated  by  laser  precision  machining  [9],  They  consist  of 
stainless  steel  plates  with  a  tetragonal  lattice  of  holes 
and  a  lattice  constant  of  15  mm.  Transmission 
measurements  show  that  periodic  crystals  exhibit  a 
cutoff  frequency  in  the  8-18  GHz  range,  below  which 
no  propagation  is  allowed.  The  cutoff  frequency  can 
be  easily  tuned  by  varying  the  interlayer  distance  or 
the  filling  fraction  of  the  metal.  Combinations  of  plates 
with  different  hole-diameters  create  defect  modes  with 
relatively  sharp  peaks,  which  are  tunable. 

In  particular,  stainless  steel  plates  of  1  mm  thick¬ 
ness  were  drilled  using  a  Nd:YAG  laser  with  a 
wavelength  of  1064  nm  and  an  output  power  of  100  W 
at  a  CNC  working  station  with  four  degrees  of  free¬ 
dom.  The  laser  was  operated  in  pulsed  mode  with  a 
repetition  rate  of  120  Hz.  Oxygen  flow  was  necessary 
to  prevent  deposition  of  debris  on  the  metal  surface. 
The  distance  of  the  laser  focusing  system  from  the 
sample  surface  was  0.5  mm  and  the  cutting  velocity 


600  mm/min.  Holes  with  diameters  of  8,  10,  12  and  14 
mm  were  drilled  on  each  plate  in  a  square  arrangement 
with  a  centre-to-centre  distance  of  1 5  mm  (10-10  array). 
Four  layers  of  each  specific  structure  were  fabricated. 
Consequently,  periodic  structures  of  up  to  4  layers,  as 
shown  in  Fig.  3,  could  be  built  and  measured  while 
combinations  of  layers  with  different  hole-diameters 
could  also  be  used  to  study  defects. 

All  periodic  structures  behave  like  high-pass  fil¬ 
ters.  They  exhibit  a  cutoff  frequency  in  the  8-18  GHz 
frequency  range,  below  which  there  are  no  propagat¬ 
ing  modes.  Within  the  band-gap,  the  rejection  rate  of 
the  transmitted  EM  signal  ranges  between  5  and  10 
dB  per  layer.  The  transmission  characteristics  of  peri¬ 
odic  structures  consisting  of  three  or  four  metallic  layers 
with  a  tetragonal  lattice  of  12  mm-diameter  holes  at  a 
distance  of  7  mm  are  presented  in  Fig.  4.  It  can  be  seen 
that  a  structure  with  three  layers  exhibits  a  cutoff  fre¬ 
quency  at  12  GHz. 

The  cutoff  frequency  of  the  MPBG  structures  can 
be  easily  tuned  by  varying  the  distance  between  the 
metallic  layers.  In  particular,  by  increasing  the 
interlayer  distance  the  gap-edge  shifts  to  lower  frequen¬ 
cies.  In  Fig.  5,  it  is  demonstrated  that  the  cutoff  fre¬ 
quency  is  linearly  related  to  the  interlayer  distance. 
These  results  are  in  good  agreement  with  theoretical 
calculations  that  are  shown  in  Fig.  6.  The  filling  frac¬ 
tion  of  the  metal  also  affects  the  observed  cutoff  fre¬ 
quency.  A  decrease  in  the  hole-diameter,  which 
corresponds  to  an  increase  in  the  filling  ratio  of  the 
metal,  shifts  the  gap-edge  to  higher  frequencies.  This 
behaviour  is  shown  in  Fig.  7,  where  it  can  be  seen  that 
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Fig.  4.  Transmission  profile  for  a  periodic  structure  with  three 
or  four  layers.  The  lattice  constant  (within  one  layer)  is  15  mm, 
the  diameter  of  the  holes  is  12  mm,  and  the  distance  between 
the  layers  is  7  mm. 


Fig.  5.  Cutoff  frequency  as  a  function  of  the  interlayer  distance 
for  a  periodic  structure  similar  with  the  one  described  in  Fig.  2. 


indeed  the  cutoff  frequency  increases  as  the 
hole-diameter  drops. 

Defects  are  created  easily  by  replacing  the  middle 
layer  of  a  periodic  structure  with  a  layer  having  larger 
holes.  Thus,  in  this  way,  metal  deficiency  is  created  in 
the  middle  layer.  All  defect  structures  show  sharp  trans¬ 
mitted  modes  in  the  prohibited  frequencies  of  the  peri¬ 
odic  counterparts.  The  quality  factor  Q  of  the  defect 
modes  is  about  70.  Additionally,  the  introduction  of  a 
defect  in  the  MPBG  structure  shifts  the  gap-edge  to 
higher  frequencies.  In  Fig.  8,  defect  structures 
consisting  of  three  layers  with  hole-diameters  of  8-12-8 
mm  and  8-14-8  mm  are  compared  to  the  fully  periodic 
one  (8-8-8  mm).  The  removal  of  more  metal  from  the 
middle  layer  (hole-diameters  of  8-14-8  mm  instead  of 
8-12-8  mm)  moves  the  defect  mode  to  lower  frequen¬ 
cies  while  the  position  of  the  gap-edge  is  not  affected. 


This  demonstrates  that  not  only  the  cutoff  frequency 
but  also  the  frequency  of  the  defect  modes  can  be  tuned. 

By  changing  the  interlayer  distance  of  the  defect 
structure  with  hole-diameters  of  8-12-8  mm  from  7  mm 
to  9  mm,  a  shift  of  the  total  frequency  spectrum  to  lower 
frequencies  occurs,  as  has  been  also  observed  for  the 
periodic  counterparts.  It  is  important  to  note  that  this 
effect  does  not  change  the  distance  between  the  defect 
peak  and  the  gap-edge.  The  above  observations  have 
been  also  demonstrated  on  different  arrangements  such 
as  the  12-12-14-12-12  mm  structure. 

However,  the  size  of  the  simple  metallic  PBGs  is 
restrictive  for  their  application  in  microwave  technol¬ 
ogy.  Therefore,  it  is  essential  to  reduce  their  size  by 
introducing  other  dielectrics  rather  than  air. 
Metallo-dielectric  PBG  structures  operating  in  the  mi¬ 
crowave  frequency  range  and  especially  between  10 


Fig.  6  .  Calculated  transmission  profiles  for  periodic  structures 
consisting  of  three  layers  with  hole-diameters  of  12-12-12  mm 
at  a  distance  of  11,  9,  and  7  mm  (solid,  dashed  and  dotted 
line  respectively).  The  lattice  constant  within  one  layer  is  15 
mm  for  all  structures. 


Fabrication  of  2-D  and  3-D  Photonic  Band-Cap  Crystals  in  the  GHz  and  THz  Regions 
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Fig.  8.  Transmission  profiles  for  defect  structures  consisting  of 
three  layers  with  hole-diameters  of  8-12-8  mm  and  8-14-8 
mm  (dashed  and  dotted  line  respectively)  as  well  as  for  the 
periodic  structure  with  hole-diameters  8-8-8  mm  (solid  line). 
In  all  cases,  the  lattice  constant  within  one  layer  is  15  mm  and 
the  interlayer  distance  is  7  mm. 

and  18  GHz  were  fabricated  according  to  the 
layer-by-layer  model.  In  particular,  a  metal  (aluminum) 
pattern  was  applied  onto  Al,03  wafers.  The  pattern 
consists  of  a  square  lattice  of  metal  grids  with  a  lattice 
constant  of  5  mm.  Transmission  measurements  show 
that  periodic  structures,  consisting  of  up  to  four  wafers 
laterally  aligned  to  each  other,  exhibit  a  cutoff  frequency 
in  the  10-18  GHz  frequency  range,  below  which  no 
propagation  is  allowed  (Fig.  9).  The  cutoff  frequency 
can  be  easily  tuned  by  varying  the  interlayer  distance 
or  the  width  of  the  metallic  grids.  Combinations  of 
wafers  with  grids  of  different  width  create  defect  modes 
with  relatively  sharp  peaks,  which  are  tunable.  The  size 
of  the  PBG  structures  was  strongly  reduced  compared 
to  the  simple  metallic  meshes,  for  the  same  frequency 
characteristics,  by  using  a  supporting  dielectric  sub¬ 
strate  for  the  metallic  grids.  The  application  of  AfO, 
wafers  as  a  supporting  material  for  a  metallic  grid  pat¬ 
tern  leads  to  a  reduction  of  the  size  of  the  simple  MPBG 
structure  by  a  factor  of  3.  This  can  be  further  improved 
by  using  high  dielectric  microwave  materials.  Dielec¬ 
tric  materials  with  tunable  properties  may  open  new 
perspectives  for  the  application  of  “clever”  PBG  struc¬ 
tures  in  microwave  technology. 

Recent  activities  are  targeting  the  fabrication  of 
monolithic  MPBG  structures  consisting  of  microwave 
dielectrics  patterned  with  metallic  grids  for  applica¬ 
tions  in  microwave  technology. 

6.  CONCLUSIONS 

Laser  micro-machining  and  LIGA  are  two  very  useful 
techniques  for  the  fabrication  of  PBG  structures.  LIGA 
has  been  selected  for  its  ability  to  produce  accurate  and 
reproducible  3-D  PBGs  up  to  the  THz  spectral  range, 
while  laser  machining  has  been  selected  for  its  versa¬ 
tility  to  producing  effective  tuneable  2-D  metallic  struc¬ 
tures  in  the  millimeter  wave  regime.  Results  presented 


Frequency  (GHz) 

Fig.  9. Transmission  profile  for  a  periodic  structure  with  four 
A1  / A1203  layers.  The  lattice  constant  of  the  structure  is  5  mm, 
and  the  width  of  the  A1  wires  is  2  mm. 


have  demonstrated  the  advantages  of  these  two  tech¬ 
niques  for  high  quality,  reproducible  and  low-cost  PBGs 
and  their  potential  for  immediate  application  in  the 
production  of  novel  components  in  quantum  optics. 
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Abstract.  The  kinetics  of  the  growth  of  thin  metal  films  from  solutions  of  photoactive  chemical  coordination 
compounds  irradiated  by  light  at  a  certain  wavelength  has  been  theoretically  investigated.  In  the  framework 
of  our  mathematical  model  the  time  dependencies  of  the  basic  characteristics  of  the  film  growth  process 
(size  distribution  function,  density  of  nuclcation,  concentration  of  metal  on  the  surface)  have  been  computed. 
They  proved  to  be  in  accordance  with  the  experimental  data.  It  is  shown  that  the  flux  of  metal  (substance) 
from  the  solution  onto  the  surface  of  the  substrate  can  be  derived  from  the  stationary  diffusion  equation 
with  the  added  term  taking  into  account  the  influence  of  the  photochemical  substance  source.  It  may  be 
inferred  that  the  variation  of  light  intensity,  in  great  range,  keeps  the  kinetics  of  the  process  of  flint  deposition 
out  of  the  drastic  changes. 

1.  INTRODUCTION 

The  processes  of  metal  films  formation  and  growth  on 
solid  substrates  are  of  great  importance  from  both  theo¬ 
retical  and  practical  viewpoints  [1].  In  the  papers  [2- 
5],  the  process  of  obtaining  Ag,  Au,  Cu  films  on  flat 
quartz  substrates  in  the  solutions  of  photoactive 
chemical  coordination  compounds  was  investigated 
experimentally.  The  essence  of  this  method  is  the  fol¬ 
lowing.  Let  us  consider  a  solution  of  the  coordination 
compound  \M:L^Xm  ,  where  M:'  is  the  central  metal 
atom  (s+  is  the  oxidation  level);  L  denotes  innersphere 
ligands;  X  denotes  outersphere  counterions.  If  ligand 
to  metal  charge  transfer  (  )  occurs  in  this 

compound,  the  irradiation  of  solution  by  light  at  the 
wavelength  X  =  ^  results  in  the  reduction  of  the 
metal  atom.  By  carrying  the  irradiation  on,  the  central 
metal  atom  can  be  reducted  to  the  zeroth  valence  state. 

Thus,  we  can  write  the  scheme  of  the  entire  process  as 

[M"LJXm  — M"  +  L+  +  X~ , 

One  should  put  the  substrate  to  be  deposited  on  the 
bottom  of  the  vessel  filled  with  the  solution.  If  we  di¬ 
rect  a  light  beam  from  above,  its  intensity  decreases 
according  to  the  Lambert-Buger-Behr  law  (Fig.  1). 

Zeroth  valence  metal  atoms  generated  by  light  will  dif¬ 
fuse  towards  the  substrate  and  be  adsorbed  on  it.  When  the 
critical  concenteration  of  metal  on  the  substrate  surface  is 
reached,  the  growth  of  new  crystal  metal  phase  starts. 
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It  was  noted  in  [2-5],  that  the  essential  feature  of 
films  obtained  in  this  way  is  the  narrow  size  distribu¬ 
tion  of  grains  in  the  film.  In  those  papers,  on  the  basis 
of  the  Volmer-  Weber  theory  [6]  the  qualitative  expla¬ 
nation  of  results  was  proposed.  The  process  of  forma¬ 
tion  and  growth  of  a  new  phase  on  the  surface  is  the 
complicated  multistage  one  [7],  and  the  Volmer- Weber 
theory  does  not  provide  any  exact  quantitative  descrip¬ 
tion. 

In  our  opinion,  the  method  described  above  is  very 
constructive  and  interesting  for  problems  of  films  grow¬ 
ing.  In  this  work,  we  investigated  theoretically  the  time 
dependencies  of  the  basic  characteristics  of  the  grow¬ 
ing  film,  on  the  basis  of  modern  film  growth  ideas, 
presented  in  the  publications  [1,  7], 

2.  FLUX  OF  METAL  ATOMS  ONTO 
THE  SUBSTRATE 

In  this  work,  the  initial  stage  of  the  film  condensation 
on  the  solid  substrate  was  investigated.  The  substrate 
is  considered  to  be  a  homogeneous  continuous  surface. 
Thermal  effects  of  crystallization  and  adsorption  are 
neglected.  Passing  through  the  solution,  the  light  in¬ 
teracts  with  molecules  of  coordination  compounds.  That 
results  in  the  formation  of  zeroth  valence  metal  atoms. 
The  power  of  metal  atoms  source  is  given  by 

g(x)  =  IB<paexp[-ax],  (1) 
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Fig.  1.  The  scheme  of  the  process  of  photochemical  film 
deposition;  /  —  the  thickness  of  the  solution  layer,  I()  —  the 
incident  light  intensity. 


In  fact,  some  average  time-dependent  concentration 
of  adsorbed  atoms  on  the  substrate  exceeding  coexists. 
It  is  the  reduced  difference  of  these  concentrations  be¬ 
ing  referred  to  as  the  supersaturation 


£(0  = 


c  ( t)-i. 


(5) 


that  is  the  driving  force  for  the  phase  transformation 
on  surface.  However,  due  to  the  growth  is  carried  out 
at  room  temperature  and  the  solubility  of  metal  is  very 
low,  the  equilibrium  concentration  of  metal  c'  is  also 
low.  This  results  in  formation  of  surface  metal  islands 
at  small  supersaturations.  Therefore,  following  [1,  7, 
8],  we  assume  that  the  application  of  the  boundary 
condition  (4)  is  enough  satisfactory.  By  differentiation 
of  the  solution  of  the  Eq.  2  with  the  boundary  conditions 
(3)  and  (4),  we  can  derive  the  form  of  the  flux  of  metal 
to  the  surface 


J,  =  IM  l-exp(-a/)). 


(6) 


where  /  is  the  incident  light  intensity;  a  is  the  absorb- 
tion  coefficient;  y  is  the  quantum  yield.  The  absorption 
coefficient  a  can  be  represented  by  a  =  2.3s  c  .om/),cv,  where 
t  is  the  extinction;  c  .  is  the  concentration  of  ab- 
sorbing  centres  in  the  solution.  Here,  we  suppose  that 
c  ,  is  constant.  If  the  height  of  solution  /  (Fig.  1)  is 
small,  the  probability  of  the  formation  of  metal  clusters 
in  the  bulk  of  solution  (homogeneous  nucleation)  van¬ 
ishes.  In  this  case,  the  flux  of  metal  onto  the  substrate 
can  be  derived  from  the  ID  stationary  diffusion  equa¬ 
tion  written  by 

Z)Ac(x)  +  g(x)  =  0,  (2) 

where  A  denotes  the  Laplacian;  D  is  the  diffusion 
coefficient;  the  power  of  source  is  given  by  the  Eq.  1. 
We  make  use  of  the  stationary  diffusion  equation  for 
the  following  reasons: 

1 .  the  source  of  metal  atoms  is  stationary 

2.  we  suppose  that  the  characteristic  time  of  photo¬ 
chemical  reaction  is  much  less  than  that  of  diffu¬ 
sion  and,  i.e.  the  photochemical  reaction  is  not  the 
ratelimiting  stage. 

The  boundary  conditions  to  the  Eq.  2  are 


3.  METAL  ISLANDS  FORMATION 
PROCESSES 

According  to  the  paper  [7],  a  process  of  new  phase 
formation  goes  through  three  basic  stages.  The  first 
stage  is  the  fluctuation  formation  of  new  phase  nuclei. 
It  is  described  by  the  Zeldovich  theory,  modified  for 
the  formation  on  the  surface  [7].  During  this  stage,  the 
stationary  flux  of  nuclei,  which  are  spherical  segments, 
has  the  form  [7]: 


/(£)  =  +  1)  In(£+  1)  expf 


V 


ln2(£+l)J 


(7) 


where  %  is  the  supersaturation;  h=(a  /(A:7))34tiv,  a  is 
the  surface  tension;  v  is  the  volume  per  atom  in  nucleus; 
K  is  the  proportionality  coefficient  depending  on  the 
duffusion  coefficient  of  adatoms.  The  time  lag  is  esti¬ 
mated  by: 


_ 1 

l'  ~  -F"(nJW(nS 


(8) 


=  0,  (3) 

,v=0 

c(/)  =  c;.  (4) 

It  means  that  the  flux  of  metal  into  the  system  at  the 
upper  boundary  of  solution  (r=0)  misses  and  the 
concentration  of  metal  on  the  surface  (.v=l)  is  constant 
and  equal  to  c*. 


dc 

dx 


where  n  is  the  number  of  atoms  in  critical  nucleus; 

cr 

F(n)  is  the  Helmholtz  energy  of  nucleus  taking  form  of 
a  spherical  segment;  W(n)  is  the  diffusion  coefficient 
in  the  size  space  taking  into  account  the  probability  of 
random  eleminating  or  joining  atoms  to  nucleus  (the 
fluctuations  of  nucleus  in  the  size  space)  [9], 


F(nr)  = 


b 

ln(£+l) 


(9) 
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At  t  >  ts  the  second  stage  starts.  This  is  the  process 
of  formation  of  new  phase  islands  in  the  system.  This 
stage  is  described  by  the  system  of  equations: 

df{n\t)  ^  d(v(tr,t) /(/?;/)) 
dl  dn 

where  v(irj  )  =  is  the  linear  rate  of  nucleus  growth, 
written  as 

v(n\t)  =  a^{t)n"' ,  (11) 

where  a  and  m  are  coefficients  depending  on  the  mecha¬ 
nism  of  nuclei  growth, 

K  '  <  * 

- — - — 1  =  f  J  dl  -  k  f  \  f  {n\t)ndndt ,  (12) 

m  t  '  it 

where  k  is  the  coefficient  depending  on  the  metal  un¬ 
der  investigation.  The  equation  above  represents  the 
matter  conservation.  The  equation  (10)  should  be 
complemented  by  the  initial  and  boundary  conditions 
[7,  10] 


/(£(/)) 

/(«;  0)  =  0,  /(«„;/)=  (13) 

K«;0  K 

The  method  of  solution  of  this  equations  system 

was  developed  in  the  paper  [10].  The  solution  of  the 
system  allows  us  to  derive  the  time  dependencies  of 
the  basic  characteristics  of  the  film  growth  process: 

1 .  the  supersaturation  £(f)  is  as 

£(0  =  £r — p— * —  (14) 

l  +  -r*(/)^(7’(0) 

r 

2.  the  frequency  of  nuclei  formation  I(t )  is  as 


lit)  =  /(£,) 


exp  [-T\t)0t(Tm 

i+V(o^(no) 


3.  the  size  distribution  function  /(»;/ )  is  as 


/(«;0  =  ■ 


k  + : 


exp 


-(no 


k 


Pin) 

atA 


x 

<  X 

■y*tm  o- 


atA,  J 


(15) 


(16) 


where  the  auxiliary  function  is  p(n)  - 


Fig.  2.  The  time  dependence  of  supersaturation  on  the  surface 
at  different  values  of  the  incident  light  intensity:  curve  1  — 
102"  quantum- m'2 -s  ’;  curve  2  —  10w  quantum- m'2- sA 3  — 
1018  quantum- m'2-*'1;  l  —  time  (s),  ^(/)  —  supersaturation 
(non-dimensional  quantity). 


The  formulae  (14-16)  contain  the  parameter  T  -  quan¬ 
tity,  which  is  inverse  with  respect  to  the  expansion  small 
parameter.  This  parameter  has  the  very  important  prop¬ 
erty:  F  ~  ncr.  All  other  symbols  in  these  formulae 
correspond  to  those  in  the  origin  paper  [10]. 

4.  RESULTS  OF  CALCULATIONS  AND 
CONCLUSIONS 

Assuming  that  the  growing  film  is  determined  by  the 
following  constants:  Ig  =  1020  quantum-m  '-s ‘',q  =  0.8, 
ot  =  50  nr',  l  =  2-10'4  m  [2-5],  D  =  10'9  m2-s'',\  =  10‘30 
m\  c*  —  1020  aloms  nr2,  a  =  1  J-nr2,  k  =  1  ,a  =  1 014,  K  = 
1032/tr=  1 0‘2  s  [11],  SxMf=  10-4  m1,  T-  298  K  . 

Thus  the  value  of  the  steady  flux  (Eq.  6)  is  2-  1018 
atoms- m  2-s~' .  ncr  «  T  =  53.  That  criterium  proves  that 
the  solution  obtained  by  the  small  parameter  method 
sufficiently  describes  the  kinetics  of  the  film  growth. 
The  time  lag  value  is  /.  =  1.8 TO'6  s. 

The  results  of  calculations  are  presented  in  Figs. 
2-5. 

It  is  seen  from  Fig.  4,  the  size  distribution  takes  the 
form  of  a  bell.  It  is  the  form  which  is  predicted  by  the 
theory  [1,7]. 

We  have  calculated  the  time  dependencies  of  the 
size  distribution  function,  the  frequency  of  nucleation 
and  the  supersaturation  on  the  surface  at  three  differ¬ 
ent  values  of  the  incident  light  intensity  of  1020,  1019 
and  1018  quantum- nr' -s’' .  The  variation  of  light  inten¬ 
sity,  in  great  range,  keeps  the  kinetics  of  the  process  of 
film  deposition  out  of  the  drastic  changes. 


m  + 1 
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Fig.  3.  The  time  dependence  of  nuclcation  frequency  on  the 
surface  at  different  values  of  the  incident  light  intensity:  curve 
1  —  10211  quanlumrn^-s^;  curve  2  —  1019  quantum- m'2-sA]  3 
—  1018  quantum-m'2-s''\  t  —  time  (s),  /(/)  —  nuclcation 
frequency  (s'1). 


Fig.  5.  The  size  distribution  function  at  the  irradiation  time 
of  40  min  and  the  incident  light  intensity:  curve  1  —  1020 
quantum-m’2-s'1-,  curve  2  —  1019  quantum-m'2-s']\  1018 
quantum-m '2,s'';  n  —  the  number  of  nuclei  (non-dimensional 
quantity),  /(n)  —  the  size  distribution  function  (non- 
dimensional  quantity). 


Fig.  4.  The  size  distribution  function  at  the  incident  light 
intensity  of  /(|  =  10z"  quanlum-m'2-sA  and  the  irradiation  time: 
curve  1  —  5  min;  curve  2  —  10  min;  curve  3  —  40  min;  n  — 
the  number  of  nuclei  (non-dimensional  quantity),  f(n )  —  the 
size  distribution  function  (non-dimensional  quantity). 
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Abstract.  A  theoretical  model  is  constructed  that  describes  diffusion  processes  enhanced  by  the  elastic 
interaction  of  the  diffusing  species  and  amorphous  grain  boundaries  in  nanocrystalline  and  polycrystalline 
solids.  A  partial  differential  equation  governing  evolution  of  the  spatially  inhomogeneous  concentration  of 
the  diffusing  species  near  amorphous  grain  boundaries  is  numerically  solved.  It  is  shown  that  the  enhanced 
diffusion  essentially  influences  the  growth  of  amorphous  grain  boundaries  in  solids  under  irradiation  and 
thermal  treatment  as  well  in  sintered  ceramics  and  is  capable  of  contributing  to  the  anomalously  fast  diffusion 
in  nanocrystallinc  solids. 


1.  INTRODUCTION 

The  process  of  the  diffusion,  where  one  species 
rearranges  its  configuration  in  a  lattice  (a  solid  in  this 
case)  in  response  to  an  external  force,  is  an  attempt  to 
approach  equilibrium.  When  a  concentration  gradient 
exists  as  the  only  driving  force,  this  irreversible  pro¬ 
cess  can  be  concisely  formulated  leading  to  a  second 
order  parabolic  partial  differential  equation  (PDE).  The 
solution  of  this  PDE  has  been  determined  for  many 
types  of  boundary  and  initial  conditions  and  geometries. 

If  additional  driving  forces  are  present  due  to  elec¬ 
trical,  magnetic,  thermal  or  elastic  fields,  the  appro¬ 
priate  PDE  can  be  formulated  subject  to  simplifying 
assumptions.  The  question  now  becomes,  whether  this 
additional  driving  force  will  enhance  or  retard  the  flow 
of  diffusing  species  thorough  the  lattice.  In  particular, 
we  are  interested  in  the  effect  of  an  amorphous  inter¬ 
granular  boundary  dislocation  stress  field  has  on  the 
diffusion  of  various  species  in  nanocrystalline  and  poly¬ 
crystalline  solids.  Such  boundaries  frequently  exist  in 
sintered  nanostructured  and  polycrystalline  ceramics 
[1-5]  as  well  in  monophase  and  composite  solids  un¬ 
der  irradiation  and  thermal  treatment,  where  diffusion- 
and  irradiation-induced  amorphization  processes  of¬ 
ten  start  to  occur  at  grain  boundaries  [6,7].  Amorphous 
grain  boundaries  in  sintered  ceramics  are  characterized 
by  a  chemical  composition  that  is  different  from  that 
of  adjacent  crystalline  grains  as  well  by  a  nano-scaled 


thickness  whose  value  is  essentially  constant  within  a 
ceramic  sample.  These  grain  boundaries  are  dependent 
on  only  the  chemical  composition  and  insensitive  to 
misorientation  parameters  of  the  boundaries  (see  [1,8] 
and  references  therein).  Amorphous  grain  boundaries 
in  monophase  and  composite  solids  under  irradiation 
and  thermal  treatment  play  the  role  as  nuclei  for 
enhanced  diffusion-  and  irradiation-induced 
amorphization,  where  the  boundaries  transform  into 
extended  amorphous  layers  with  the  meso-  and  macro¬ 
scale  thickness  growing  parallel  with  the  diffusion  and 
irradiation  processes  [6,7].  Recently,  a  model  [9,10] 
has  been  suggested  describing  the  role  of  elastic  dis¬ 
tortions  in  the  formation  of  amorphous  grain  bound¬ 
aries.  This  model,  however,  does  not  quantitatively  take 
into  account  the  effect  of  diffusion  on  evolution  of 
amorphous  grain  boundaries  after  they  have  been 
formed.  At  the  same  time,  this  effect  is  important  for 
an  adequate  description  of  both  the  structural  features 
and  the  properties  of  nanocrystalline  and  polycrystal¬ 
line  solids  containing  amorphous  grain  boundaries. 

The  model  presented  here  is  an  initial  attempt  to 
understand  the  role  of  enhanced  diffusion  in  evolution 
of  amorphous  grain  boundaries  in  (i)  sintered  ceram¬ 
ics  and  (ii)  monophase  and  composite  solids  under 
irradiation  and  thermal  treatment.  The  model  also  is 
concerned  with  the  contribution  of  amorphous  grain 
boundaries  to  enhanced  diffusion  in  nanostructured 
solids  (often  exhibiting  the  outstanding  diffusional 
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properties  [1 1,12]),  which  is  related  to  elastic  interac¬ 
tion  of  amorphous  grain  boundaries  and  the  diffusing 
species.  As  such,  we  have  assumed  that  the  crystalline 
grains  are  separated  by  an  amorphous  grain  boundary 
layer  that  provides  the  source  of  the  diffusing  species. 
To  retain  the  flavor  of  dislocation  model  of  a  simple 
tilt  boundary,  we  will  consider  that  a  double  wall  of 
dislocations  separate  the  amorphous  layer.  We  further 
assume  that  dislocations  do  not  interact  with  each  other. 

We  will  discuss  the  development  of  the  appropriate 
PDE  and  the  elastic  interaction  from  a  double  wall  of 
dislocations.  This  will  lead  to  the  formulation  of  the 
interaction  energy  that  will  be  used  in  the  PDE  whose 
solution  is  determined  numerically.  Finally,  we  will 
discuss  the  results  of  the  numerical  analysis. 


2.  MATHEMATICAL  FORMULATION 

The  flux  or  current  density  of  a  diffusing  species  in  the 
presence  of  a  concentration  gradient  is  given  as, 

j  =  -DVC,  0) 

where  D  is  the  diffusivity  and  C  is  the  concentration. 
If  there  is  an  elastic  interaction  between  the  diffusing 
species  and  a  stress  field  in  a  homogenous,  linear  and 
isotropic  solid,  Eq.  1  is  modified  as  [13] 


DVC  + 


(2) 


where  e  is  the  interaction  energy,  k  is  Boltzmann's 
constant  and  T  is  the  absolute  temperature  of  the  diffu¬ 
sion  process.  It  is  assumed  that  the  mobility  p,  caused 
by  the  force  from  the  strain  gradient  is  the  same  as  in 
the  chemical  potential  force, 


D 


M  = 


kT' 


(3) 


It  will  be  assumed  that  the  diffusivity  is  constant  and 
independent  of  both  concentration  and  strain. 

Using  the  continuity  equation,  Eq.  2  can  be  recast  as 


—  =  d(Vc  +  —  VC- 

a  V  kT 


(4) 


since 


due  to  the  absence  of  any  internal  body  forces.  Eq.  4 
along  with  the  appropriate  initial  and  boundary 
conditions  will  be  solved  to  determine  how  the 
concentration  is  affected  by  the  elastic  interaction. 


A  formulation  for  the  interaction  energy  based  upon 
a  grain  boundary  model  is  required.  For  nanocrystalline 
and  polycrystalline  materials,  an  appropriate  model  is 
the  one  proposed  by  Kolesnikova  et  al.  [9]  and  Ovid'ko 
and  Reizis  [10]  where  an  amorphous  grain  boundary 
layer  occurs  as  described  below  between  the  nanosized 
grains. 

3.  DOUBLE  WALL  MODEL 

In  the  framework  of  dislocation-structural-unit  mod¬ 
els  [14,15],  a  conventional  tilt  grain  boundary  in  a 
crystal  can  be  treated  as  a  periodic  wall  array  of  edge 
grain  boundary  dislocations  whose  spatial  positions  and 
Burgers  vectors  are  strictly  determined  by  parameters 
(misorientation,  etc.)  of  the  boundary.  In  this  model 
representation,  the  elastic  fields  and  the  elastic  energy 
density  of  a  tilt  boundary  are  those  of  the  corresponding 
array  of  grain  boundary  dislocations  [16,17],  Follow¬ 
ing  [9,10],  transformation  of  a  conventional  tilt  grain 
boundary  into  an  amorphous  grain  boundary  is  associ¬ 
ated  with  a  re-arrangement  of  its  grain  boundary  dis¬ 
location  ensemble.  In  doing  so,  due  to  the  fact  that  the 
amorphous  disordered  structures  do  not  impose  strict 
limitations  on  both  spatial  positions  and  Burgers  vec¬ 
tors  of  grain  boundary  dislocations,  such  dislocations 
in  a  new,  amorphous  boundary  are  arranged  in  the  way 
decreasing  their  elastic  energy  density.  The  discussed 
decrease  of  the  elastic  energy  density  causes  the  basic 
driving  force  F  t  for  the  transformation  of  a 
conventional  tilt  boundary  into  an  amorphous  bound¬ 
ary.  This  force  is  calculated  in  the  model  [9,10]  which 
to  a  first  approximation  describes  the  amorphous 
boundary  to  be  of  finite  width,  X,  with  dislocation  arrays 
at  each  of  the  interfaces  between  amorphous  and  the 
crystalline  grain  as  shown  in  Fig.  1.  Following  the 
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Fig.  1.  (a)  Single  dislocation  array  and  (b)  double  wall  with 
amorphous  grain  boundary  of  width  X. 
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calculations  of  references  [9]  y.d  [10],  the  driving  force 
F  is  large  at  X  <  /7,  and  =  0  at  X>  h„  where  h ,  is  the 
period  of  the  dislocation  walls  (Fig.  1).  The  elastic  in¬ 
teraction  between  the  two  dislocation  walls  exists  in 
only  the  range  of  X  <  /?,;  the  walls,  in  general,  do  not 
interact  elastically  at  X  >  h,.  As  a  corollary,  evolution 
of  the  amorphous  grain  boundary  at  X  -  h2  is  caused  by 
the  thermodynamic  force  F'  related  to  other  factors 
other  than  the  decrease  of  the  elastic  energy  density. 
More  precisely,  the  force  F’  is  as  follows  [10]: 


F'=  -ea_c(X)-(2i<  +  X) 


dsn  c 
dX  ' 


(6) 


where  za<  denotes  the  difference  between  the  free  energy 
densities  of  the  amorphous  and  crystalline  phases,  and 
k  the  characteristic  spatial  scale  of  structural 
inhomogenities  in  the  amorphous  phase.  The  term  e;t 
strongly  depends  on  the  chemical  composition  which, 
in  general,  is  spatially  inhomogeneous  in  the  vicinity 
of  the  amorphous  boundary  and  can  evolve  in  time  due 
to  diffusion.  In  these  circumstances,  the  F'-induced  evo¬ 
lution  of  the  amorphous  grain  boundary  is  sensitive  to 
the  diffusion  processes  in  its  vicinity.  On  the  other  hand, 
the  stress  field  of  the  boundary  influences  the  diffusion 
processes  and,  as  a  corollary,  evolution  of  the  bound¬ 
ary.  The  elastic  interaction  between  the  diffusing  spe¬ 
cies  and  the  stress  field  is  related  to  the  sum  of  the 
normal  stresses.  For  a  double  wall  of  non-interacting 
dislocation  infinite  arrays  a  distance  X  apart  the  sum 
of  the  normal  stresses  is, 


a,.  = 


Gb2{\  +  v) 
2(1  +  v)h1 


co(  2n^) 


where 


(7) 


®(?7,£> 


sin(2  jt£) 

- + 

X 

cosh  {2  ^-(^  - — )}  -  cos(2;r£) 

2 

sin(2/z2f) 


cosh{2/r(7/  +  — )}  -  cos(2/r£) 


(8) 


and  letting  q  =  x/h2,  =  )’/ h2  and  X'  =  XI h  for  the 
coordinate  system  given  in  Fig.  1  where  the  disloca¬ 
tion  walls  are  symmetrical  located  at  x  =  ±X/2-  The 
shear  modulus  and  Poisson's  ratio  are  denoted  by  G 
and  v,  respectively. 


4.  INTERACTION  ENERGY 

We  assume  the  volume  of  the  diffusing  and  host  spe¬ 
cies  are  different  and  thus  giving  rise  to  an  elastic  in¬ 
teraction.  On  the  average,  letAF  be  the  volume  change 


and  this  volume  difference  interacts  with  a  homogenous 
applied  stress  field.  It  can  be  shown  that  the  elastic 
interaction  is  given  as 

e  =  (AV)Q,  (9) 

where  ©  =-o«/  3,  the  hydrostatic  component  of  the  stress 
field.  This  interaction  will  alter  the  total  energy  of  the 
system  depending  upon  the  sign  of  the  product. 

In  Fig.  2a  is  plotted  the  contours  of  constant 
normalized  interaction  energy  (i.e.,  co(p,^))  for  a  strip 
|£,|<0.5.  As  expected,  a  large  energy  gradient  exists  near 
the  dislocation,  E,-0  due  to  the  elastic  singularity  of 
the  dislocation  core.  Depending  upon  the  volume 
change,  this  elastic  interaction  will  enhance  diffusion. 
For  illustrative  purposes,  it  will  be  assumed  that  AK>0. 
In  Fig.  2a,  from  a  qualitative  viewpoint,  the  diffusing 
species  would  to  migrate  below  the  extra  dislocation 
half  plane  to  lower  their  total  energy.  A  spherical  point 
defect  interacting  with  a  single  edge  dislocation  also 
exhibits  this  behavior  [18]. 

The  force,/]  on  the  migrating  specie  is  determined 
from  -Vs  and  the  results  are  shown  in  Fig.  2b  for  the 
corresponding  interaction  energy  of  Fig.  2a.  The  dif¬ 
fusing  specie  will  tend  to  move  away  from  above  the 
dislocation  and  migrate  to  the  edge  of  the  extra  half 
plane  of  the  dislocation.  This  periodic  case  of  a  tilt 
boundary  is  different  from  that  for  an  isolated  edge 
dislocation  in  that  the  diffusing  specie  in  the  upper 
one  fourth  of  the  strip  diffuses  to  the  next  strip  above. 
This  is  only  a  qualitative  assessment  and  should  serve 
as  guide  when  a  concentration  gradient  is  not 
considered. 


5.  NUMERICAL  SOLUTION 


Before  solving  Eq.  4,  the  initial  and  boundary  conditions 
must  be  specified.  We  need  to  only  consider  one  strip,  |^| 
<  1/2,  since  we  assume  the  grain  boundary  to  be  mod¬ 
eled  by  an  infinite  array  of  dislocations  and  can  utilize  a 
periodic  boundary  condition,  C(q, -1/2,0  =  C(r|,  1/2,0 
for  r|  >  0.  We  will  assume  that  the  amorphous  layer  in 
between  the  nanosized  grains  will  be  an  undiminished 


source  of  the  diffusing  specie,  thus,  C(0,^,/)  =  1.0  for 
|^|  <  1/2.  Since  the  dislocation  stresses  from  a  regular 
array  are  short  range  in  nature,  it  is  reasonable  to  uti¬ 
lize  C(co,^,0=0  in  each  and  every  strip.  The  above  bound¬ 
ary  conditions  are  valid  for  t  >  0.  For  an  initial  condition, 
we  will  assume  that  each  strip  is  devoid  of  any  of  the 
diffusing  specie  at  t  =  0  leading  to  C(r),£,,0)  =  0. 

Eq.  4  is  normalized  for  computational  purposes  by 
rescaling  time  and  letting  the  material  coefficients  become 
DGb2{  1  +  v) 

(10) 


t  =  Dt  and  y  = 


2kt{\  -  v)h2 
The  normalized  computational  form  becomes 
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Fig.  2(a)  shows  the  normalized  interaction  energy,  (0.  Above 
the  dislocation,  CO>0  (solid)  while  below  CO<0  (dashed)  for 
A V>0.  The  area  to  the  left  of  the  dislocation  is  one  half  of  the 
amorphous  layer. 


—  -  V2C  +  /VC-Vcu  (11) 

dr 

with  the  boundary  conditions  for  t  >  0, 

C( 77-1/2,  r)  =  C(  >7,1/2,  r)  for  >7  >  0,  (12) 

C(0, £,  r)  =  1.0  for  |  £\  <  1/2 ,  (13) 

C{oo,4,t)  =  0  for  1 41  <  1/2 ,  (14) 

and  the  initial  condition  at  x  =  0, 

C(? j,%,  r)  =  0  for  i)  >  0  and  |£|<  1/2.  (15) 


An  alternating-direction-implicit  (ADI)  method  is  uti¬ 
lized  to  obtain  a  numerical  solution  of  Eq.  11.  Each 
time  step  is  divided  equally  in  half.  In  the  first  half,  a 
partial  solution  is  determined  by  reducing  the  PDE  to 
a  one-dimensional  problem  in  one  of  the  orthogonal 
coordinates  while  keeping  the  other  fixed.  In  the  sec¬ 
ond  half,  the  PDE  is  now  solved  in  the  other  coordinate 
direction  using  the  partial  results  obtained  in  the  first 
half  time  step.  This  algorithm  has  the  advantage  of 
requiring  a  tridiagonal  equation  solver  and  is  relatively 
stable  (unconditionally  if  the  coefficients  of  the  PDE 
are  constant).  The  q-direction  partial  solution  requires 
a  slight  modification  to  the  tridiagonal  solver  to  ac¬ 
count  for  periodic  nature  of  the  boundary  conditions 
but  it  still  retains  the  usual  banded  form.  We  have 
chosen  to  do  the  analysis  first  in  the  q-direction  since 
the  elastic  fields  are  short-ranged  in 

The  finite  difference  equations  are  based  upon  the 
usual  five  atom  centered-computational  molecule  to 
approximate  the  Laplacian  (spatial  coordinates)  and  a 
two  point  approximation  for  the  temporal  first  deriva¬ 
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Fig.  2(b)  indicates  the  flow  of  diffusant  (A V>0) 
corresponding  to  interaction  energy  shown  in  Fig.  2a  as 
determined  from  -Vo).  For  clarity,  each  vector  is  normalized 
to  one.  The  forces  away  from  the  dislocation  arc  greatly 
attenuated.  The  area  to  the  left  of  the  dislocation  is  one  half  of 
the  amorphous  layer. 


tive.  Using  equal  increments  for  the  spatial  coordinates, 
A£,  =  Aq  and  utilizing  a  normalized  time  increment 
At,  and  rescaling  y  for  computation  clarity, 

/=4y(A£),  (16) 

Eq.  12  can  be  approximated  by  finite  difference  form 
amenable  for  numerical  computation. 

For  the  actual  computation,  a  strip  |£  |  <  1/2  high 
and  0  <  q  <  3  was  used.  The  grid  for  this  rectangle  was 
divided  into  100  and  300  equal  increments  for  the  £- 
and  q-directions,  respectively.  The  singularity  of  stress 
field  was  avoided  by  placing  the  dislocation  in  between 
the  grid  points.  For  each  time  step,  the  mass  accumu¬ 
lation,  M{t),  was  determined  by  numerically  integrat¬ 
ing  the  normalized  concentration  over  the  computa¬ 
tional  rectangle. 

To  check  the  accuracy  of  the  numerical  scheme,  an 
initial  run  was  made  with  the  interaction  set  to  zero. 
This  case  corresponds  to  diffusion  between  two  planes 
(one-dimensional  problem).  These  results  agreed  with 
the  analytical  solution  lending  confidence  to  the 
numerical  algorithm.  Runs  were  made  to  explore  the 
effect  of  the  elastic  interaction  between  grain  bound¬ 
ary  dislocations  and  diffusing  species.  First  was  to  de¬ 
termine  the  role  of  the  “strength”  of  the  interaction  as 
manifested  through  y*  and  second  was  to  determine 
the  effect  of  spacing  between  the  dislocations  (width 
of  amorphous  layer)  in  the  wall,  A,  on  the  diffusion 
process. 

A  typical  run  for  values  of  y*  =  3  and  A  =  1  is  shown 
in  Fig.  3  for  various  normalized  times.  The  concentra¬ 
tion  is  larger  above  the  dislocation  since  the  flow  lines 
from  Fig.  2a  indicate  diffusion  away  from  this  area 
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Fig.  3.  Normalized  concentration  contours  for  various  normalized  times  with  and  without  elastic  interaction.  The  values  of  y  — 
3,  X  =  1  and  y*  =  0,  X  =  1  were  used.  The  normalized  contours  have  the  same  range  of  values  for  each  time  slice.  The  amorphous 
layer  and  other  dislocation  wall  are  not  shown.  The  zero  interaction  case  is  the  smaller  segment  with  the  uniform  contours. 


while  the  concentration  is  greatly  diminished  below 
the  dislocation.  The  flow  pattern  generated  from  the 
static  case,  Fig.  2b  indicates  flow  away  from  the  inter¬ 
face  above  the  dislocation  and  flow  to  the  interface 
below  the  dislocation  due  to  the  elastic  interaction.  The 
concentration  gradient  would  prefer  to  have  a  uniform 
flow  (no  x  variation)  but  cannot  overcome  the  elastic 
interaction  particularly  near  the  singularity  at  the  dis¬ 
location.  The  diffusing  specie  accumulates  below  the 
extra  half  plane  (recall,  AF  >  0)  and  diffuses  very 
rapidly  to  this  area  due  to  the  mathematical  singular¬ 
ity.  Some  the  diffusant  originates  from  the  strip  below 
as  can  be  seen  from  the  Fig.  2b.  However  the  diffusing 
specie  attempts  decreases  its  concentration  in  the  region 
above  the  extra  half  plane  but  cannot  do  so  due  to  its 
source  at  the  amorphous-crystalline  interface.  Thus,  it 
flows  over  a  larger  area  at  nearly  the  maximum 
concentration  in  response  to  the  elastic  interaction 


enhancing  the  diffusion  process.  A  fraction  of  the 
diffusant  also  flows  to  the  strip  above  this  one. 

At  later  times,  the  diffusing  specie  is  not  effected 
by  the  elastic  interaction  due  to  the  short-ranged  effect 
of  the  grain  boundary  dislocation  elastic  fields.  Away 
from  the  dislocation,  the  concentration  profile  becomes 
flat  indicating  a  decreased  influence  ofthc  elastic  fields. 
Note  however,  at  even  these  later  times,  the  elastic  in¬ 
teraction  near  the  dislocation  remains  large  as  can  be 
seen  in  Fig.  3.  (In  Fig.  3,  the  iso-contour  lines  have  the 
same  range  of  values  for  the  various  time  slices.)  For 
comparison,  the  zero  interaction  case,  y*  =  0,  is  shown 
in  Fig.  3  in  the  small  segment.  As  expected,  the  iso¬ 
contours  are  planar  for  this  one-dimensional  case.  As 
can  been  seen  from  this  figure,  the  iso-contours  lines 
when  y*  =  0  lag  behind  the  interaction  cases  and  hence 
affect  the  mass  accumulation  of  the  diffusant. 

The  presence  of  the  stress  field  from  a  grain  bound¬ 
ary  dislocation  affects  the  mass  accumulation.  In  Fig.  4 
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Fig.  4.  Normalized  mass  accumulation  for  several  values  of  y 
and  A.. 


r,  normalized  time 

Fig.  5.  Effect  of  width  of  dislocation  wall  spacing  on  mass 
accumulation. 


is  plotted  the  mass  accumulation  for  several 
combinations  of  y*  =  0,  1,3  and  5  and  X  =  0.5  and  1. 
There  are  several  salient  features  to  note.  First,  as  ex¬ 
pected,  the  amount  of  mass  diffusion  is  increases  with 
increasing  elastic  interaction,  i.e.,  larger  y*.  Second 
there  is  very  little  discernable  difference  between  X  = 
0.5  and  X  =  1  indicating  that  the  short-range  Field  ef¬ 
fects  are  limited  at  these  widths.  Third,  there  is  an 
enhanced  diffusion  over  the  base  line  case  of  y*  =  0  of 
the  order  of  20%  at  the  later  times.  Given  the  little 
effect  of  1  on  the  mass  diffusion,  the  baseline  results 
for  y*  =  0.5  or  1  should  be  about  the  same.  For  that 
reason,  we  have  shown  only  X  =  1  for  y*  =  0. 

The  width  effect,  e.g.,  the  spacing  between  the  dis¬ 
location  walls  has  a  marginal  effect  as  the  spacing,  X, 
becomes  smaller.  In  Fig.  5,  X  is  varied  from  0.1,  0.5 
and  1  for  two  values  of  y  =  1  and  5.  There  is  a  slight 
increase  in  mass  accumulation  as  the  width  is  decreased 
to  a  value  of  0. 1  (recall  this  means  that  the  spacing  of 
the  dislocation  in  each  wall  or  array  is  10  time  larger 
that  the  separation  distance  between  the  walls).  This 
analysis  should  be  treated  only  as  an  indication  of  ef¬ 
fect  of  the  elastic  interaction  because  this  analysis  uses 
the  assumption  of  non-interacting  dislocations. 

4.  DISCUSSION 

We  discuss  the  evolution  of  amorphous  grain  bound¬ 
aries  in  (i)  sintered  ceramics  and  (ii)  solids  under 
irradiation  and  thermal  treatment  based  on  the  role  of 
enhanced  diffusion.  To  do  so,  we  assume  that  the  spe¬ 
cies  diffusing  from  a  pre-existent  amorphous  grain 
boundary  into  adjacent  crystalline  regions  are  capable 
of  inducing  the  amorphization  of  these  regions.  This 
corresponds  to  experimentally  observed  situations  (i) 
[1-5]  and  (ii)  [6,7]  where  diffusion  mixing  causes 
growth  of  the  amorphous  phase  at  grain  boundaries. 
With  the  above  taken  into  account,  the  following  sce¬ 
nario  of  the  evolution  is  suggested  which  takes  into 


consideration  the  role  of  enhanced  diffusion  in 
amorphizing  transformations  at  grain  boundaries.  Pen¬ 
etration  of  the  “amorphizing”  diffusing  species  from 
an  amorphous  grain  boundary  to  adjacent  crystalline 
regions  occurs  in  a  spatially  inhomogeneous  way  (see 
previous  sections  and  Fig.  6a)  and  leads  to  change  of 
e  in  the  regions  with  a  high  density  of  the 
amorphizing  diffusing  species.  This  induces 
amorphization  (crystal-to-glass  transition)  in  these 
regions  (Fig.  6b).  The  spatially  inhomogeneous  growth 
of  the  amorphous  phase  (Fig.  6b)  forces  crystalline- 
amorphous  interface  to  be  curved  and,  therefore,  results 
in  an  increase  of  its  contribution  to  the  total  energy 
density  of  the  amorphous  boundary.  The  crystalline- 
amorphous  interface  can  straighten  out  and,  as  a 
corollary,  decrease  the  total  energy  density  in  the  fol¬ 
lowing  two-step  way.  First,  grain  boundary  dislocations, 
due  to  absence  of  structure-imposed  limitations  on  their 
spatial  positions,  move  to  new  positions  (Fig.  6c).  The 
dislocations  located  in  new  positions  create  a  new  stress 


Fig.  6.  Increased  width  of  layer,  8k,  by  amorphization  due  to 
enhanced  diffusion. 
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field  distribution  which,  in  its  turn,  enhances  penetra¬ 
tion  of  the  amorphizing  diffusing  species  from  the 
amorphous  grain  boundary  to  new  crystalline  regions 
located  between  the  “bugles"  of  the  boundary  (Fig.  6d). 
The  diffusion  is  followed  by  amorphization  of  these 
new  regions,  in  which  case  the  crystal/glass  interface 
becomes  planar  (Fig.  6e). 

If  there  is  a  permanent  “supply”  of  the  amorphizing 
diffusing  species  to  the  amorphous  grain  boundary  due 
to  grain  boundary  diffusion  (which  quickly  transfers 
the  amorphizing  diffusing  species  from  some  remote 
permanent  source,  e.g.  a  second  phase  region  consisting 
of  chemical  elements  that  serve  as  the  amorphizing 
diffusing  species)  or  irradiation  (see  situation  (ii)),  then 
the  boundary  plays  the  role  as  a  permanent  source  of 
the  amorphizing  diffusing  species  that  penetrate  into 
the  adjacent  crystalline  grains.  In  these  circumstances, 
the  process  shown  in  Fig.  6  is  permanently  repeated 
resulting  in  diffusion-  or  irradiation-induced  perma¬ 
nent  growth  of  the  amorphous  grain  boundary.  This 
corresponds  to  experimentally  observed  [6,7]  situation 
(ii)  where  the  diffusion-  and  irradiation-induced  growth 
of  the  amorphous  phase  starts  occurring  at  grain  bound¬ 
aries,  which  then  gradually  extend  to  amorphous 
regions  of  macro-scale  extent. 

if  the  quantity  of  the  amorphizing  diffusing  species 
is  limited,  as  in  situation  (i)  with  sintered  ceramics, 
the  process  shown  in  Fig.  6  results  in  formation  of  an 
amorphous  boundary  of  finite  width.  The  boundary 
width  is  determined  only  by  the  sum  (finite)  quantity 
of  the  amorphizing  diffusing  in  a  sample.  Since  grain 
boundary  diffusion  is  faster  by  several  orders  than  bulk 
diffusion  (that  controls  the  growth  of  amorphous  grain 
boundaries;  see  Fig.  6),  basically,  the  same  quantities 
of  the  amorphizing  diffusing  species  are  distributed 
along  all  the  grain  boundaries  within  a  sample.  This 
explains  the  experimentally  observed  (see  [1,8]  and 
references  therein)  specific  feature  of  amorphous  grain 
boundaries  in  sintered  ceramics,  which  is  a  finite  bound¬ 
ary  width  whose  value  is  tentatively  the  same  for  all 
amorphous  grain  boundaries  within  a  ceramic  sample 
and  depends  on  only  the  chemical  composition. 

The  diffusion  processes  enhanced  due  to  the  elastic 
interaction  ofthe  diffusing  species  and  amorphous  grain 
boundaries  represent  an  example  of  diffusion  processes 
dependent  on  stress  field  distributions  in  polycrystal¬ 
line  and  nanocrystalline  solids.  The  contribution  of  such 
processes  to  the  macroscopic  diffusional  properties  is 
significant  in  solids  with  evolving  defect  structures  and, 
therefore,  evolving  stress  field  distributions.  In  particu¬ 
lar,  the  enhanced  diffusion  which  accompanies  trans¬ 
formations  of  amorphous  grain  boundaries  (Fig.  6)  is 
capable  of  essentially  contributing  to  the  anomalously 
fast  diffusion  in  nanocrystalline  solids  where  evolving 


(amorphous  and  conventional)  grain  boundaries  are 
inherent  structural  elements. 

4.  CONCLUSIONS 

As  an  initial  to  attempt  to  under  the  role  of  grain  bound¬ 
aries  and  their  effect  in  nanosized  materials,  we  have 
modeled  a  simple  case  of  two  non-interacting  grain 
boundary  arrays  or  walls  separating  an  amorphous  layer 
that  provides  a  source  for  the  diffusing  species.  A  simple 
elastic  interaction  is  considered  and  inserted  into  the 
diffusion  equation  that  is  solved  for  the  time  depen¬ 
dent  behavior  of  the  concentration.  The  mass  accumu¬ 
lation  is  enhanced  due  to  the  elastic  interaction  and 
that  until  the  spacing  becomes  small,  the  second  dislo¬ 
cation  array  or  wall  has  a  minimal  effect.  And  lastly, 
the  static  information  derived  from  the  interaction 
energy  and  its  gradient,  e.g.  force  on  the  diffusing 
specie,  will  only  yield  qualitative  information  and  that 
time  dependent  analysis  is  required.  The  enhanced  dif¬ 
fusion  essentially  influences  the  growth  of  amorphous 
grain  boundaries  in  solids  under  irradiation  and  ther¬ 
mal  treatment  as  well  in  sintered  ceramics  and  is 
capable  of  contributing  to  the  anomalously  fast  diffu¬ 
sion  in  nanocrystalline  materials. 
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Abstract.Thc  superconducting  transport  properties  (/-If  I-B  and  /  T)  were  correlated  with  a  systematic 
characterization  of  the  grain  boundary  microstructurc,  via  TEM  observation  and  local  chemical  analysis  in 
STEM  by  EDX  on  YBa2Cu1sOv  (D  samples).  YBa,Cu3  P;/Ags  (S  samples)  and  DyBa2Cu3  xOy/lwt%  Pt 
(0<x<0.4)  ceramics.  It  was  shown  that  in  all  samples  studied  the  critical  current  /  is  controlled  by  weak 
links  at  grain  boundaries.  The  measurements  of  1(7)  indicated  that  in  copper-deficient  yttrium  ceramics, 
Ag-doping  leads  to  change  of  weak  link  character  from  supcrconductor-insulator-superconductor  (SIS)  to 
superconductor-normal  mctal-supcrconductor  (SNS)  type.  This  SNS  behaviour  seems  to  be  a  result  of 
percolation  path  of  supercurrent  through  clean  boundaries  with  extremely  narrow  range  (=1  nm)  of  Ag 
segregation  on  them.  Specific  Ag  precipitates  of  2-5  nm  size  on  or  near  clean  boundaries  were  found  in  S 
samples  for  x=0.4.  This  is  most  likely  to  be  a  reason  for  the  improvement  of  /  in  these  samples  by  a  factor 
of  3.  In  copper-deficient  dysprosium  ceramics  doped  with  Pt,  Pt-doping  leads  to  a  substantial  increase  in  / 
and  pronounced  'fish-tail’  effect  in  magnetic  field  at  T<11  K  in  non-stoichiometric  composition  with  x=0.2. 
Better  behaviour  of  /  in  magnetic  fields  in  these  samples  may  be  associated  with  the  substantial  increase  in 
the  proportion  of  clean  boundaries  and  decrease  in  twin  spacing. 


1.  INTRODUCTION 

In  the  high-7)  granular  superconductors,  transport 
properties  are  mainly  controlled  by  the  grain  boundary 
microstructure  unlike  T  which  is  determined  by  the 
crystal  structure  and  oxygen  content.  Depending  on 
crystallography  and  structure  of  boundaries,  they  can 
be  favourable  for  the  passage  of  supercurrent  or  can  be 
weak  links  (reduced  /  regions)  in  the  superconducting 
current  path  of  a  ceramic  sample.  Recently  many  stud¬ 
ies  have  been  directed  at  investigating  the  dissipation 
at  a  separate  boundary  of  different  type  [1],  Although 
studies  of  separate  boundary  behaviour  highly 
contribute  to  understanding  transport  in  superconduct¬ 
ing  polycrystalline  materials,  they  are  rather  idealistic 
models  of  true  polycrystalline  materials  which  contain 
many  boundaries  of  different  types,  in  parallel  and  se¬ 
ries,  simultaneously.  So,  along  with  the  study  of  a  sepa¬ 
rate  boundary,  it  is  also  very  important  to  investigate 
the  relation  between  transport  properties  and  grain 
boundary  behaviour  as  a  network  of  boundaries. 


The  main  aim  of  this  work  was  to  find  possible 
correlation  between  microstructure  features  and  super¬ 
conducting  transport  properties.  One  way  for  changing 
intergrain  critical  current  density  of  polycrystalline 
materials  is  doping  with  different  elements.  Small  de¬ 
viation  from  the  stoichiometric  composition  can  also 
lead  to  a  change  in  structure  and  distribution  of  bound¬ 
aries  in  a  ceramic  sample. 

So,  in  this  work  the  superconducting  transport  prop¬ 
erties  {I-V,  I-B  and  /  -  7)  were  correlated  with  a  sys¬ 
tematic  characterization  of  the  grain  boundary  micro¬ 
structure,  via  TEM  observation  and  local  chemical 
analysis  in  STEM  by  EDX.  The  following  supercon¬ 
ducting  ceramics  were  studied. 

1 .  YBa2Cu3  Oy  (0<x<0.4)  -  Cu-deficient  yttrium  ce¬ 
ramics  (D  samples);  YBa2Cu3_xO  /Agx  (0<x<0.4)  ce¬ 
ramics  with  silver  incorporated  in  amounts  equal 
to  the  copper  deficiency  (S  samples). 

2.  DyBa,Cu3  xOy/l wt%  Pt  (0<x<0.4)  -  Cu-deficient 
dysprosium  ceramics  doped  with  1  wt%  Pt. 
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2.  EXPERIMENTAL 

The  YBCO  ceramic  system  was  fabricated  via  the  modi¬ 
fied  citrate  gel  process  [2].  The  sol-gel  method  allows 
one  to  achieve  highly  homogeneous  mixing  of  cations 
on  an  atomic  scale  and  produce  high-quality  powder. 
Two  types  of  YBCO  samples  were  prepared.  In  samples 
of  the  first  type  YBa2Cu3  Oy  (0<x<0.4)  (D  type)  Cu- 
deficient  compositions  were  obtained.  In  samples  of  the 
second  type,  YBa2Cu,  xOy/Agx  (0<x<0.4)  (substituted  S 
type),  Ag  was  considered  as  a  substituent  for  Cu. 

The  copper  deficient  ceramics,  DyBa2Cu3  xOy,  doped 
by  lwt.%  Pt,  were  prepared  by  the  usual  solid  state 
reaction.  For  comparison  two  compositions  with  x=0  and 
x=0.2  of  the  copper  deficient  ceramics  DyBa2Cu3  xOv 
without  Pt  were  also  fabricated  using  the  same  method. 

The  phase  were  characterised  by  powder  X-ray  dif¬ 
fraction  (XRD)  using  a  Philips  automated  X-ray 
diffractometer.  Microstructure  and  chemical 
compositions  were  analysed  with  a  Jeol  2000  FX  TEM 
equipped  with  a  Link  EXL  X-ray  selective  analyser 
and  with  an  SEM  equipped  with  an  EDAX  X-ray  se¬ 
lective  analyser. 

The  superconducting  properties  were  studied  by 
measuring  the  resistivity-temperature  (R-T)  and 
current-voltage  (/-I7)  characteristics  using  the  four- 
probe  technique.  The  critical  current  /  was  determined 
at  various  temperatures  using  the  1  pV  nmr1  criterion, 
with  and  without  external  magnetic  field  B.  The  exter¬ 
nal  magnetic  field  ranged  from  0  to  75  G  and  was  per¬ 
pendicular  to  the  transport  current. 

3.  RESULTS  AND  DISCUSSION 

3.1  COPPER  DEFICIENT  YTTRIUM  SYSTEM 

WITH  AND  WITHOUT  AG-DOP1NG 

First  we  will  consider  123-yttrium  system  with  copper 
deficiency  with  and  without  Ag-doping.  According  to 
the  data  in  the  literature  [3-5]  distribution  of  Ag  in 
ceramics  and  its  influence  on  ./  are  subject  to  contra¬ 
dictory  conclusions  and  strongly  depend  on  composition 
of  ceramics,  processing  technology  and  thermal  treat¬ 
ment.  Our  ceramic  samples  were  prepared  via  the  modi¬ 
fied  citrate  gel  process  [2]  which  allows  one  to  achieve 
highly  homogeneous  mixing  of  cations  on  an  atomic 
scale  and  produce  high-quality  homogeneous  ceram¬ 
ics. 

EDX  analysis  in  TEM  showed  (Table  1)  that  the 
Ag  content  (x)  of  YBCO  grains  was  lower  than  the 
nominal  amount  of  Ag,  in  spite  of  the  sol-gel  method. 
Ag  mainly  precipitated  as  small  particles  (1-2  pm) 
homogeneously  scattered  throughout  the  sample. 


Table  1.  Average  concentration  of  silver  in  the  grains 
measured  by  EDX  in  TEM  and  STEM  for 
YBa2Cu3  xOy/Agx  for  all  doping  levels. 


Nomina 1  Ag  concentration  x 

Real  Ag  concentration 

0.1 

<  limit  of  detection 

0.2 

<0.02 

0.4 

0.03 

X-ray  diffraction  and  SEM  analysis  showed  that 
all  samples  had  an  good  orthorhombic  structure.  For 
samples  with  x>0.2,  along  with  Ag  phase  (S  samples) 
several  nonsuperconducting  phases  such  as  Y2BaCu05 
(211),  BaCu02  (Oil)  and  CuO  (D  and  S  samples)  were 
registrated.  The  amount  of  these  nonsuperconducting 
phases  (in  D  and  S  samples)  and  Ag-phase  (S  samples) 
increased  with  increasing  x.  Second  phases  such  as 
BaCuO,  and  CuO  were  mainly  located  at  grain  bound¬ 
aries  and  triple  junctions.  These  impurity  phases  were 
about  1-3  pm  in  size,  the  mean  distance  between  them 
was  comparable  with  respective  grain  size.  The  aver¬ 
age  grain  size  dramatically  decreased  with  x  (Fig.  1). 
The  nonsuperconducting  phases  and  Ag  precipitates 
appear  to  stop  grain  growth  during  sintering  resulting 
in  fme-grain  microstructure  for  x>0.2. 

For  all  samples  of  both  type  (D  and  S)  the  critical 
temperature  was  practically  the  same  (7)=91-92  K), 
but  the  transport  properties  were  quite  different.  Fig.  2 
shows  the  measured  values  of  critical  current  density 
7  at  77  K  versus  x.  For  S  samples,  J  increased  only 
slightly  with  increase  in  Ag  content  from  x=0  to  x=0.2, 
whereas  J  increased  more  than  threefold  (from  .7=50 
A/cm2  for  x=0  to  7  =  170  A/cm2)  for  x=0.4.  For  D 
samples  with  x>0.2,  7  is  =  twice  that  obtained  for  the 
stoichiometric  composition. 

It  should  be  specially  noted  that  in  samples  of  both 
type  (D  and  S),  the  obtained  /- V  curves  at  77  K  exhib¬ 
ited  a  power-low  behaviour.  /  was  strongly  depressed 


Fig.  1.  Dependence  of  the  average  grain  size  of  YBa2CuJxOy 
(D)  and  YBa2Cu3  Oy/Agx  (S)  samples  on  x. 
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Fig.  2.  Variation  of  the  critical  current  with  x  for  D  and  S 
samples. 

by  low  magnetic  fields.  Such  behaviour  is  characteristic 
of  granular  materials  which  can  be  described  as  a 
collection  of  superconducting  grains  connected  by  weak 
links  [6], 

Thus,  in  all  samples  /.  was  controlled  by  grain 
boundary  weak  links  which  were  characterized  from 
7  (7)  curves.  It  is  known  that  for  a  Josephson  tunnel 
(SIS)  junction  the  dependence  of  Ic  on  temperature 
agrees  with  the  equation  7.=const  (1-7/7)”)  for  tempera¬ 
ture  close  to  the  superconducting  transition  tempera¬ 
ture  T'  of  the  intergrain  junction  [7].  For  the  SNS- 
type  proximity  junctions,  close  to  T'  the  7t(7)  depen¬ 
dence  can  be  written  in  the  form  /  =const  (1-777)*)2 
[8],  To  identify  the  character  of  weak  links  we  tried  to 
find  the  dependence  of  /.  on  temperature  for  T close  to 
the  critical  temperature  Tc  of  a  sample.  The  obtained 
dependence  of  7  on  (1-777’)  were  plotted  on  a 
logarithmic  scale  in  order  to  determine  the  value  of  p 
in  the  equation  7c=const  (l-777j)p  (Fig.  3).  As  can  be 
seen,  for  D  samples  7  (7)  corresponds  to  an  SIS  junc¬ 
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Fig.  3.  Critical  current  versus  (1-T /T  )  plotted  for  D  and  S 
samples  in  a  logarithmic  scale  to  determine  the  value  of  b  in 
the  equation  7  =const  (1  -T/T  )F  P=1  and  |3=2  cases  are 
shown  by  the  dashed  lines. 


tion  (P~l).  On  the  other  hand,  for  S  samples  7(7)  agrees 
with  an  SNS  type  junction  (p«2)  for  any  x  value. 

In  this  work  a  systematic  statistical  study  of  grain 
boundaries  was  made  to  correlate  their  distribution  to 
the  observed  change  in  J  when  x  varied  from  0  to  0.4. 
So,  40-50  consecutive  boundaries  in  each  sample  were 
characterized  by  bright  field  tilting  imaging  in  TEM. 
The  grain  boundaries  were  distributed  into  3  types. 
‘Clean’  boundaries  correspond  to  boundaries  in  which 
no  intergranular  phase  was  detected  in  TEM  bright 
field.  In  contrast,  ‘dirty’  boundaries  exhibit  an  inter¬ 
granular  phase  (thickness  >2  nm).  At  last,  ‘thin-film- 
coated’  boundaries  display  a  thin  phase  with  a  thick¬ 
ness  <2  nm.  For  clean  boundaries,  we  paid  attention 
on  peculiarities  such  as  dislocations,  small  facets  or 
fine  precipitates  on  them.  One  should  take  into  account 
the  fact  that  clean  boundaries  are  favourable  for  trans¬ 
mitting  the  supercurrent  only  when  the  orientation 
crystallographic  relationships  between  adjacent  grains 
are  suitable.  Usually,  no  more  than  50%  of  clean  bound¬ 
aries  are  considered  as  favourable.  According  to  the 
data  in  [9],  the  percolation  threshold  is  about  12%. 

The  histogram  in  Fig.  4  shows  the  statistical  study 
results.  In  the  D  samples,  the  ratio  of  clean  boundaries 
was  found  to  remain  between  27  and  45%.  The  pro¬ 
portion  of  clean  boundaries  decreased  with  x.  The  ob¬ 
served  increase  in  J  with  x  in  D  samples  can  be  ex¬ 
plained  by  the  increase  in  the  proportion  of  the  faceted 
clean  boundaries,  which  can  act  as  pinning  centres  (Fig. 
4).  Ag  doping  resulted  in  a  higher  proportion  of  clean 
boundaries  forx>0.2,  but  there  were  only  few  faceted 
boundaries  among  them.  For  the  composition  with 
x=0.4,  which  exhibited  3  times  increase  in  J ,  a  pecu¬ 
liarity  in  the  structure  of  clean  boundaries  was  found, 
namely,  about  half  of  the  clean  boundaries  contained 


Fig.  4.  Histogram  showing  the  distribution  of  clean,  thin- 
film -coated  and  dirty  grain  boundaries  for  a  sample  of  the 
stoichiometric  composition  (x=0),  D  and  S  samples.  The 
dashed  lines  show  the  fraction  of  faceted  clean  boundaries 
among  all  clean  boundaries. 
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Fig.  5.  Micrograph  of  a  ‘clean’  boundary  with  Ag  precipitates 
(shown  by  arrows)  in  TEM  for  an  S  sample  (.\=0.4). 


small  precipitates  of  2-5  nrn  in  size  on  the  boundary 
itself,  or  close  to  it  (up  to  5  nm).  The  average  distance 
between  these  precipitates  was  about  30  nm.  An  ex¬ 
ample  of  such  a  ‘clean’  boundary  with  Ag  precipitates 
is  presented  in  Fig.  5.  These  small  Ag  precipitates  can 
act  as  very  efficient  pinning  centres  resulting  most  likely 
in  the  3-fold  increase  in  .7  in  this  composition. 

To  detect  possible  Ag  segregation  at  clean  grain 
boundaries,  EDX  analysis  in  STEM  with  1  nm  probe 
was  used.  The  detected  concentration  of  Ag  on  some 
clean  boundaries  was  twice  higher  than  the  obtained 
solubility  limit  of  Ag  in  the  grains  (x=0.03).  This  in¬ 
creased  Ag  segregation  on  grain  boundaries  most  prob¬ 
ably  leads  to  SNS  behaviour  of  weak  links  in  S  samples 
as  result  of  a  percolation  path  of  supercurrent  through 
such  ‘clean’  boundaries  with  an  extremely  narrow  Ag 
segregation  (si  nm)  on  them. 


3.2  Cu-DEFICIENT  DYSPROSIUM  CERAMICS 

DOPED  WITH  1  wt.%  Pt 

The  next  system  studied  was  Cu-deficient  dysprosium 
ceramics  doped  by  Pt.  Such  a  ceramic  system  was 
chosen  due  to  the  following  reason.  As  we  saw  above, 
the  deviation  from  stoichiometry  leads  to  the  forma¬ 
tion  of  secondary  phases  in  sintered  ceramics.  Accord¬ 
ing  to  the  literature  data  [10.  11],  in  melt-textured  ce¬ 
ramics  the  addition  of  Pt  can  result  in  decrease  in  size 
of  2 1 1  secondary  phase  inclusions,  improving  thereby 
their  pinning  possibility  and.  then  the  critical  current. 
In  our  work,  for  Cu-deficient  dysprosium  sitered  ce¬ 
ramics  doped  with  Pt.  we  tried  again  to  identify  micro- 
structural  elements,  which  can  be  responsible  for  the 
differences  in  superconducting  properties. 

In  all  compositions,  Pt-doping  resulted  in  the  nucle- 
ation  of  intragranular  secondary  phases  of  submicron 


Fig.  6.  Critical  current  density  at  77  K  versus  applied  magnetic 
field  for  DyBa2Cu}  xOy/lwt.%Pt  and  DyBa2Cu3  samples. 


size,  which  contain  Pt  and  were  homogeneously  dis¬ 
tributed  into  123  grains.  In  addition,  in  the  non-sto- 
ichiometric  composition,  Pt  was  also  located  in  in¬ 
tergranular  BaCu02  secondary  phases.  Unlike  the  melt- 
textured  123  ceramics  [10.  11],  no  influence  of  Pt  on 
size  and  distribution  of  2 1 1  secondary  phase  was  found. 
Traces  of  Pt  in  the  123  grains  were  detected  in  all 
compositions  doped  by  Pt.  Distribution  of  Pt  in  the 
matrix  was  inhomogeneous.  The  maximum  content  of 
Pt  in  a  grain  was  0.06  molar  content. 

Fig.  6  displays  the  dependence  of  critical  current 
density  J  on  applied  magnetic  field.  As  seen,  at  zero- 
magnetic  field  doping  with  platinum  practically  does 
not  change  J  in  stoichiometric  composition,  but  leads 
to  its  substantial  increase  for  non-stoichiometric 
composition  with  x=0.2.  The  stoichiometric  compo¬ 
sitions  exhibit  a  sharp  drop  of  J,  values  even  in  very 
low  magnetic  field.  Non-stoichiometric  compositions 
show  a  much  smoother  decrease  in  J  with  magnetic 
field.  Doping  with  Pt  strongly  contributed  to  the  so- 
called  ‘fish-tail’  effect  which  means  stabilization  of  J 
in  magnetic  field.  For  this  composition  this  ‘fish-tail’ 
behaviour  became  more  pronounced  when  the  tempera¬ 
ture  was  decreased  (Fig.  7).  So,  the  copper-deficient 
samples,  especially  the  one  doped  with  Pt  seem  to  have 
better  possibilities  for  pinning  magnetic  vortices  in 
applied  magnetic  fields. 

Similar  to  the  system  with  Ag,  for  these  ceramics 
with  Pt  we  also  tried  to  identify  the  character  of  the 
intergranular  junction  network  from  /  versus  tempera¬ 
ture  dependence  (Fig.  8).  Grain  boundary  network  was 
found  to  behave  predominantly  as  SIS  junctions,  even 
in  Pt-doped  samples  so  that  Pt  seems  not  to  segregate 
in  sufficient  amount,  at  least,  at  grain  boundaries. 
Within  the  accuracy  of  our  EDX  analysis  we  did  not 
reveal  higher  concentration  of  Pt  in  boundaries 
compared  with  grains. 
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Fig.  7.  Critical  current  density  at  77  K  versus  applied  magnetic 
field  for  DyBa2Cu2gOy/lwt.%Pt  at  different  temperatures. 

The  above  superconducting  properties  were 
correlated  with  the  microstructural  features  of  the  stud¬ 
ied  samples.  In  all  compositions  the  123  matrix  exhib¬ 
its  well  crystallized  and  well  twinned  grains,  but  with 
different  twin  density.  It  is  known  that  twin  bound¬ 
aries  in  the  123  materials  may  act  as  flux  pinning  cen¬ 
tres  [12,  13].  In  granular  ceramics,  where  critical 
current  is  limited  by  weak  links  at  grain  boundaries,  it 
is  reasonable  to  consider  twins  near  the  boundaries. 
The  higher  the  twin  density  the  more  efficient  the  pin¬ 
ning  of  vortices  along  the  boundary.  We  compared  the 
twin  density  in  our  compositions.  As  seen  in  the  Table 
2,  in  non-stoichiometric  composition,  the  average  spac¬ 
ing  d  between  twins  is  twice  smaller  than  in  the  sto¬ 
ichiometric  one.  In  addition,  in  non-stoichiometric 
composition  doped  with  Pt  we  locally  observed  very 
fine  twins  with  d=20  nm,  whereas  in  such  composition 
only  without  Pt  such  a  high  density  of  twins  was  not 
found,  at  all.  The  presence  of  these  highly  dense  twins 
can  be  caused  by  partial  substitution  of  Cu  by  Pt  in 
grains.  A  decrease  of  twin  spacing  due  to  Pt  doping 
was  observed  in  melt-textured  123  ceramics  in  [13, 14], 
The  higher  density  of  twin  boundaries  in  the  vicinity 
of  grain  boundaries  can  improve  the  pinning  proper¬ 
ties  under  magnetic  field.  The  highest  local  density  of 
twins  was  observed  in  the  composition,  which  exhib¬ 
ited  the  best  ‘fish-tail’  effect  in  magnetic  field. 


Fig.  8.  Critical  current  versus  1-T/T  plotted  on  logarithmic 
scale  to  determine  the  value  [3  in  the  equation  /  =const  (1-T / 
T)p  for  DyBa2Cu3  xO  /lwt.%Pt  and  DyBa2Cu3  Oy  samples. 


We  also  performed  comparative  statistical  study  of 
grain  boundaries,  Table  3.  It  was  found  that  in  non- 
stoichiometric  composition,  unlike  the  stoichiometric 
one,  grain  boundaries  are  mostly  clean  or  thin-film- 
coated.  Doping  with  Pt  resulted  in  a  higher  proportion 
of  clean  boundaries.  Such  statistical  results  can  explain 
the  sharp  drop  of  Ic  under  very  low  magnetic  fields 
(<10  G)  in  DyBa2Cu3Oy/lwt.%Pt  samples.  As  already 
noted  above,  no  more  than  a  half  of  clean  boundaries 
can  be  considered  as  favourable  [9],  So  the  supercur¬ 
rent  is  likely  to  pass  through  a  few  thin-film-coated 
GBs  resulting  in  the  sharp  drop  in  /  in  low  magnetic 
fields.  The  non-stoichiometric  compositions  with  x=0.2 
have  enough  clean  boundaries  to  provide  the  percola¬ 
tion  path  for  supercurrent.  A  possibility  for 
redistribution  of  the  supercurrent  through  clean  bound¬ 
aries  in  a  magnetic  field  and  high  dense  twins  locally 
observed  can  be  reasons  for  better  I  behaviour  in  a 

c 

magnetic  field  in  non-stoichiometric  composition  doped 
with  Pt. 

CONCLUSION 

(i)  In  all  considered  ceramic  systems  the  critical  current 
/  is  controlled  by  weak  links  at  grain  boundaries. 


Table  2.  The  spacing  between  twins  in  DyBaCuO(Pt)  samples. 


DyBcpCu  3Oy/lwt.  %Pl 

DyBa7Cu,sOy/lwt.  %Pt 

DyBafu:/)y 

Average  twin  spacing  d 

110  nm 

50  nm  (fine  twins  with 
d*20  nm  locally  observed) 

60  nm 
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Table  3.  Tlic  distribution  of  clean,  thin-film-coated  and  dirty  grain  boundaries  in  DyBa2Cu3 .,Oy/lwt%  Pt  samples. 


DyBafnp/lwt%  Pt 

DyBa,Cu2  £) J Iwt%  Pt 

DyBa,Cu,sOy 

Clean 

24  % 

58  % 

44% 

Thin-film 

16  % 

29% 

37% 

Dirty 

56  % 

5  % 

15  % 

Non-characterized 

4% 

8  % 

4  % 

(ii)  In  Cu-deficient  yttrium  ceramics  doped  with  Ag  in 
amount  equal  the  deficiency,  Ag  doping  leads  to 
change  of  type  of  the  weak  links  from  SIS  to  SNS- 
type.  Such  SNS  behaviour  is  a  result  of  percolation 
path  of  supercurrent  through  clean  boundaries  with 
extremely  narrow  (lnm)  Ag-segregation  on  them. 
For  x=0.4,  decoration  of  clean  boundaries  by  very 
fine  Ag-precipitates  (2-5  mn)  is  the  most  likely 
reason  for  3  fold  increase  of  critical  current, 

(iii)  In  non-stoichiometric  composition  with  x=0.2,  Pt- 
doping  leads  to  pronounced  increase  in  Ic  and  fish¬ 
tail  effect  in  magnetic  field.  This  may  be  associ¬ 
ated  with  the  substantial  increase  in  proportion  of 
clean  boundaries  and  highly  dense  twins  locally 
observed  (due  to  Pt  segregation  in  123  grains).  No 
noticeable  segregation  of  Pt  in  grain  boundaries  was 
revealed. 
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Abstract.  The  size  effect  in  thin  film  sandwich  structures  is  considered.  Three  types  of  boundary  conditions 
for  dynamic  polarization  at  interface  electrode-ferroelectric  layer  are  formulated.  The  type  of  boundary 
conditions  depends  on  matching  crystal  lattices  of  electrodes  and  ferroelectric  layer.  A  model  describing  the 
dependence  of  dielectric  permittivity  on  biasing  field,  temperature,  and  thickness  of  the  thin  ferroelectric 
film  is  proposed.  The  results  of  the  simulation  are  in  good  agreement  with  experiments. 


The  dielectric  nonlinearity  of  ferroelectrics  allows  to 
use  these  materials  as  the  basis  of  electrically  tunable 
devices.  Thin  film  ferroelectrics  are  used  at  microwaves. 
If  the  thickness  of  the  film  is  comparable  with  the 
correlation  radius  the  size  effect  appears  in  the  film, 
i.e.  the  dielectric  permittivity  of  film  depends  on  its 
thickness.  The  size  effect  is  related  with  the  spatial 
distribution  of  polarization  inside  the  ferroelectric  film 
and  type  of  the  boundary  conditions  for  polarization. 
The  size  effect  results  in  a  decrease  in  the  polarization 
and  therefore,  in  reduced  dielectric  permittivity  and 
tunability.  Through  the  selection  of  electrode  material 
one  can  control  the  size  effect  on  dielectric  response  of 
the  thin  ferroelectric  films. 

The  dielectric  permittivity  of  ferroelectric  bulk  ma¬ 
terial  is  a  function  of  biasing  field  and  temperature.  The 
model  of  dependence  of  dielectric  permittivity  on  bias¬ 
ing  field  and  temperature  is  based  on  Ginsburg- 
Devonshire  expansion  of  free  energy  density  in  the  power 
series  over  the  order  parameter.  For  ferroelectrics  the 
order  parameter  is  a  spontaneous  polarization  [1-4], 

The  phenomenological  model  of  the  dielectric  per¬ 
mittivity  of  a  bulk  sample  is  presented  by  the  follow¬ 
ing  equation: 

e~'  + 

[(£2  +V)U1  -9k]  -7}  (1) 


where 

=  >/£  +  £’  (2) 


n(T)  = 


0, 


Tc  \ 


1  1 

f  T\ 

—  + 

16  1 

-1, 


(3) 


Eh ' 


£„  =• 


2  D„ 


*0(3  O’ 


(4) 


ew  is  the  permittivity  of  a  free  space,  h  is  the  thickness 
of  ferroelectric  film  included  in  sandwich  capacitor 
represented  in  Fig.  1,  UB  is  a  biasing  voltage. 

Model  parameters:  sQ0  is  the  analogue  to  Curie- 
Weiss  constant;  T(.  is  the  effective  Curie  temperature; 
9  is  the  effective  Debye  temperature  of  the  sublattice 
oscillations  causing  ferroelectric  polarization;  EN  is  the 


Fig.  1.  Sandwich  capacitor  as  a  simplified  general  structure. 
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Table  i.  Mode!  parameters  for  single  crystal  SrTiOy 


7'r(K) 

Eoo 

£v(kV/cm) 

6,(K) 

$.v 

single  crystal 

42 

2081 

19.3 

175 

0.018 

SrTi03  [5] 

normalizing  electric  field;  lJX  is  the  statistical  disper¬ 
sion  of  the  biasing  field  characterizing  the  quality  of 
material.  The  numerical  values  of  model  parameters 
for  bulk  material  [5]  are  presented  in  Table  1. 

The  dependence  of  dielectric  permittivity  of  the  film 
on  its  thickness  should  be  included  in  the  model.  In 
order  to  take  into  account  the  size  effect,  it  is  necessary 
to  solve  a  second-order  differential  equation  with 
respect  to  polarization.  This  equation  is  a  consequence 
of  Ginsburg-Devonshire  expansion  [6]: 


(I)  Zero  boundary  conditions: 


(II)  Intermediate  boundary  conditions: 


f 

Pic(x)  ±  b 
V 


dPJx)'' 
dx  j 


=  0 


h  <  h  <  a) 


(III)  Free  boundary  conditions: 


dlP(x)  D(x)  D\x) 

-2X,  :  t  +  :  -£„£(*),  (5) 


dx 


e{T)  D\ 


dP  (x) 

ac  v  /  | 

dx 
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h  — >  co 

b-ir  ^ 
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e(T)  = 

V  (T) 


(6) 


where  P(x),  D(x),  E(x)  are  polarization,  displacement, 
and  electric  field;  X,  is  a  correlation  parameter;  x  axis 
is  directed  normally  to  the  electrode-ferroelectric  film 
interface. 

The  numerical  values  of  the  correlation  parameter 
determined  from  the  experimental  data  on  inelastic 
neutron  scattering  on  soft  mode  of  ferroelectric  crystal 
are  presented  in  Table  2  [7,8],  The  film  thickness  above 
which  the  size  effect  appears  is  a  correlation  radius  of 
ferroelectric  polarization  which  can  be  extracted  from 
the  equation  (5). 


Table  2.  Numerical  values  of  correlation  parameter. 


SrTi03 

KTa03 

BaTi03* 

pm2 

1.15  1 0'7 

5.371 0'7 

io-9 

*  Preliminary  estimation 


r(T)  =  ^2Al£(T). 


(7) 


The  boundary  conditions  of  the  type  (I)  are  refered 
to  a  sandwich  capacitor  with  electrodes  made  of  a 
normal  metal.  The  polarization  distribution  provided 
by  these  boundary  conditions  causes  the  size  effect.  The 
boundary  conditions  of  the  type  (II)  are  refered  in  a 
sandwich  capacitor  with  YBa,Cu307>.  (YBCO)  elec¬ 
trodes.  In  this  case  polarization  can  partially  penetrate 
inside  the  electrodes.  To  allow  the  better  understand¬ 
ing  of  that,  the  effective  film  thickness  h^f  was  intro¬ 
duced.  Free  boundary  conditions  (type  III)  are  realized 
in  structures  with  SrRuO,  (SRO)  electrodes.  The  spa¬ 
tial  distribution  of  polarization  inside  the  ferroelectric 
layer  for  three  types  of  boundary  conditions  is  shown 
in  Fig.  2. 

Having  solved  equation  (5)  with  zero  boundary 
conditions,  one  can  derive  the  expression  for  dynamic 
polarization  taking  into  account  the  size  effect: 


PJx)  =  Pm  (0)  1 


ch(ax) 
ch(a  h  /  2) 


s(T)  V  DyS  ) 


For  different  crystal  structures  of  electrodes  three 
different  types  of  boundary  conditions  at  electrode-fer¬ 
roelectric  interface  can  be  realized.  The  spatial  distri¬ 
bution  of  polarization  corresponding  to  a  certain  type 
of  boundary  conditions  induces  the  size  effect  in  the 
film. 

There  are  three  types  of  boundary  conditions  for 
dynamic  polarization  [9]: 


where  Qa  and  Qk  are  the  alternating  and  direct  current 
components  of  the  charge  at  the  electrodes  respectively; 
S  is  the  area  of  the  electrodes. 

The  inverse  effective  dielectric  permittivity  of  a  thin 
ferroelectric  film  is: 
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Fig.  2.  Spatial  distribution  of  dynamic  polarization  inside  ferroelectric  layer  in  the  case:  a  —  zero  boundary  conditions;  b  — 
intermediate  boundary  conditions;  c  —  free  boundary  conditions. 


“3 
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(-2  ,  KW  , 

[(£  +  n)  + 


+[(?  +  n)'l2  +  z]2 


-77  + a2}  (8) 

where  the  parameter  of  size  effect  is  a  =  ^2em  /  ahejr ; 
a-\! 

The  numerical  values  of  the  model  parameters  ob¬ 
tained  for  various  sandwich  structures  are  presented 
in  Table  3.  Fig.  3  shows  good  agreement  between  ex¬ 
perimental  and  model  dependencies  calculated  using 


model  parameters  from  Table  3.  Curves  in  Fig.  3  illus¬ 
trate  influence  of  size  effect  on  dielectric  permittivity. 
In  the  ferroelectric  capacitors  with  YBCO  electrodes 
dynamic  polarization  partially  penetrates  inside  elec¬ 
trodes  and  influence  of  the  size  effect  in  this  case  is 
attenuated. 

Experimental  [13]  and  model  dependencies  of  di¬ 
electric  permittivity  of  thin  film  Ba0  |2Sr0  8STi03  (BSTO) 
included  in  sandwich  capacitors  with  Pt  and  SRO  elec¬ 
trodes  are  presented  in  Fig.  4.  At  room  temperature 


Tabic  3.  Model  parameters  for  thin  ferroelectric  films. 


Capacitor 

hfom) 

hj^m) 

sna 

TC(K) 

Q,(K) 

EN(kV/cm) 

5, 

a 

YBCO/STO/Au  [10] 

0.250 

0.230 

4265 

34 

152 

6.0 

2.3 

16 

(I) 

Pt/STO/Ni  [11] 

0.046 

0.037 

4143 

40 

152 

8.7 

1.5 

7 

(I) 

YBCO/STO/YBCO  [12] 

0.800 

2.300 

3090 

35 

152 

6.5 

0.8 

1.2 

(II) 

SRO/BSTO/SRO  [13] 

0.020 

OO 

4400 

42 

175 

12 

1.0 

0 

(III) 

SRO/BSTO/Pt  [13] 

0.020 

0.37 

4400 

42 

175 

12 

1.0 

3 

(II) 
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effective  permittivity  of  the  film  in  capacitor  with  SRO 
electrodes  is  higher  than  permittivity  of  single  crystal 
in  about  2.5  times.  One  can  conclude  that  in  the  case 
of  good  matching  of  crystal  lattices  of  electrodes  and 
ferroelectric  film  the  size  effect  is  suppressed  and  free 
boundary  conditions  for  dynamic  polarization  are 
realized.  Tunability  of  such  capacitor  is  more  better 
than  tunability  of  the  capacitor  with  Pt  electrode. 

Thus  interface  between  thin  ferroelectric  film  and 
electrodes  determines  distribution  of  polarization  inside 
ferroelectric  layer  and  therefore  variation  of  dielectric 
characteristics  of  sandwich  capacitor.  By  chosen  elec¬ 
trodes  one  can  suppress  size  effect  in  thin  film  struc¬ 
tures. 
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Abstract.  A  theoretical  model  is  suggested  which  describes  the  effect  of  tilt  boundaries  with  chaotically 
arranged  dislocations  on  the  critical  current  density  in  high-transition-temperature  (7))  superconductors. 
Stress  fields  of  such  boundaries,  that  suppress  high-7'  superconductivity,  are  revealed  to  be  long-range  as 
compared  to  those  of  tilt  boundaries  with  periodically  ordered  grain  boundary  dislocations.  With  this  factor 
taken  into  consideration,  the  dependence  of  the  critical  current  density  accross  tilt  boundaries  on  boundary 
misorientation  0  is  calculated  and  compared  with  the  experimental  data  [1]. 


1.  INTRODUCTION 

The  transport  properties  of  grain  boundaries  in  high  - 
Tc  superconducting  materials  are  the  subject  of  inten¬ 
sive  theoretical  and  experimental  studies  (e.g.,  see  [1- 
14])  motivated  by  extremely  large  interest  to 
applications  of  such  materials.  It  has  been  experimen¬ 
tally  revealed  that  grain  boundaries  drastically  suppress 
the  superconducting  critical  current  [1-3],  In  particu¬ 
lar,  the  critical  current  density,  /.,  a.cross  tilt  bound¬ 
aries  in  YBa2Cu,07S  superconductors  decreases  with 
boundary  misorientation  angle  9  as  follows  :  jc( 9)*  jc 
(0°)exp(-9/8°),  for  9  ranging  from  0°  to  1 5°.  The  critical 
current  density  yc(9)  across  high-angle  tilt  boundaries 
with  0  >  15°  is  tentatively  constant  and  low:  yc(0  > 
15°)  *  0.02  y)  (0°)  [1-3]. 

Microscopic  mechanism  of  the  grain  boundary  ef¬ 
fect  on  high-7’  superconductivity  is  not  understood  yet, 
though  several  models  of  the  effect  in  question  have 
been  suggested  [5-14].  In  particular,  theoretical  mod¬ 
els  have  been  elaborated  describing  the  grain  bound¬ 
ary  effect  as  that  related  to  the  following  factors  :  (i)  a 
decrease  of  electron  free  path  in  vicinities  of  grain 
boundaries  [6];  (/'/')  stress-field-induced  structural  dis¬ 
order  and  pinning  of  magnetic  vortices  at  grain  bound¬ 
aries  [7-9];  (Hi)  formation  of  the  antiferromagnetic 
phase  within  grain  boundary  cores  [10];  (/v)  faceting 
of  grain  boundaries  and  d  -symmetry  of  the  supercon¬ 
ducting  order  parameter  [11];  (v)  deviations  from  bulk 
stoichiometry  in  vicinities  of  grain  boundaries  [12-14]. 


Most  models  discussed  treat  stress  fields  of  grain  bound¬ 
ary  dislocations  as  those  responsible  for  suppression 
of  the  superconducting  order  parameter  in  vicinities  of 
grain  boundaries  in  high-7),  superconductors  [7-10]. 
Though  the  concrete  mechanism  of  the  suppression  is 
not  unambigiously  recognized  yet,  many  researchers 
assume  the  existence  of  a  critical  stress,  oc,  which 
characterizes  superconductor-to-non-superconductor 
transition  in  stressed  regions  [7-10].  More  precisely, 
superconductivity  is  supposed  to  be  completely  sup¬ 
pressed  in  regions  where  some  component(s)  of  stress 
tensor  is  (are)  larger  than  ac.  However,  this  idea  is  in 
contradiction  with  experimental  data  in  case  of  low- 
angle  boundaries  consisting  of  periodically  arranged 
lattice  dislocations  (Fig. la).  Actually,  for  any  value  of 
c  ,  the  volume  fraction  of  “critically  stressed”  regions 
(regions  where  some  component(s)  of  stress  tensor  is 
(are)  larger  than  o)  in  vicinities  of  such  low-angle 
boundaries  (that  create  short-  or  intermediate-range 
stress  fields  which  drop  as  exp(-x/hB),  where  x  denotes 
the  distance  from  a  boundary  plane,  and  hQ  the 
interspacing  (period)  between  periodically  arranged 
dislocations,  see  (Figs,  la  and  2a)  decreases  with 
boundary  misorientation  0. 

This,  in  context  of  the  idea  on  existence  of  the 
critical  stress  sc,  results  in  increase  of  the  critical  current 
density,  j ,  with  boundary  misorientation  0,  while  ex¬ 
periments  [1-3]  are  indicative  of  the  fact  that  jc  drasti¬ 
cally  decreases  with  9  in  the  range  of  9  from  0°  to 
tentatively  15°.  In  order  to  avoid  the  contradiction  dis- 
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(a)  :  (b)  • 

Fig.  1.  Ordered  and  disordered,  low-angle  tilt  boundaries  are 
represented  as  (a)  periodically  ordered  and  (b)  disordered  walls 
of  crystal  lattice  dislocations  of  the  edge  type,  respectively. 


(b) 


Fig.  2.  High  -angle  tilt  boundaries  with  (a)  periodically  ordered 
and  (b)  disordered  ensembles  of  grain  boundary  dislocations. 


cussed,  Gurevich  and  Pashitskii  [10]  assume  that  the 
non-superconducting,  critically  stressed  regions  consist 
of  two,  normal  metal  and  antiferromagnetic,  phases, 
in  which  case  the  antiferromagnetic  phase  effectively 
destroys  the  coupling  of  Cooper  electrons  in  vicinities 
of  grain  boundaries. 

In  this  paper  we  suggest  an  alternative  model  of 
the  grain  boundary  effect  on  high-Tc  superconductiv¬ 
ity,  based  on  both  the  idea  on  critical  stress  a.  and  that 
on  disorder  in  spatial  arrangement  of  grain  boundary 
dislocations  in  high- 7  superconductors.  Actually,  high- 
T  superconducting  materials  commonly  are  synthe¬ 
sized  at  non-equilibrium  conditions  at  which  grain 
boundary  dislocations  can  form  “non-equilibrium” 
structures  (Figs,  lb  and  2b)  that  are  different  from 
“equilibrium”,  periodic  ones  (Figs,  la  and  2a).  It  this 
situation,  grain  boundaries  in  high- 7],  superconduct¬ 
ing  materials,  as  with  boundaries  in  conventional  ma¬ 
terials  synthesized  at  non-equilibrium  conditions 
[15,16],  create  long-range  stress  fields  which  drop  as 
rw  or  x'm.  The  main  aim  of  this  paper  is  to  calculate 
the  angular  dependence  of  the  critical  current  density, 
j  (9),  in  high-7  superconductors  with  tilt  boundaries 
containing  disorderedly  distributed  dislocations. 

2.  MODELS  OF  LOW-  AND  HIGH- 
ANGLE  TILT  BOUNDARIES  IN 
HIGH-Tc  SUPERCONDUCTORS 

Stress  fields  of  grain  boundaries,  in  fact,  are  created  by 
boundary  defects  (first  of  all,  dislocations)  whose  pa¬ 
rameters  are  caused  by  the  structural  geometry  of 
boundaries.  So,  low-angle  tilt  boundaries  (with 
misorientation  G  ranging  from  0°  to  0c  *  15°)  are  ef¬ 


fectively  represented  as  walls  of  lattice  dislocations  (of 
the  edge  type)  with  Burgers  vectors  being  those  of  the 
crystalline  lattice  [17,18]  (Fig.  1 ).  For  0  >  0c,  cores  of 
lattice  dislocations  overlap  which  compose  grain 
boundaries.  As  a  corollary,  high-angle  tilt  boundaries 
(with  9  >  0 )  have  layer-like  cores  (characterized  by 
thickness  «  1  nm)  (Fig.2)  with  the  structure  being  dif¬ 
ferent  from  the  crystalline  structure  of  adjacent  grains. 
The  layer-like  cores  of  high-angle  boundaries 
commonly  contain  grain  boundary  dislocations  (Fig.2) 
with  Burgers  vectors  being  those  of  displacement-shift- 
complete  lattices  of  high-angle  boundaries  [18].  Moduli 
of  such  Burgers  vectors  commonly  are  lower  than 
moduli  of  Burgers  vectors  of  crystal  lattice  dislocations 
composing  low-angle  boundaries. 

Grain  boundaries  at  equilibrium  conditions  form 
either  periodic  (Figs.  1  a  and  2a)  [  1 7, 1 8]  or  quasiperiodic 
[18-20]  ensembles  that  correspond  to  minimum  elas¬ 
tic  energy  density  of  such  boundaries.  Grain  bound¬ 
aries  in  poiycrystalline  materials  synthesized  at  highly 
non-equilibrium  conditions  are  non-equilibrium  in  the 
sense  that  they  are  characterized  by  disorder  in  spatial 
arrangement  of  boundary  dislocations  (see  Figs,  lb  and 
2b).  In  doing  so,  ensembles  of  grain  boundary  disloca¬ 
tions  at  high-angle  boundaries  (Fig.  2b),  in  general, 
are  more  disorderedly  distributed  than  lattice  disloca¬ 
tions  composing  low-angle  boundaries  (Fig.  lb).  The 
aforesaid  is  related  to  the  fact  that  dislocations  at  high- 
angle  boundaries  are  characterized  by  small  Burgers 
vectors  as  compared  to  those  at  low-angle  boundaries; 
in  this  case  displacements  of  dislocations  from  their 
equilibrium  positions  at  high-angle  boundaries  are  ac¬ 
companied  by  low  changes  in  the  elastic  energy  den- 
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sity  as  compared  to  those  of  dislocations  composing 
low-angle  boundaries. 

Now  let  us  consider  a  low-angle  tilt  boundary  of 
infinite  extent,  which  consists  of  lattice  dislocations 
periodically  spaced  in  plane  yz  in  coordinate  system 
shown  in  (Fig.  la).  Such  dislocations  are  characterized 
by  ordinates  yn  =  rth0  and  Burgers  vectors  b  =  (b, 0,0) 
with  h  and  b  being  in  the  Frank’s  relationship  [18]: 
b=  2h  sin(0/2).  Dislocations  at  “non-equilibrium”  tilt 
boundary  (Fig.  lb)  are  assumed  to  be  displaced  by 
hQ  5 n(-fl  <  8n  <  a)  in  a  randomly  uniform  way  from 
their  “equilibrium”  positions  corresponding  to  the  pe¬ 
riodic  arrangement  shown  in  Fig.  la.  In  other  words, 
the  probability  of  finding  the  /7th  dislocation  in  an 
infinitesimal  interval  (y,y  +  dy)  belonging  to  the  inter¬ 
val  y-a  <y  <yn  +  a  is  equal  to  dy/2a.  In  the  situation 
discussed,  ordinates  of  the  dislocations  composing  the 
“non-equilibrium”  tilt  boundary  (Fig.  lb)  are  given  as: 

Kn  +  SJ 

y"  2sin(0/2)  ’  (1) 

where 


<8  >=  0,  <  52  >=  —  ,<85  >=  0 , 

i,  ’  ti  j  " 


(2) 


with  <...>  being  the  averaging  on  the  non-equilibrium 
dislocation  ensemble.  The  dispersion  of  grain  bound¬ 
ary  stress  is  equal  to  the  sum  of  the  dispersions  of 
stresses  induced  by  the  boundary  dislocations: 


Daafl(x’y)  =  Z  >  - 

11  = -re 

<Wall{x,y)]>}  ( a,/3={x,y )) 


(3) 


with  cr (“2  (on  the  r.h.s.  of  formula  (3))  being  the  stress 
field  of  the  /7th  dislocation,  which,  following  [17],  can 
be  written  as: 


Here  G  denotes  the  shear  modulus,  and  p  the  Poisson 
ratio. 

After  some  algebra  based  on  the  calculation  scheme 
[15],  from  (l)-(7)  we  have  found  the  averaged  (on  the 
ordinate  y)  dispersion,  Doap(x),  of  stress  tensor 
components  in  vicinity  of  a  non-equilibrium  tilt  bound¬ 
ary.  Following  the  results  of  our  calculations,  the  largest 
dispersion  is  that  of  cr  : 
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Der  v(x)  = 
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[  2rrhn(\- /u)J 
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+  a 
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For  any  x,  values  of  Do  and  Doyy  are  lower  by  factor 
5,  while  values  of  Do__  (for  p=0.2;  see  experimental 
data  [21])  are  lower  by  factor  ~  4  than  those  of  Doxx. 
Therefore,  in  context  of  the  idea  [7-10]  on  the  crucial 
effect  of  the  largest  stress  tensor  component  on  high-7] 
superconductivity,  here  and  in  the  following  we  will 
focus  on  analysis  of  the  dispersion  Doxx.  In  doing  so, 
we  find  the  mean  value  of  modulus  of  oxx  in  vicinity  of 
a  non-equilibrium  tilt  boundary  of  infinite  extent  to  be 
as  follows  : 


(10) 


because  <ox>  =  0. 

The  dispersion  of  stress  tensor  components  in  the 
case  of  grain  boundaries  with  partly  relaxed  disloca¬ 
tion  distributions  can  be  calculated  by  the  same  methods 
as  with  previously  examined  non-relaxed  disordered 
dislocation  ensembles.  So,  in  paper  [16]  it  has  been 
shown  that  the  dispersion  in  question  can  be  written  as 
A  2  Da  with  A  (0  <  A  <1)  being  the  factor  of  relaxation. 
In  these  circumstances,  the  mean  value  of  modulus  of 
c  depends  on  the  distance,  x,  from  a  partly  relaxed 
tilt  boundary  as  follows: 
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According  to  formula  (10),  o{x)  is  strongly  affected  by 
the  two  non-dimensional  parameters,  a  and  A:  a  oc  x'3/2, 
if  x  <1  ab/2  sin(9/2);  and  a  «  x‘l/2,  if  x  <1  ab/2  sin(0/2) 

[15]. 


(7) 
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Fig. 3.  Theoretical  (solid  curve)  and  experimental  [1]  (dots) 
dependences  of  the  averaged  density  (j  )  of  the  superconducting 
critical  current  accross  “non-equilibrium”  tilt  boundaries  on 
boundary  misorientation  (9). 

3.  CRITICAL  CURRENT  DENSITY 
ACROSS  DISORDERED  TILT 
BOUNDARIES 

In  spirit  of  models  [7-10],  let  us  assume  that  high-7^ 
superconductivity  disappears  in  stressed  regions  in 
which  o(x)  >  ctc  =  a  G,  where  o.  <s  1.  Following  [8], 
a  ~  1  O'2.  In  the  situation  discussed,  we  find  that  the 

c 

non-superconducting  phase  in  vicinity  of  a  tilt  bound¬ 
ary  with  a  partly  relaxed  dislocation  distribution  lies 
in  the  layer:  -x  <  x  <  xc,  where  xc  is  extracted  from 
equation  c(x  =  xc)  =  a  G.  For  definiteness,  let  us  as¬ 
sume  that  the  non-superconducting  phase  exhibits  the 
normal  metal  properties.  Then,  according  to  the  theory 
of  Cooper  pairs  tunneling  [22,23],  we  find  the  critical 
current  density  across  a  normal  metal  layer  of  thick¬ 
ness  2x  in  vicinity  of  the  tilt  boundary  to  be  as  fol¬ 
lows: 

j  (#) 

— - =  exp(-2x  (0)  /  £),  (11) 

.7,(0°) 

where  E,  is  the  electron  coherence  length  in  the  normal 
metal  phase.  Formula  (10)  and  (11)  allow  one  to 
calculate  xc(0)  and,  therefore,  the  dependence  of  the 
critical  current  density,  jc,  on  boundary  misorientation. 

In  the  framework  of  the  suggested  model,  the 
angular  dependence  of  jc  was  calculated  as  that  opti¬ 
mally  corresponding  to  the  experimental  data  [1],  in 
which  case  parameters,  a  and  A,  was  used  as  adjusting 
parameters.  In  doing  so,  for  characteristic  values  of  p 
=  0.2  [21],  a  =  10'2  and  EJb=  5,  we  find  the  optimal 
correspondence  between  our  theoretical  calculations  of 
y  ,(0)  and  the  experiments  [1]  (Fig.3)  to  come  into  play 
at  a  «  4.13  and  A  =0.3 13.  Fig.3  is  indicative  of  a 
satisfactory  agreement  between  the  model  under 
consideration  and  the  experimental  data  [1], 


Analysis  of  the  theoretical  dependence y](0)  shows 
that  this  dependence  is  weakly  affected  by  parameters, 
Ejb  and  ac,  because  a  ~  0.83 EJb,  and  Ao/))  runs  paral¬ 
lel  with  a.  at  constant  a.  When  Ijb  grows,  value  of  Aiipl 
weakly  increases.  With  this  taken  into  conside-ration, 
we  find  that  “structural  transition”  from  low-angle  tilt 
boundaries  (weakly  disordered  walls  of  dislocation  with 
large  Burgers  vectors;  see  Fig. lb)  to  high-angle  tilt 
boundaries  (containing  highly  disordered  ensembles 
of  dislocations  with  small  Burgers  vectors;  see  Fig.2b) 
occuring  at  0  «  0c  does  not  result  in  any  jump-like 
changes  of  the  theoretical  dependence  jc  (9)  at  0  »  0c 
(see  Fig.3). 

4.  CONCLUSION 

In  conclusion,  in  this  paper  we  have  developed  a  theo¬ 
retical  model  which  describes  the  tilt  boundary  effect 
on  high-7]  superconductivity  as  that  related  to  long- 
range  stresses  (which  suppress  the  superconducting 
order  parameter)  created  by  tilt  boundaries  with 
disorderdly  arranged  dislocations.  In  this  context,  the 
disorder  of  grain  boundary  dislocation  ensembles  is 
caused  by  highly  non-equilibrium  conditions  of  syn¬ 
thesis  of  polycrystalline  high-7]  superconductors.  The 
angular  dependence,  jc  (9),  calculated  within  the  frame¬ 
work  of  the  suggested  model,  is  in  a  satisfactory  agree¬ 
ment  with  the  experimental  data  [1]  (see  Fig.3). 
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Abstract.  After  a  short  review  of  recent  TEM  measurements  to  characterize  partial  disci  inations  in  fragment 
boundary  triple  junctions  in  cold-rolled  copper,  STM  results  from  the  same  specimens  are  presented  which 
support  a  disclinalion  interpretation  of  the  observed  fragment  structure.  The  role  of  the  fragment  boundary 
mosaic  or  of  the  corresponding  immobile  partial  disclination  network  in  work-hardening  is  discussed. 
Although  the  orientational  mismatches  around  the  fragment  boundary  triple  junctions  are  only  small,  the 
disclination  contribution  dominates  work-hardening  behaviour  at  large  strains. 


1.  INTRODUCTION 

The  substructure  development  of  crystalline  metals 
under  cold  deformation  up  to  large  strains  is 
characterized  by  the  formation  of  a  fragment  structure 
on  the  mesoscopic  scale  [1],  Typically,  with  increasing 
strain  the  mean  fragment  size  decreases  monotonously 
starting  from  the  micrometer  range,  whereas  the  mean 
misorientation  across  the  fragment  boundaries  in¬ 
creases  monotonously  reaching  several  degrees  [1-3]'. 
In  contrast  to  the  cell  structure,  the  fragment  structure 
does  not  saturate1 2  and  its  formation  and  evolution  are 
quite  universal  phenomena  which  are  observed  in  fee, 
hdp,  and  bcc  single-  as  well  as  polycrystals  [1]. 

While  the  coexistence  of  the  cell  and  fragment  struc¬ 
tures  and  the  evolutions  of  their  characteristic  param¬ 
eters  are  well  investigated  [1-3],  the  defect  composition 
of  the  fragment  boundary  mosaic  and  its  influence  on 
work-hardening  are  not  yet  completely  understood. 
Especially,  to  the  authors’  knowledge,  there  is  no  work¬ 
hardening  model  up  to  now  which  takes  the  fragment 
structure  explicitly  into  account. 

Following  the  Russian  school  of  plasticity,  the  au¬ 
thors  propose  the  triple  junctions  of  the  fragment 

1  Flanscn  and  coworkcrs  use  the  term  «ccll  b!ock»  instead  of 
fragment. 

2  Saturation  means  in  this  context  that  average  cell  structure 
parameters  like  cell  size,  cell  wall  width,  or  misorientation  do 
not  change  any  more. 

Corresponding  author:  M. 


boundaries  to  have  a  disclination  character  [4-6]  and, 
thus,  to  be  sources  of  long-range  stresses  hindering 
dislocation  motion  and,  thereby,  contributing  to  work¬ 
hardening  at  large  strains. 

2.  NONCOMPENSATED  NODES  OF 
FRAGMENT  BOUNDARIES 

According  to  experimental  observations,  a  model  de¬ 
scribing  the  fragment  structure  development  as  well 
as  its  coupling  to  the  macroscopic  mechanical 
behaviour  requires  tools  which  take  into  account  such 
collective  behaviour  of  dislocations  that  results  into  the 
formation,  spread-out  and  enlargement  of 
misorientations  and  into  the  presence  of  long-range 
stresses  in  the  fragment  interiors.  A  model  on  the  sub¬ 
structure  development  and  work-hardening  under  cold 
deformation  up  to  large  strains  suggested  by  the  present 
authors  [7,8]  uses  disclinations  [9,10]  as  such  a  tool. 
Especially,  the  formation  of  partial  disclination  dipoles 
(PDD)  is  modelled  through  incidental  non-balanced 
trapping  of  neighbouring  mobile  dislocations  of  the 
same  Burgers  vector  into  cell  walls  and  by  incidental 
trapping  of  mobile  dislocations  with  “matching” 
Burgers  vector  into  preexisting  fragment  boundaries 
[8].  The  first  PDD  configuration  corresponds  to  a  both- 
side  terminated  excess  dislocation  wall,  the  second  one 
to  a  tilt  or  twist  boundary  segment  with  enlarged 
misorientation.  To  describe  the  non-conservative  propa- 
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Fig.  1.  A  propagating  partial  disclination  approaches  a  preexisting  fragment  boundary,  is  locked  in  this  boundary  and  leaves  an 
orientational  mismatch  behind. 


gation  of  the  partial  disclinations  of  such  dipoles  along 
cell  walls  or  fragment  boundaries,  respectively,  the 
capturing  mechanism  proposed  by  Vladimirov  and 
Romanov  [10,1 1]  is  used. 

If  a  partial  disclination  propagates  like  that  and 
approaches  a  preexisting  fragment  boundary,  then  it 
will  rotate  the  crystal  regions  to  the  left  and  to  the  right 
of  the  growing  backward  fragment  boundary  against 
each  other  (Fig.  lb).  When  it  is  finally  locked  in  the 
preexisting  fragment  boundary,  an  orientational 
mismatch  is  left  behind  (Fig.  lc).  This  mismatch  can 
be  accomodated  in  two  different  ways:  either  “plasti¬ 
cally”  by  generation  of  an  additional  dislocation  bound¬ 
ary  or  “elastically”  by  rotational  distortion.  In  the  latter 
case,  the  triple  junction  is  a  non-compensated  node 
containing  a  line  defect  of  the  disclination  type. 

3.  TEM  RESULTS 

Since  partial  disclinations  locked  in  the  triple  junc¬ 
tions  of  the  fragment  boundaries  should  -  at  least  at 
low  homologous  temperatures  -  be  stable  under  load 
relaxation,  they  are  suitable  objects  for  an  experimen¬ 
tal  confirmation  of  the  disclination  approach  used  in 
the  above  mentioned  model.  TEM  microdiffraction 
measurements  around  fragment  boundary  triple  junc¬ 
tions  in  copper  single  and  polycrystals  rolled  down  to 
70%  and  50%  thickness  reduction  at  room  tempera¬ 
ture,  respectively,  showed  [12] 

-  that  the  product  of  the  three  misorientation  matri¬ 
ces  does  not  give  the  identity  and 

-  that  the  local  orientation  varies  continuously  when 
passing  with  the  electron  beam  by  the  node. 

Whereas  the  first  observation  proves  an  orientational 
mismatch,  the  second  one  points  towards  its 
accomodation  by  elastic  distortion  rather  than  by  an 
additional  dislocation  boundary.  Both  items  together 
indicate  that  the  triple  junction  includes  an  additional 
line  defect  which  has  -  probably  beside  other 


contributions  [13,14]  -  a  disclination  nature.  The 
measured  orientational  mismatches  correspond  to 
Frank  vectors  of 

a  =  [-0.0292,0.0094,0.0052]  with \m\=  1.78° ,  (1) 

for  a  node  in  the  single  crystal  specimen  and 

a  =  [-0.0046  -0.0316,0.0004]  with  \a\  =  1.93°,  (2) 

for  a  node  in  the  polycrystal  specimen.  These  PD  pow¬ 
ers  are  in  good  agreement  with  the  theoretical  expec¬ 
tation  of  |<5|«l-3°  given  by  Romanov  and  Vladimirov 
in  [10].  At  present,  PDD  groups  in  a  copper  single 
crystal  rolled  down  to  60%  thickness  reduction  as 
shown  in  the  TEM  micrographs  in  Figs.  2  and  3  are 
analysed. 

First  results  indicate  that  the  PDD  groups  are  not 
fully  compensating  dipole  (Fig.  2)  and  quadrupole  (Fig. 
3)  configurations,  respectively  [15]. 

4.  STM  RESULTS 

Since  it  is  very  difficult  to  find  sufficiently  large  per¬ 
fect  crystal  volumes  to  obtain  Kikuchi  patterns  in  the 
surroundings  of  the  nodes  in  heavily  deformed  metals, 
an  additional  local  characterization  method  which  does 
not  make  use  of  diffraction  is  desirable.  Scanning  tun¬ 
nel  microscopy  (STM)  is  such  a  method.  In  their  pio¬ 
neering  works,  Vettegren’  et  al.  [16-20]  have  studied 
the  profiles  of  polished  surfaces  of  loaded  copper,  gold, 
and  molybdenum  specimens  and  found  clusters  of  etch 
pits.  The  pits  have  been  explained  by  the  exit  of  dislo¬ 
cation  groups  after  bursting  through  barriers  on  inter¬ 
secting  glide  planes,  whereas  the  clusters  have  been 
interpreted  as  an  early  stage  of  crack  formation  [21]. 

Unfortunately,  the  STM  contrast  theory  is  not  yet 
very  well  developed.  It  is  known  [22]  that  the  strong 
electric  field  between  the  tip  and  the  surface  has  an 
etching  effect.  It  is  generally  admitted  that  this  “field 
etching”  should  also  result  into  preferred  evaporation 
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Fig.  2.  TEM  micrograph  of  a  partial  disclination  dipole 
configuration.  Copper  single  crystal  rolled  down  to  60% 
thickness  reduction  at  room  temperature.  2.61x1.90  pm2. 


Fig.  3.  TEM  micrograph  of  a  partial  disclination  quadrupolc 
configuration.  Copper  single  crystal  rolled  down  to  60% 
thickness  reduction  at  room  temperature.  1.59x1.16  pm2. 


and/or  surface  diffusion  at  high-energy  surface  regions. 
From  an  analogy  to  Mullins’  thermal  etching  theory 
[23],  one  can  conclude  that  etch  grooves  should  arise 
and  deepen  at  the  lines  where  planar  defects  like  grain 
boundaries  or  fragment  boundaries  reach  the  surface. 
Depth  and  shape  z(x,t)  of  the  grooves  scale  with  the 
plane  energy  y  of  the  planar  defect,  e.g.  for  thermal 
etching  with  dominating  evaporation  and  condensation 
through  [23] 


z(x,i)  =  -2  tan 


(3) 


with  p  defined  by  the  mechanical  equilibrium  between 
the  planar  defect  (e.g.  grain  boundary)  tension  and  the 
surface  tensions,  2y(/sin|3  =y>(),  and  A  given  by 


A  =  - 


p„r,,Ql 


(2  7tM)m(kJ)1 


(4) 


with  vapour  pressure  p0  for  plain  surface,  atomic  vol¬ 
ume  fi.  atomic  mass  M  and  temperature  T.  Fragment 
boundaries  can  be  approximated  as  low-angle  grain 
boundaries,  so  that  their  plane  energy  can  be  expressed 
through  the  misorientation,  e.g.  for  a  tilt  boundary  [24] 


Gb 


4^(1  -  v) 


In- - In  (p 

v  j 


(5) 


with  misorientation  cp  and  dislocation  core  radius  rQ. 
Thus,  at  least  in  the  thermal  etching  analogy,  the  groove 
depth  and  shape  are  directly  related  to  the  fragment 
boundary  misorientation.  Correspondingly,  “field  etch¬ 
ing”  should  also  result  into  notching  of  etch  pits  at  the 
points  where  line  defects  like  dislocations  or 
disci  inations  reach  the  surface.  Depth  and  shape  of  the 
pits  should  then  scale  with  the  line  energy,  that  means 


with  the  powers  or  with  the  moduli  of  the  Burgers  or 
Frank  vectors,  respectively.  Although  Romanov  has 
given  the  elastic  fields  around  PDDs  reaching  the  sur¬ 
face  [25],  there  is  no  etch  pit  theory  up  to  now.  not 
even  for  the  thermal  analogy.  Therefore,  the  following 
discussion  has  to  remain  a  qualitative  one:  If  the  triple 
junctions  of  fragment  boundaries  include  additional 
high-energy  line  defects,  then  the  corresponding  etch 
pits  should  be  much  deeper  than  the  etch  grooves 
corresponding  to  the  joining  fragment  boundaries. 

For  the  present  work,  cut  surfaces  of  copper  poly¬ 
crystals  rolled  at  room  temperature  have  been  abraded, 
polished  and  treated  with  the  basic  alumina  OPS  be- 


Fig.  4.  STM  surface  profile  with  fragment  boundary  grooves. 
Copper  polycrystal  rolled  down  to  70%  thickness  reduction  at 
room  temperature.  Arrow  “A"  marks  a  splitting  point  of  a 
fragment  boundary,  arrow  “B”  a  triple  junction,  arrows  “C”  a 
chain  of  nearly  equidistant  shallow  etch  pits.  5x5  pm2.  Grey 
scale  from  white  to  black  corresponding  to  15  nm  depth. 
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fore  being  investigated  with  a  Nanoscope  II  STM.  Fig. 
4  shows  a  5x5  pm2  surface  profile  from  a  copper  poly¬ 
crystal  rolled  down  to  71%  thickness  reduction.  The 
parallelity  of  the  two  grooves  on  the  left  (to  each  other 
and  also  to  the  very  much  different  structure  of  a  chain 
of  nearly  equidistant  small  pits  marked  with  the  arrows 
“C”!),  the  similarity  of  the  defect  structure  with  the 
one  observed  in  TEM  (confer  e.g.  [2,26]),  especially 
the  band  width  in  the  micrometer  range  and  the  split¬ 
ting  up  of  one  groove  into  two  neighbouring  ones  at 
the  point  “A”  (confer  splitting  up  of  dense  dislocation 
walls  (DDW)  into  first  generation  microbands  (MB1), 


e.g.  [2,26])  point  to  a  deformation-induced  rather  than 
to  a  polishing  effect. 

Fig.  5  gives  the  scanning  profiles  for  the  three  lines 
marked  in  Fig.  4,  line  1  cutting  through  the  two 
neighbouring  ordinary  grooves,  line  2,  in  addition  to 
this,  cutting  through  the  triple  junction  marked  with 
“B”,  and  line  3  cutting  through  the  splitting  point 
marked  with  “A”.  Obviously,  the  pits  at  the  triple  junc¬ 
tions  are  deeper  than  the  ordinary  grooves  by  about  a 
factor  2  which  might  indicate  the  presence  of  an  addi¬ 
tional  high-energy  line  defect,  e.g.  a  partial  disclination. 


Fig.  5.  STM  surface  profiles  along  the  lines  1-3  in  Fig.  4.  400  scaning  points  correspond  to  5  pm. 
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5.  EFFECT  ON  WORK-HARDENING 

In  a  disclination-based  model  for  the  fragment  struc¬ 
ture  development  like  the  one  proposed  in  [7,8],  propa¬ 
gating  PDDs  are  immobilized  as  discussed  in  section 
2  and  Fig.  1  -  resulting  in  a  mosaic  of  fragment  bound¬ 
aries  or,  correspondingly,  in  a  network  of  immobile 
partial  disclinations  in  the  node  lines  of  this  mosaic. 
The  evolutions  of  the  cell  wall  dislocation  density  and 
of  the  propagating  and  immobile  disclination  densi¬ 
ties  can  be  described  with  the  equations  [8] 


calculated  by  distributing  the  available  excess  disloca¬ 
tion  density  which  develops  according  to  (£-term  de¬ 
scribes  capturing  of  mobile  dislocations  by  propagat¬ 
ing  partial  disclinations) 

—  =Kt4pt+KfJe,+E0i,  (11) 

ds 

on  the  total  fragment  boundary  plane  area,  finally  giv¬ 
ing  [7] 

(p  *  bpca  /  \2\Jo.  (12) 


dp 

ds 
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including  immobilization  (/-term)  of  mobile  and  dy¬ 
namic  recovery  (/?-term)  of  mobile  with  cell  wall  dis¬ 
locations,  PDD  generation  by  incidental  non-balanced 
trapping  of  N  neighbouring  mobile  dislocations  into 
cell  walls  (//  -term)  and  by  incidental  trapping  of  mo¬ 
bile  dislocations  with  the  “right”  Burgers  vector  into 
fragment  boundaries  (//-term),  and  immobilization  of 
propagating  PDDs  in  the  network  of  immobile  ones 
(./-term).  The  coefficients  can  be  estimated  from  physi¬ 
cal  models  for  the  elementary  processes  [7,8].  Whereas 
the  cell  wall  dislocation  density  saturates  according  to 
the  Kocks-type  equation  (6),  the  immobile  disclination 
density  increases  monotonously  according  to  equation 
(8)  -  corresponding  to  a  monotonously  decreasing  mean 
fragment  size. 

The  cell  wall  dislocation  density  is  coupled  to  a 
flow  stress  contribution  according  to 


A<7 


r 


%aGb  Vp~t 

m  \  £ 


(9) 


where  ms.  is  the  Schmid  factor  and  denotes  the  vol¬ 
ume  fraction  of  the  cell  walls.  The  effect  of  the  immo¬ 
bile  disclination  density  on  the  flow  stress  is  described 
by  Romanov’s  and  Vladimirov’s  law  [10] 

Acrtf  =  pG\o)\.  (10) 

The  mean  disclination  power  |tb|  is  equivalent  to  the 
average  orientational  mismatch  around  the  node  lines 
and  may  -  following  the  discussion  in  section  2  and 
Fig.  1  -  be  approximated  by  the  mean  misorientation 
cp.  Assuming  a  Poisson-Voronoi  geometry  of  the  frag¬ 
ment  boundary  mosaic,  the  mean  misorientation  is 


Integration  of  the  evolution  equations  (6-8)  and  the 
flow  stress  contributions  (9-10)  shows  that  the 
disclination  contribution  to  flow  stress  dominates  the 
dislocation  one  starting  from  a  von  Mises  strain  of  about 
0.9  (parameters  for  copper  at  room  temperature).  Stage 
111  of  single  crystal  plastic  deformation  can  be  attrib¬ 
uted  to  cell  structure  or  dislocation  work-hardening, 
stage  IV  to  fragment  or  disclination  work-hardening 
[7,8]. 

6.  CONCLUSIONS 

TEM  microdiffraction  measurements  indicate  the  pres¬ 
ence  of  line  defects  of  disclination  nature  in  the  triple 
junctions  of  fragment  boundaries.  STM  shows  deep  etch 
pits  at  the  points  where  the  triple  junction  lines  reach 
the  surface  —  pointing  to  high-energy  line  defects  in 
the  triple  junctions.  Although  the  orientational 
mismatches  remain  small,  the  disclination  long-range 
stress  fields  do  have  a  significant  influence  on  work¬ 
hardening  and  even  dominate  the  dislocation 
contribution  at  large  strains. 
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